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Radiating  atoms  in  dense  plasmas  are  highly  affected  by  their  local  plasma 
environment.  We  examine  the  effect  of  this  environment  on  hydrogenic  radi- 
ators by  using  an  atomic  wave  function  basis  set  that  is  itself  a  function  of 
the  plasma  electric  microfield  for  the  calculation  of  spectral  line  profiles.  Our 
theoretical  development  includes  the  effect  of  static  ion  perturbers  up  to  the 
quadrupole  term  in  the  multipole  expansion  of  the  radiator-perturbing  ion  in- 
teraction, as  well  as  dynamic  electron  broadening  up  to  second  order  in  the 
radiator-perturbing  electron  interaction.  Effects  due  to  the  presence  of  the  uni- 
form ion-electric  field  are  treated  exactly  within  this  framework,  and  the  field 
gradient  is  then  treated  as  a  perturbation  by  inclusion  of  the  ion- quadrupole 
term.  We  employ  a  basis  set  for  the  representation  of  the  upper  state  of  the 
transition  that  includes  states  of  principal  quantum  number  n  as  well  as  n+1 
to  allow  for  the  field  mixing  of  these  two  manifolds.  The  ion-quadrupole  effect 
is  treated  in  a  new  way  that  includes  ion  correlations  by  using  a  renormalized 

vi 


independent-particle  model  for  the  perturbing  ions.  We  also  present  a  pre- 
liminary assessment  of  the  importance  of  field  ionization  on  the  spectral  line 
shapes.  Our  results  are  compared  with  previously  reported  calculations  using 
field  independent  matrix  elements,  neglect  of  possible  field  ionization  and  the 
independent-particle  model  of  the  ion-quadrupole  effect.  Hence,  we  are  able  to 
identify  the  most  important  effects  that  arise  from  the  use  of  the  field  depen- 
dent basis  set.  These  comparisons  also  allow  us  to  examine  the  importance  of 
ion  correlations  when  determining  the  ion-quadrupole  effect. 
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CHAPTER  I 
INTRODUCTION 

In  this  dissertation  we,  that  is  you  and  I,  will  be  examining  the  theory  of 

spectral  line  radiation  from  highly  ionized  atoms  in  dense  plasmas.  Before  we 

go  into  this  in  detail  it  will  be  helpful  and  illuminating  to  stand  back  and  look 

at  the  relevance  of  this  from  a  wider  perspective;  to  take  a  look  at  the  big 

picture. 

1.1  The  Big  Picture 
In  general,  a  plasma  can  be  defined  as  a  statistical  system  containing  mo- 
bile charges.  We  will  be  concerned  with  a  subcategory  of  this  definition  and 
will  limit  ourselves  to  discussing  electrically  neutral  plasmas  associated  with 
hot  highly  ionized  atoms.  This  is  no  great  limitation  since  over  99  percent 
of  the  known  universe  is  covered  in  this  category.2  With  the  aid  of  sophisti- 
cated astronomical  detection  equipment,  most  of  what  we  see  of  the  universe 
is  electromagnetic  radiation  emanating  from  a  variety  of  objects  ranging  from 
ordinary  stars  to  extremely  luminescent  quasistellar  objects  located  at  the  edge 
of  the  known  universe.  The  matter  in  all  of  these  entities  exists  almost  entirely 
in  the  plasma  state;  plasmas  whose  constituents  range  from  ordinary  ionized 
atoms  to  exotic  matter-antimatter  electron  pairs.  The  radiation  reaching  the 
earth  from  these  distant  bodies  runs  the  gamut  from  radio  waves,  through  the 
visible  spectrum  and  into  the  realm  of  X-rays  and  7-rays.  It  is  easy  to  see  that 
an  understanding  of  the  behavior  of  these  plasmas  can  and  will  give  us  great 
insight  into  the  structure  and  workings  of  the  universe  in  which  we  live. 
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Of  more  mundane  concerns  are  the  plasmas  encountered  here  on  earth. 
These  are  evident  in  such  natural  phenomena  as  lightning  and  fire.  Man- 
made  plasmas  are  also  numerous  in  example,  ranging  from  ordinary  fluorescent 
lighting  to  terrifying  large  scale  thermonuclear  explosions  that  result  from  the 
fusion  of  hydrogen  isotopes.  Since  the  early  1950s,  one  of  the  primary  goals  of 
plasma  physics  has  been  to  bring  under  control  these  tremendous  quantities  of 
energy  liberated  by  thermonuclear  fusion.  When  this  goal  is  achieved,  we  will 
have  succeeded  in  obtaining  access  to  a  nearly  unlimited  supply  of  relatively 
clean  and  safe  energy  that  will  last  into  the  foreseeable  future.** 

Due  to  the  extremes  of  physical  conditions  under  which  many  plasmas  ex- 
ist and  due  to  their  often  ephemeral  nature,  making  direct  measurements  of 
the  composition,  density,  temperature  and  other  physical  parameters  that  are 
necessary  for  a  clear  understanding  of  their  inner  workings  is  often  impossi- 
ble. What  is  needed  is  an  indirect  nonintrusive  probe  of  the  plasma  interior. 
The  radiation  emitted  or  absorbed  by  these  plasmas  is  the  ideal  answer  to  this 
problem.  Isolated,  partially  ionized  atoms  can  emit  or  absorb  radiation  that 
results  from  atomic  electrons  passing  from  one  energy  level  to  another.  This 
radiation  is  emitted  or  absorbed  in  discrete  amounts  leading  to  distinct  and 
identifiable  sharp  spectral  lines.  These  lines  have  a  slight  width,  called  the 
natural  width,  due  to  quantum  fluctuations  in  the  atom's  own  electromagnetic 
field.  As  long  as  the  atoms  remain  isolated,  the  spectral  lines  remain  rela- 
tively sharp.  When  linked  to  their  plasma  environment,  however,  the  situation 
changes  radically.  The  amounts  of  energy  emitted  or  absorbed  by  the  atoms 
are  now  affected  by  the  plasma,  often  revealing  pertinent  information  about 
this  environment.  Thus  we  may  explore  the  atom's  surroundings  indirectly 
through  changes  in  the  profile  of  its  spectral  lines.  This  technique  is  employed 


in  many  areas  of  science  to  reveal  information  about  the  internal  structure 
of  matter  in  a  noninvasive  way  by  examining  emitted  or  absorbed  radiation. 
Some  examples  of  this  application  are  nuclear  magnetic  resonance  (NMR)  and 
Mossbauer  spectroscopy.  In  this  dissertation  we  will  be  concerned  with  spec- 
tral lines  of  emitted  radiation  from  plasmas  although  the  formalism  applies 
equally  well  to  spectral  lines  originating  from  the  absorption  of  radiation.  The 
line  shapes  will  reveal  information  about  the  density  and  temperature  in  the 
vicinity  of  the  radiating  atoms,  regardless  of  whether  they  are  located  in  a 
simple  laboratory  experiment  or  in  a  distant  star. 


1.2  Focus  of  This  Work 
The  main  focus  of  this  work  will  be  on  the  shapes  of  line  spectra  from 
dense,  high  temperature  plasmas  consisting  of  atoms  of  moderate  atomic  num- 
ber, Z.  These  terms  "dense"  and  "high"  are  quite  relative  and  historically  have 
meant  the  densest  and  hottest  plasmas  obtainable  in  laboratory  experiments 
for  a  particular  element  at  that  particular  point  in  time.  As  an  example,  in 
the  mid  1960's  these  terms  referred  to  electron  number  densities,  ne,  of  around 
10  '  cm-"*  and  temperatures  corresponding  to  a  few  eV  for  elements  such  as 
hydrogen  and  helium.5  Today,  high  density  refers  to  ne  =  1024  to  1025  cm-3 
and  temperatures  corresponding  to  600  to  1000  eV  for  elements  such  as  argon 
and  krypton  with  Z  =  18  and  36  respectively.6,7 >8>9>10  One  method  of  creat- 
ing these  dense  plasma  conditions  in  the  laboratory  is  by  laser  driven  inertial 
confinement.  Current  experiments  using  this  method  have  created  conditions 
of  many  times  solid  density  for  periods  of  a  few  hundred  picoseconds.9  The 
mechanism  of  compression  employs  high  intensity  laser  light  to  symmetrically 
irradiate  spherical  targets  containing  the  material  under  investigation.  *  When 


this  energy  is  rapidly  deposited  in  the  outer  layers  of  the  target  (mostly  by  in- 
verse Bremsstrahlung),  these  layers  are  quickly  vaporized  and  blown  outwards 
away  from  the  target.  This  propels  the  remainder  of  the  target  material  in- 
ward. This  is  in  close  analogy  to  the  way  a  rocket  works  by  combusting  fuel  in 
an  engine  and  expelling  it  in  one  direction  in  order  to  be  propelled  in  the  other 
direction.  If  our  laser  irradiation  is  sufficiently  symmetrical,  the  target  will 
be  reduced  in  size  many  times  as  it  is  compressed,  and  resulting  shock  waves 
will  heat  the  core  to  a  high  temperature.  This  compression  and  heating  can 
produce  the  plasmas  of  high  density  and  temperature  that  we  wish  to  examine 
in  this  dissertation. 

Approximate  descriptions  of  spectral  line  shapes  emitted  from  dense,  hot 
plasmas  are  usually  derived  from  theories  that  use  an  atomic  physics  descrip- 
tion of  atoms  that  are  isolated  from  their  environment.  In  reality,  the  atomic 
wave  functions  for  the  radiating  atoms  can  be  strongly  perturbed  by  surround- 
ing charged  particles,  and  in  extreme  cases  it  is  not  realistic  to  think  of  the 
atoms  as  separate  entities  at  all  but  rather  as  clusters  of  ions  surrounded  by 
clouds  of  electrons.  Radiation  from  atoms  in  this  highly  interacting  environ- 
ment will  run  the  gamut  from  continuum  to  discrete.  Therefore,  in  order  to 
construct  a  tractable  method  for  the  calculation  of  spectral  line  shapes,  we 
must  often,  somewhat  arbitrarily,  distinguish  between  the  radiator  and  its  per- 
turbers.  However,  the  more  influence  of  the  surrounding  particles  that  we  can 
include  in  the  description  of  our  radiating  atoms  at  the  outset  of  our  develop- 
ment, the  better  off  we  will  be  when  considering  further  perturbations.  Our 
goal  in  this  work  will  be  to  include  some  of  the  higher  order  effects  of  the 
plasma  microfield  on  atomic  spectral  line  emission.  To  this  end,  we  will  use 
for  our  zero-order  atomic  wave  functions,  numerical  solutions  to  Schrodinger's 


equation  for  an  atom  in  a  uniform  electric  field.  These  field  dependent  solutions 
are  supplemented  by  solutions  obtained  from  a  systematic  perturbation  theory 
treatment  for  small  field  values.  The  zero  order  atomic  Hamiltonian  will  thus 
contain  the  perturbing  ion-dipole  interaction  and  higher  order  corrections  will 
describe  nonuniformities  in  the  plasma  electric  field.  We  will  include  the  first 
of  these  higher  order  corrections,  the  ion-quadrupole  term,  to  account  for  the 
effects  of  the  field  gradient.12'13'14,15  Griem16  has  pointed  out  as  early  as  1954 
the  possible  importance  of  this  field  dependence,  as  well  as  the  ion  quadrupole 
effect,  as  a  source  of  spectral  line  asymmetry. 

The  plasma  electric  microfield  can  produce  additional  effects  on  the  radi- 
ating atoms  beyond  those  given  by  the  simple  field  dependence  of  their  atomic 
wave  functions  and  energy  levels.  Atomic  bound  state  energies  which  were  dis- 
crete in  the  absence  of  the  plasma  electric  field  now  become  resonance  states 
with  finite  widths.  The  formerly  bound  electrons  are  now  free  to  quantum 
mechanically  tunnel  through  the  potential  barrier  created  by  the  electric  field. 
Although  atomic  electrons  that  lie  deep  in  the  radiator  potential  well  are  little 
affected  by  the  plasma  electric  field,  high  lying  electron  levels  are  perturbed 
to  the  point  where  they  are,  for  all  practical  purposes,  no  longer  bound  to  the 
atom.  This  field  ionization  phenomenon  could  have  important  consequences 
for  the  relative  populations  of  the  radiator  excited  states. 

1.3  Outline 
Chapter  2  focuses  on  the  development  of  a  general  plasma  line  broaden- 
ing theory  and  the  accompanying  approximations.    We  discuss  fundamental 
length  and  time  scales  for  the  interactions  of  the  two  main  components  of  the 
plasma  with  the  radiator:    the  highly  charged  ions  and  the  unbound  plasma 
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electrons.  We  also  discuss  under  what  conditions  and  where  in  the  line  profile 
the  phenomenon  of  perturber  motion  will  be  important  (both  ion  and  elec- 
tron). Next  follows  a  discussion  of  the  strength  and  nature  of  the  various 
interactions  between  the  plasma  particles.  From  this  consideration,  we  decide 
which  interactions  are  important  and  which  can  be  ignored. 

A  brief  account  of  the  historical  development  of  the  main  ideas  of  spectral 
line  broadening  theory  is  presented.  We  find  the  historical  roots  for  the  ideas 
of  Stark,  Doppler  and  electron  broadening  in  the  late  1800's  and  early  1900's. 

Next  follows  a  detailed  derivation  of  a  calculable  form  of  the  line  shape 
function.  We  first  present  a  general  formalism  appropriate  for  multi-electron 
radiators  that  includes  field  dependent  wave  functions,  an  exact  treatment  of 
the  perturbing  ion  radiator-dipole  interaction,  field  dependent  NLTE  level  pop- 
ulations and  field  dependent  electron  broadening.  This  general  formalism  is 
then  restricted  to  the  case  of  hydrogenic  radiators  with  no  lower  state  broad- 
ening. We  also  discuss  a  method  for  approximately  including  field  induced 
resonance  widths  in  the  line  shape.  Chapter  2  ends  with  a  brief  note  about 
other  line  broadening  phenomenon  that  may  need  to  be  included  before  com- 
parisons with  experiment  can  be  carried  out. 

Chapter  3  begins  with  the  development  of  a  theory  for  the  approximate 
evaluation  of  the  ion-quadrupole  effect.  We  present  approximations  for  the 
probability  function  for  the  field  gradient  term  as  well  as  an  approximation  for 
the  average  field  gradient  constrained  to  have  a  given  field  value.  These  ap- 
proximations include  ion-ion  correlations  and  represent  an  improvement  over 
previous  theories.  We  then  compare  our  calculation  of  this  field  gradient  with 
previous  calculations  and  simple  approximations.  Next,  we  discuss  the  methods 


of  calculation  of  field  dependent  wave  functions  and  matrix  elements  by  pertur- 
bation theory  and  by  direct  numerical  solution  of  Schrodinger's  equation.  The 
numerical  solution  technique  is  also  used  to  calculate  resonance  widths  that 
give  an  approximate  measure  of  the  lifetimes  of  the  field  dependent  resonance 
states  produced  by  the  presence  of  the  plasma  microfield.  These  lifetimes  are 
then  used  to  study  the  effect  of  the  resulting  possible  level  depletion  on  the 
line  shape.  We  discuss  a  simple  population  kinetics  model  for  this  purpose. 

In  chapter  4,  we  present  some  theoretical  spectral  line  shapes  resulting 
from  our  work.  We  calculate  line  shapes  for  the  La  and  Lp  lines  of  hydrogenic 
argon  at  kT  =  800  eV  and  for  an  electron  density  range  of  1  x  1024  to  1  x 
1025cm-^.  We  then  examine  the  importance  of  the  various  effects  studied  in 
this  dissertation  on  the  spectral  line  shape,  width  and  asymmetry.  We  close 
with  some  conclusions  and  suggestions  for  further  work. 

Appendices  discuss  several  useful  topics  and  give  some  details  of  calcula- 
tions presented  in  the  main  text.  We  have  a  discussion  of  physical  units,  some 
details  of  the  model  used  to  calculate  the  average  field  gradient  required  for  the 
treatment  of  the  ion-quadrupole  effect,  a  discussion  of  parabolic  coordinates 
and  the  details  of  the  perturbation  theory  calculation  of  the  field  dependent 
fine  structure  corrections  to  the  radiator  Hamiltonian. 


CHAPTER  II 

THEORETICAL  BACKGROUND  TO  THE 

LINE  BROADENING  PROBLEM 

We  will  now  develop  a  theory  ''°  for  the  description  of  spectral  line  shapes 
in  dense  plasmas.  The  physical  conditions  we  will  be  interested  in  considering 
will  be  pertinent  to  inertia!  confinement  fusion  plasmas.  These  plasmas  will 
be  at,  or  near  solid  density  and  at  relatively  hot  temperatures  (i.e.  hot  enough 
that  the  plasma  is  not  quantum  mechanically  degenerate).  We  will  also  be 
interested  in  systems  where  the  plasma  ions  have  much  greater  mass  than 
the  electrons.  This  difference  will  allow  us  to  consider  different  domains:  one 
defined  by  the  fast  moving  or  dynamic  electrons  and  one  appropriate  for  the 
description  of  slow  moving  or  quasi-static  ions.  This  will  allow  us  to  proceed 
with  the  development  of  a  simplified  formalism. 

The  physical  systems  and  conditions  that  we  will  examine  in  this  disser- 
tation are  restricted  to  argon  radiators  immersed  in  a  pure  argon  plasma  at 
a  temperature  corresponding  to  800  eV  and  an  electron  number  density  rang- 
ing from  10  cm-**  to  10  ^  cm  .  These  densities  are  high  enough  to  bring 
out  the  higher  order  field  effects  but  not  high  enough  for  the  plasma  to  be  in 
the  degenerate  electron  regime.  In  local  thermodynamic  equilibrium  (LTE), 
the  line  shape  is  not  very  sensitive  to  small  variations  in  temperature.  Con- 
sequently, we  will  not  examine  the  temperature  dependence  of  the  line  shape. 
These  physical  conditions  are  also  relevant  to  current  experiments9  that  employ 
argon  radiators. 
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2.1  Dense  Plasma  Fundamentals 

In  dense  plasmas  there  are  several  fundamental  length  and  time  scales  that 
are  important  for  understanding  the  interactions  among  the  plasma  particles. 
Since  our  plasmas  are  dynamic  systems,  the  particles  are  in  constant  motion 
and  travel,  on  the  average,  a  characteristic  distance  in  a  given  relevant  time 
interval.  This  fact  allows  us  to  establish  some  relations  between  length  and 
time  scales.  The  use  of  these  scales  will  allow  us  to  classify  the  various  plasma 
particles  into  categories  that  are  governed  by  interactions  of  differing  strengths. 
This  will  be  important  when  deciding  what  approximations  are  suitable  and 
consistent  in  our  description  of  the  plasma.  Without  approximations  we  would 
be  hopelessly  lost  in  a  morass  of  ~  10"  coupled  equations  of  motion.  A  fate, 
preferably,  to  be  avoided. 

We  will  now  discuss  some  relevant  length  scales1'19'^  (see  Appendix  A). 
The  ion-sphere  radius  r0  is  an  estimate  of  the  average  distance  between  two 
ions  of  species  j .  It  is  given  by 


o,j 
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(2.1) 


where  ny  is  the  number  density  of  species  j. 

Another  important  quantity  is  the  Debye  length  \q.  It  is  an  indication 
of  the  distance  from  an  ion  beyond  which  screening  of  its  charge  becomes 
significant.  It  is  given  by 


(2.2) 


_  (       kT       xl/2 

where  k  is  Boltzmann's  constant,  T  is  the  temperature  and  Ze  is  the  charge  of 
species  j  that  forms  the  screening  cloud.  For  more  than  one  component,  the 


10 
resultant  screening  length  is  \Zr  =  £\-  A  «  .  where  the  sum  is  over  the  different 
species. 

The  thermal  de  Broglie  wavelength  A  gives  an  estimate  of  the  quantum 
mechanical  wave  nature  of  the  plasma  particles.  It  is  given  by 
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(2.3) 


where  m  is  the  particle  mass.  Except  for  the  factor  of  2k  this  is  just  the 
quantum  mechanical  wave  length  of  plasma  particles  possessing  average  kinetic 
energy  kT.  It  gives  us  an  estimate  of  the  spatial  extent  of  the  particle  wave 
packet. 

The  radial  dimension  of  an  atomic  ion  can  be  estimated  by  considering  the 
radius  of  the  Bohr  orbit  for  the  most  probable  quantum  state,  and  writing  the 
radius  as. 


rn  =  a0n  jZ  ; 


(2.4) 


n  is  the  principal  quantum  number  for  the  atomic  ion,  a0  is  the  radius  of  the 
first  Bohr  orbit  for  hydrogen  and  Z  is  the  nuclear  charge. 

In  order  to  relate  time  scales  to  characteristic  lengths  we  need  an  estimate 
of  the  velocity  for  particles  of  species  j.  Using  the  magnitude  of  the  most 
probable  velocity  for  classical  motion,  we  have 


1/2 


vmP,j  =  \2kT/mjj 


(2.5) 


We  have  now  defined  the  basic  quantities  which  will  allow  us  to  make 
estimates  of  several  important  time  scales  of  interest.  Additionally,  the  length 
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scales  will  be  used  in  determining  appropriate  scale  factors  for  perturbation 
expansions. 

We  will  now  look  at  the  time  scales  relevant  to  motion  of  the  plasma  parti- 
cles, specifically  the  perturbing  electrons  and  ions.  It  will  be  shown  in  section 
2.3  that  the  line  shape  function  /(Aw)  can  be  represented  by  the  transform  of 
a  dipole  autocorrelation  function  C(t )  such  that 

I(Au)=fexp(iAut)C(t)dt,  (2.6) 

where  Aw  is  the  separation  from  line  center  and  C(t)  is  a  decreasing  function 
of  t.  When  t  >  1/Aw,  the  exponential  begins  to  oscillate  rapidly  causing  the 
contribution  to  the  integral  from  the  integrand  in  that  region  to  be  small.  So, 
for  a  given  value  of  Aw,  the  integral  is  determined,  for  the  most  part,  by  C(t) 
such  that  0  <  t  <  1/Aw.  Hence,  we  can  define  the  time  interval  of  interest 
for  the  line  shape  at  the  distance  from  line  center  Aw  as  r  =  1/Aw.  If  the 
duration  of  a  perturber's  collision  with  the  radiator  is  significantly  greater  than 
the  time  of  interest  r  corresponding  to  the  part  of  the  spectral  line  Aw,  we  may 
regard  the  perturber  as  being  stationary  during  the  time  of  interest  and  treat 
its  perturbation  as  static.  This  is  known  as  the  quasi-static  approximation. 
If  the  duration  of  collision  does  not  meet  this  criterion,  we  must  regard  its 
perturbation  as  dynamic  and  explicitly  take  its  time  dependence  into  account. 
For  a  particular  separation  from  line  center  Aw,  the  effect  of  perturbing 
ion  motion  will  be  negligibly  small  if17 


F(t) 


F(t) 


«  \Auitf(F)l  (2.7) 

where  F(t)  is  the  perturbing  ion  electric  field  strength  at  the  radiator,  F(t)  is  its 
time  derivative  and  Awj  t(F)  represents  the  separation  from  line  center,  Aw, 


12 
due  to  the  static  linear  Stark  effect.  Note  that  Aujjt(F)  =  u^  t(F)  —  u?j  y(0), 
where  i  and  /  refer  to  the  initial  and  final  state  of  the  transition,  respectively, 
and  ujjf  is  the  transition  energy  u^  j  =  Ui—ut.  The  linear  Stark  shifted  energy 
for  level  i  is  u>;(0)  +  ^niqiUF /me.  We  use  the  parabolic  representation  (see 
Appendix  C)  and  take  F  =  Zje/r  where  r  is  the  distance  from  the  radiator 
to  the  perturbing  ion,  and  obtain  on  substitution  into  Eq.  (2.7) 


Alo  = 


_  3hZi(riiqi  -  rifqf)         v 
2Zmer2  >  r ' 


(2.8) 


where  n^  t  and  q^  r  are  the  parabolic  quantum  numbers  for  the  initial  and  final 
states  and  v  is  the  perturbing  ion  velocity.  Eliminating  r  from  both  sides  of 
the  inequality  gives 


Aw> 


2mev2Z 


3hZi(niqi  -  rifqf)  ' 
We  can  use  the  most  probable  velocity  as  an  estimate  for  v  to  obtain 


(2.9) 


Au;> 


AmekTZ 


=  AuJi  , 


(2.10) 


fihZi(niqi  -  rifqf) 

where  we  have  used  the  perturbing  ion-radiator  reduced  mass  \i  in  the  expres- 
sion for  v  because  we  are  only  interested  in  the  relative  motion.  Here,  we  have 
defined  Au;j  as  the  characteristic  shift  from  line  center  due  to  ion  motion.  The 
criterion  for  the  quasistatic  ion  approximation  to  be  valid  is  thus:  Au  ^>  Aw,'. 
For  the  Lyman  a  line  in  pure  argon  plasma  with  kT  =  800  eV,  qi  =  1  and 
qf  =  0  we  have  hAujj  =  0.044  eV.  For  higher  series  members  UAuJi  will  be  less. 
We  conclude  that  estimated  shifts  due  to  ion  dynamics  will  be  much  smaller 
than  the  estimated  shifts  due  to  the  static-ion  Stark  effect.  The  latter  are  on 
the  order  of  tens  to  hundreds  of  electron  volts  for  the  physical  conditions  we 
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are  interested  in  in  this  work.  Consequently,  we  can  confidently  use  the  qua- 
sistatic  approximation  for  the  perturbing  ions.   For  low  mass  perturbing  ions 
the  situation  is  quite  different.  In  the  case  of  proton  perturbers  ^Ao;j  =  16  eV 
and  ion  motion  will  become  a  significant  component  of  the  line  width. 

The  story  for  the  perturbing  electrons  is,  however,  another  matter.  Due 
to  their  much  lighter  mass,  the  electrons  will  be  traveling  much  faster  when 
they  have  the  same  kinetic  energy;  their  dynamic  influence  on  the  radiator 
will  need  to  be  considered.  If  we  treat  the  plasma  ions  as  static,  the  plasma 
electrons  will  form  shielding  charge  clouds  around  them,  screening  the  ionic 
charge  at  a  distance  of  roughly  the  Debye  length.20  This  is  given  by  Xp  e  = 
(kT/4ire2ne)1'2.  If  we  regard  the  time  of  interest  for  the  interaction  of  a 
perturbing  electron  with  the  screened  radiator  as  roughly  the  time  it  takes  to 
cross  the  Debye  sphere  radius  Xp  e,  we  have 


ve 


D,e 


vmp,e 


(2.11) 


where  we  have  estimated  ve  by  its  most  probable  value.  This  gives  us 


t71  =  V2lu 


p,e  , 


(2.12) 


where  u>Pje  is  the  electron  plasma  frequency  given  by 


/47rnee2\  V2 

P>e  =  (  j 

V     nip     J 


(2.13) 


The  portion  of  the  line  profile  that  corresponds  to  electron  dynamics  is 
then  given  by 


Au;  <  uPte  , 


(2.14) 


14 
with  a  quasistatic  electron  approximation  being  applicable  for  Au;  ^>  u}p>e.  For 
our  physical  conditions  up<e  is  in  the  neighborhood  of  10— 150  eV.  Examination 
of  the  statics  and  dynamics  of  the  perturbing  ions  and  electrons  allows  us  to 
conclude  that  the  quasistatic  approximation  for  ions  is  good  throughout  the  line 
profile  except  at  the  very  center  for  Aw  less  than  1  eV  while  electrons  behave 
dynamically  inside  u>p,e-  This  will  allow  us  to  separate  the  problem  of  plasma 
perturbations  of  the  radiator  into  two  regimes  coupled  only  by  the  static  plasma 
microfield:  that  of  purely  static  screened  ions  and  purely  dynamic  electrons. 
This  will  greatly  simplify  the  line  shape  theory  that  we  wish  to  formulate. 

For  our  plasma  particles  to  behave  non-degenerately  we  require  that  the 
interparticle  spacing  be  much  greater  than  the  quantum  mechanical  wavelength 
of  the  moving  particles.  For  the  plasma  electrons  this  condition  is 


Ae  <  r0,e  • 


(2.15) 


This  requirement,  except  for  a  slight  difference  in  the  numerical  constant,  is 
the  same  as  Ep/kT  <C  1,  where  the  electron  Fermi  energy,  Ep,  is  given  by 


^ 


2me 


For  our  plasma  at  kTe  —  800  eV  and  ne  —  1  x  1025  cm-3  the  above  ratio  is 
0.21.  However,  for  a  density  of  ne  =  1  x  10  "cm-3  the  ratio  becomes  0.98 
and  it  is  no  longer  a  reasonable  approximation  to  treat  the  plasma  electrons 
non-degenerately.  We  will  take  an  electron  density  upper  bound  of  ne  =  1  x 
10  5cm_<J  to  avoid  the  problem  of  electron  degeneracy  in  this  work.  The 
much  heavier  plasma  ions,  on  the  other  hand,  will  behave  classically  under  any 
conditions  we  will  be  considering  here. 
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To  judge  the  importance  of  correlations  between  particles  in  our  system  we 
can  look  at  the  Coulomb  coupling  constant  T.  This  is  the  ratio  of  a  particle's 
potential  energy  of  interaction  with  a  neighboring  particle  to  its  kinetic  energy. 
If  r  >  1  then  interactions  with  other  particles  will  dominate  ideal  gas  behavior. 
If  T  <  1  then  the  kinetic  energy  will  be  most  important  in  describing  the 
particle  behavior. 

First  we  look  at  the  coupling  between  electrons.  For  the  potential  energy  we 
use  the  potential  produced  by  an  electron  at  the  average  interparticle  spacing 
r0)e.  For  the  kinetic  energy  we  use  kT  to  give 


For  our  highest  density,  ne  =  1  x  1025cm-3,  this  gives  re,e  =  0.063.  We  con- 
clude that  the  plasma  electrons  are  weakly  coupled  to  each  other  and  can  be 
treated,  in  a  first  approximation,  as  an  ideal  Boltzmann  electron  gas.  Correla- 
tions between  electrons  will  be  treated  in  an  approximate  way. 
For  the  ions  we  use  r0  j  and  obtain 


*  i,i  — 


Tn  iK  1 


(2.18) 


but  since  r0j  =  Z^'^r0)e  we  have  T,-  j  =  Z5/3re,e-  For  our  fully  stripped  argon 
plasma,  Z=18  and  for  our  highest  density  we  have  T,-  ,•  =  7.8.  Clearly,  then,  we 
must  not  ignore  ion-ion  correlations.  In  particular  the  calculation  of  the  ion 
microfield  will  require  a  careful  treatment  of  these  correlations. 

For  ion-electron  correlations  we  take  as  the  interparticle  spacing  r0   = 
(r0,e  +  r0ji)/2  to  give 


r  ,-*l 

e'1       r0kT 


(2.19) 
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Noting  that  r0  =  (Z1/3  +  l)r0)e/2,  we  can  write 


2Z 


(2.20) 


For  ne  =  1  x  1025 cm-3  this  gives  re?-  =  0.62.  So  ion-electron  coupling  is  of 
order  one  and,  although  not  as  strong  as  ion-ion  correlations,  will  need  to  be 
taken  into  account.  This  will  be  accomplished  by  a  mean  field  approximation 
using  screened  ion  potentials. 

In  order  to  examine  the  effect  of  ion  radiator  interactions  it  will  be  use- 
ful to  make  a  multipole  expansion  of  the  radiator-perturbing  ion  interaction 
potential.  To  do  this  we  define  an  expansion  parameter  8,  where 


ra 


0,1 


(2.21) 


Here,  ra  =  (?/,n|r|V,n)  ~  a0vr / Z  is  the  most  probable  orbital  radius  for  a 
radiator  in  the  state  with  principal  quantum  number  n  and  nuclear  charge  Z. 
So 


This  is  8  =  8.531  x  10_9n2ne1/3/^4/3  with  ng  given  in  cm-3  At  ne  =  1  x 
1025  cm-3  for  n  =  2  we  have  6  =  0.16  and  for  n  =  3  we  have  6  =  0.35.  The 
dipole  term  of  the  multipole  expansion  is  second  order  in  8  and  the  quadrupole 
term  is  third  order  in  8.  For  n  =  3  this  gives  £2  =  0.12  and  <S3  =  0.043.  So  our 
expansion  is  quite  legitimate. 
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2.2  Historical  Background 

We  will  now  take  a  brief  look  at  the  historical  development  of  spectral  line 
shape  theory  with  the  goal  of  identifying  the  origin  of  the  main  ideas  of  line 
broadening. 

In  1889,  Lord  Rayleigh2"*  correctly  explained  the  broadening  of  spectral 
lines  due  to  the  Doppler  effect.  He  used  the  not  completely  accepted  idea  at 
that  time,  of  a  collection  of  randomly  moving  atomic  radiators  whose  velocity 
was  governed  by  Maxwell's  distribution.  Interest  was  such  that  by  1895  A.  A. 
Michelson24  was  able  to  summarize  the  possible  mechanisms  of  line  broadening 
as  he  saw  them.  Among  the  proposed  mechanisms  were  the  following: 

1)  The  change  in  the  emitted  radiator  frequency  due  to  neighboring  molecules. 
In  those  days  the  term  "molecule"  was  used  to  refer  to  any  atomic  sized 
entity;22 

2)  Doppler  broadening  due  to  radiator  motion; 

3)  Broadening  of  radiator  emission  lines  due  to  radiation  interrupting  colli- 
sions. Michelson  developed  a  theory  for  this  effect  and  showed  that  the 
resulting  line  shape  was  that  of  a  Lorentzian.  These  early  theories  consid- 
ered the  radiators  to  be  classical  oscillators.  We  will  see  that  collisional 
broadening  is  often  proportional  to  the  density,  and  hence,  to  the  pressure, 
assuming  the  temperature  is  held  constant.  This  is  where  the  terminology 
pressure  broadening  originated; 

4)  The  natural  line  width  due  to  radiation  damping  which  was  viewed  as  the 
leaking  away  of  the  oscillator's  energy  to  the  environment. 

In  1901  C.  Godfrey  5  combined  Doppler  and  collisional  broadening  in  a 
unified  description.  At  that  time  these  two  effects  were  viewed  as  the  most  sig- 
nificant causes  of  observed  line  broadening.  With  the  advent  of  the  quantum 


18 
theory  the  time  was  ripe  for  a  reworking  of  the  old  classical  line  broadening 
theory.  In  1924,  W.  Lenz26  applied  some  elements  of  the  new  quantum  theory 
to  collisional  broadening  but  still  considered  the  radiator  as  a  classical  oscil- 
lator. Returning  to  the  problem  in  193327,  Lenz  incorporated  line  shifts  and 
asymmetries  into  his  collisional  line  broadening  theory  by  using  finite  collision 
times.  Previously  all  radiation  interrupting  collisions  had  been  assumed  to 
occur  instantaneously. 

Michelson's  point  about  the  radiation  from  atoms  being  affected  by  atoms 
in  their  vicinity  was  taken  up  by  Holtsmark28,29  in  1919  and  developed  into 
the  theory  of  statistical  Stark  broadening.  Holtsmark  considered  the  effect 
of  a  particular  configuration  of  static  perturbers  surrounding  the  radiator  on 
the  shift  of  the  radiator's  characteristic  frequency.  He  then  averaged  over  the 
possible  configurations  to  produce  a  broadened  line  profile.  If  the  perturbers 
are  electric  monopoles  this  theory  is  referred  to  as  Stark  broadening  theory  after 
an  idea  originally  suggested  by  J.  Stark30  in  1914.  Holtsmark  also  developed 
similar  theories  for  perturbing  dipoles  and  quadrupoles. 

This  brings  us  to  a  point  in  the  historical  development  of  line  broadening 
theory  where  we  see  the  emergence  of  the  primary  theoretical  points  of  view 
as  to  the  physical  causes  of  line  broadening  in  plasmas.  These  are  Doppler 
broadening,  collisional  broadening  and  statistical  Stark  broadening.  Many  fur- 
ther developments  and  refinements17'18  have  occurred  since  these  "early  days" 
of  line  broadening  physics  but  we  will  not  follow  the  historical  development 
further.  We  will  jump  to  the  point  of  present  day  developments. 
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2.3  Theoretical  Formulation 

In  order  to  examine  the  line  radiation  from  a  plasma,  we  must  look  at  the 
power  emitted  by  dipole  radiation  from  an  emitter  immersed  in  the  plasma. 
At  the  outset  we  assume  that  the  plasma  consists  of  ionic  emitters  surrounded 
by  a  plasma  consisting  of  perturbing  ions  and  electrons.  This  will  allow  us 
to  model  the  plasma  line  radiation  by  examining  the  radiation  from  one  of 
the  emitters  and  averaging  over  all  possible  configurations  of  the  surrounding 
plasma.  This  is  equivalent  to  looking  at  a  large  number  of  plasma  radiators  in 
differing  environmental  conditions. 

The  resulting  power  radiated  by  dipole  radiation  from  an  emitter  immersed 
in  a  plasma,  as  a  function  of  the  frequency,  is  given  by^'"^ 


P(w)  =  |^I(u>)  (2.23) 

where  I(w),  the  line  shape  function,  is  defined  by 


I(u)  =  Y/^-^ab)\(b\de-i^\a)\2pa   .  (2.24) 

a,b 

Here,  a  and  b  refer  to  the  initial  and  final  states  of  the  total  radiator-plasma 
system,  respectively.  The  energy  difference  between  these  initial  and  final 
states  can  be  written  as  u>aj  =  (Ea  —  E^)/h.  The  dipole  moment  of  the  radiator 
is  given  by  d,  and  the  population  of  the  initial  states  is  given  by  pa  which  is 
an  eigenvalue  of  the  density  matrix  p,  where  p  —  Yli  Pi  I*)  (*!•  The  radiator 
center-of-mass  coordinate  is  given  by  R  and  the  wave  vector  of  the  emitted 
radiation  is  given  by  k  where  k  =  u/c.  The  eigenvectors  \i)  are  eigenfunctions 
of  the  system  Hamiltonian  H.  Free-free  and  free-bound  transitions  will  not  be 
examined  in  this  dissertation  so  we  consider  only  dipole  radiation  from  atomic 
ions. 
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Since  a  similar  formalism  can  be  used  to  describe  absorption  line  profiles, 
we  will  concern  ourselves  only  with  the  line  shape  function  I(u>).  Our  system 
will  consist  of  plasma  electrons,  ions  and  the  radiator  which  will  always  be 
referred  to  in  this  section  by  the  subscripts  "e",  "i",  and  "r",  respectively.  For 
simplicity,  we  also  set  h  =  1  for  the  remainder  of  this  section. 

We  can  express  the  Dirac  delta  function  as  an  integral  representation  that 
will  allow  us  to  incorporate  the  time  development  of  the  system  into  the  dipole 
moment  operators  (i.e.  the  Heisenberg  representation).  If  we  work  with  this 
representation,  we  have  the  advantage  that  all  of  our  approximations  must 
be  made  explicitly  on  the  system  Hamiltonian  which  will  be  contained  in  the 
expression  for  /(w).  This  gives 


1      f°° 

a,b 

=  -Re  / 
*       JO 


<fte'(w-(w« -<"»))< 


ik-R 


(b\de-tlcK\a) 


Pa 


(2.25) 


dtelwt4>{t) 


where 


4(t)  =  J2  H  eilfid\b)  ■  (b\  e^de-^^e-^pa  \a) 
a,b 

=  Y,(*\dk(0)\b)-(b\d_k(t)p\a)   . 


(2.26) 


0.6 


Here  we  have  used  the  fact  that  p  \a)  =  pa  \a)  and  e    lHi  \a)  =  e    Wl'0<  \a)  and  the 
definition  of  the  time  development  of  the  dipole  and  radiator  position  operator, 


d_k(t)  =  etHid(0)e 


(2.27) 


It  follows  that  df.  =  delkR.   We  will  consider  our  sum  over  states  a  and  b  to 
now  be  a  sum  over  a  complete  set  of  states  and  replace  the  sum  with  a  trace. 
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We  are  free  to  do  this  at  the  outset  and  only  look  at  the  elements  of  the  sum 
that  correspond  to  the  spectral  transitions  we  are  interested  in.  Now  we  have 


=  (4(0)  •  d_k(t)) 


(2.28) 


and  thus 


1  f°° 

I(u)  =  -Re  /      dtelU)t  (dk(0)  ■  d_k(t))  .  (2.29) 

*       JO 

To  address  the  problem  of  radiator  motion,  we  will  assume  that  the  radiator 
velocity  is  statistically  independent  of  perturbing  ion  and  electron  effects  as 
well  as  of  the  internal  state  of  the  radiator.  Consequently,  we  can  factor  the 
radiator  translational  motion  from  the  density  matrix.  This  gives 


P  =  PtrPr,i,e  ■  (2.30) 

Here,  ptr  is  the  density  matrix  for  the  radiator  translational  motion.  Assuming 
ideal  gas  behavior  for  the  radiators,  we  have 

/   mr   \  3/2       mrv2 

"Host)    e~w>  <2-31) 

where  mr  is  the  radiator  mass  and  v  is  its  velocity.  We  may  now  separate  out 
the  radiator's  translational  motion  to  give 

/(w)  =  ^Re  f°°dteiut   fdvptre^^-^id-  d(t))  .  (2.32) 

If  we  assume  that  the  radiator's  motion  is  constant  during  the  time  of  interest, 
then  we  can  write 

R(t )  =  R  +  vt  .  (2.33) 
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This  gives 

1  too  r 

I(lo)  =  -Re  /      dteiut   /  dvptre-ik-vt(d-  d(t))  .  (2.34) 

The  integral  over  v  may  be  performed  to  give 

I  TOO  /'in  \22tT.>2 

I(u>)  =  -Re  /      dteiujte~y-c)  ^  (d-  d(t))  .  (2.35) 

It  can  be  shown  with  the  aid  of  the  Fourier  convolution  theorem3"*  that  this 
expression  is  equivalent  to 

/oo 
dw'lD(u)-u>')Is(u')  ,  (2.36) 

-oo 

where  i£)(u;  —  u/  )  is  the  Doppler  distribution  given  by 


'rf-^-^-'PTSrfy  (2'37) 


where3 


*»  =  4  (^)   •  (2.38) 

Here,  mr  is  the  mass  of  the  radiator  and  c  is  the  speed  of  light.  The  line  shape 
function  for  the  static  radiator  is  given  by 


1  f°° 

Is{u)  =  -Re         dtelut(d-d(t))  .  (2.39) 

7T  JO 

This  result  is  the  usual  method  for  treating  the  Doppler  effect  but  we  see 
that  it  follows,  with  suitable  approximations,  directly  from  the  expression  for 
the  power  emitted  by  a  radiator.  The  interaction  of  radiator  motion  with 
plasma  perturbations  can  be  treated32'35  as  an  extension  of  the  above  devel- 
opment  by  not  making  the  approximation:    R(t)  =  R  +  vt.    We  will  employ 
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the  Doppler  convolution  given  by  Eq.    (2.36)  to  account  for  radiator  motion 
and  continue  with  the  development  of  /s(u>).  We  will  drop  the  subscript  from 
Isfo)  for  convenience  in  what  follows.  This  gives 


I(u)  =  -Re  [ 


00 


jut 


dt  elut  (d(0)  ■  d(t))  . 


(2.40) 


We  see  that  the  line  shape  problem  is  reduced  to  the  calculation  of  the 
system's  dipole-dipole  autocorrelation  function.  This  is  the  starting  point  for 
some  numerical  simulation  approaches.36  We,  however,  will  continue  with  the 
formal  development  of  the  theory  along  the  general  lines  of  the  work  of  Smith 
and  Hooper. 

We  will  see  below  that  it  is  possible  to  approximately  factor  the  density 
matrix  into  product  form  so  that 


P  =  PiPepr 


(2.41) 


Here  the  individual  density  matrices  are  for  the  three  subsystems:  the  plasma 
ions,  the  plasma  electrons  and  the  radiator.    We  will  examine  the  LTE  case 
where  the  density  matrix  becomes  the  simple  Boltzmann  factor  p  =  e~@H  /Tr  e~@H 
Consider  our  system  Hamiltonian 

H  =  H?+  tfr°  +  ffj  +  Viyr  +  Vij  +  Vfy  •  (2.42) 

Here  the  iP's  refer  to  the  kinetic  and  potential  energies  of  the  respective 
subsystems  so  that  Hf  =  T{  +  V{  and  H%  =  Te  +  Ve.  The  V^'s  refer  to  the 
potential  energy  of  interaction  between  the  subsystems  j  and  k.  We  make  the 
quasistatic  approximation  for  the  ions  so  we  assume  that  they  do  not  move 
appreciatively  during  the  atomic  radiation  time.    The  use  of  the  quasistatic 
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approximation  allows  the  interaction  of  perturbing  ions  with  the  radiator  to 
be  statistically  screened  by  the  rapidly  moving  plasma  electrons.  Consequently, 
we  will  replace  Vj  e  +  V{  with  a  screened  ionic  potential  and  as  an  approximation 
obtain 

V(  w  Vi  +  Vitt  .  (2.43) 

This  gives  us  the  effective  ionic  Hamiltonian  Hf'  =  T,-  +  V-.  What  precise 
form  this  screened  potential  should  take  has  been  the  subject  of  considerable 
work.  8™|4°  Here,  we  will  use  the  lowest  order  screening  approximation  for 
Vj:  the  Debye-Hiickel  potential.20  In  order  to  make  the  density  matrix  factor- 
ization, we  will  neglect  Vrie,  but  only  in  the  density  matrix.  This  amounts  to 
neglecting  radiator-perturbing  electron  correlations  at  the  initial  time.  This  is 
a  reasonable  approximation  as  long  as  kT  >  Vr<e.  Now  we  have 


rO/ 


HKHr  +  H»  +  B»  +  Vijr. 


(2.44) 


The  ion-radiator  interaction,  Vj  r,  will  be  expressed  as  a  multipole  expansion 
which  we  will  truncate.  For  the  density  matrix  we  will  keep  only  terms  through 
the  dipole.  The  monopole  term  has  only  ion  coordinates  and  is  incorporated 
into  Hj  to  give  H'^.  The  dipole  term  contains  both  ion  and  radiator  coor- 
dinates but  we  will  use  a  mathematical  device  to  eliminate  those  of  the  ions 
and  combine  it  with  H®  to  give  Hr(e).  Here  e*is  the  plasma  microfield  to  be 
discussed  below.  If  we  assume  that  the  remaining  effective  Hamiltonians  for 
the  three  subsystems  commute  amongst  themselves  the  density  matrix  is  now 
factorable  with  the  Hamiltonian  H  w  H[  +  Hr(T)  +  H®.  We  will  see  what  these 
terms  are  in  detail  shortly. 
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Turning  from  the  approximate  form  for  the  Hamiltonian  in  the  density 
matrix,  we  consider  the  terms  in  the  Hamiltonian  of  the  dipole  moment  time 
development  operator.  We  make  a  multipole  expansion41  of  the  radiator-ion 
interaction 


Vr^xm-d-m-l^Qij^1- 


•j 


dxj 


(2.45) 


The  electric  fields  and  potentials  refer  to  those  produced  by  the  plasma  ions 
at  the  radiator  which  we  take  to  be  at  the  origin.  The  effective  charge  of  the 
radiator  is  given  by  \  =  Z  —  a  where  Z  is  the  radiator  charge  and  a  is  the 
number  of  bound  radiator  electrons.  The  radiator  dipole  moment  operator 
is  given  by  d  and  the  components  of  the  radiator  quadrupole  moment  tensor 
are  given  by  Qij.  One  effect  of  the  quadrupole  term  is  the  production  of  an 
asymmetry  in  the  line  by  intensification  of  the  blue  portion  of  the  line  relative 
to  the  red.12  The  octupole  term  is  fourth  order  in  the  expansion  parameter  6 
compared  to  the  third  order  of  the  quadrupole  term.  For  the  special  case  of 
hydrogen-like  ions,  Joyce42  has  shown  that  the  octupole  term,  to  first  order  in 
perturbation  theory,  produces  symmetric  shifts  of  the  hydrogenic  energy  levels 
and  thus  will  not  contribute  significantly  to  the  asymmetry  of  the  spectral 
line.  Since  we  are  using  a  basis  set  for  this  calculation  that  is  made  up  of 
eigenfunctions  of  the  radiator  Hamiltonian  that  includes  the  field  dependent 
dipole  term,  we  take  the  quadrupole  term  as  the  first  order  correction  to  this 
Hamiltonian,  and  neglect  all  higher  order  terms. 

We  will  symbolically  represent  all  field  gradient  terms  dE{(0)/dxj  by  the 
generic  term  E^v  =  d^Ev.  Whenever  we  use  this  term  we  will  mean  all  rele- 
vant partial  derivatives.  Our  device  for  removing  the  ion  coordinates  from  the 
interaction  potential  VJ  r  involves  the  introduction  of  a  delta  function  in  the 
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manner  shown  below.  This  will  alow  us  to  separate  the  ion  coordinates  from 
the  radiator-perturbing  electron  subsystem  and  proceed  with  the  development 
as  shown  below.  Introducing  the  delta  functions  into  Eq.  (2.28)  gives 


cj>{t)  =  j  dej 


de^v  Trre{ 


PiS(e-  E)S(eflv  -  E^)d(0)  •  d(t)pre     .      (2.46) 


Here,  the  integration  variable  e^u  is  used  in  the  same  way  as  Euv  above. 

The  time  dependence  of  the  radiator  dipole  operator  is  given  in  terms  of  the 
total  system  Hamiltonian  H.  The  ion  coordinates  in  H  that  are  not  removed 
by  the  delta  function  will  commute  with  d(0)  in  the  static  ion  approximation 
and  cancel.43,44  This  leaves  the  Hamiltonian  JTe,r  that  has  only  radiator  and 
perturbing  electron  coordinates.  This  is  illustrated  by 


iH„tvr\\-iHert  _  JLtji 


d{t)  =  el"^d(0)e 


=  ethtd(0)  , 


where 


Her  =  H°  +  Ve,r  +  Hrtf-±Y,Qi,J^-r 


(2.47) 


(2.48) 


•J 


and  L  is  the  Liouville  superoperator  expressed  in  terms  of  a  commutator  acting 
on  an  arbitrary  operator  O  as 


LO  =  [Her,0)  . 


(2.49) 


We  have  now  successfully  separated  the  system  ion  coordinates  from  the 
radiator  and  perturbing  electron  subsystems.  This  is  expressed  by 


I(u>)  =   /  de      dtpVW(t,epV)J(uXtpV)  , 
where  W(?t€/tv)  is  the  ion  probability  function,  and 


(2.50) 
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J(a>,e,e^) 


7T     yo 


-JW< 


^eJWI(40)-J(0)e,r 


d<  eiw<Tr, 


«T-  e-*^ 4]  (2-51) 


1 


ImTr 


IT 


e,r 


d  ■  Re,r(u)(pe,rd) 


We  have  defined  the  resolvent  operator  as 


R, 


Jroo 
f      dteiute-iLt  =  [u-L}-1  . 
0 


(2.52) 


Re  and  Jm  refer  to  the  real  and  imaginary  parts  respectively.  The  ion  proba- 
bility function  W(e,  e^)  is  given  by 


(2.53) 


Wfcefu,)  =  Tn  [piS(e-  M)S(efU/  -  E^) 

=  (S(e-  E)6(e^  -  E^))i  . 
The  field  E  and  field  gradient  E^u  are  sums  of  contributions  from  all  of  the 
individual  system  ions.  We  will  return  to  the  evaluation  of  W(e,  t^v)  in  chapter 
3.   Recalling  that  pe<r  =  pepr,  we  can  perform  the  trace  over  plasma  electron 
coordinates  to  give 


J(u,etefw) 


1 


ImTrr 


TT 


d-  (Rer(u))e(prd) 


(2.54) 


We  see  that  the  effects  due  to  the  perturbing  electrons  have  been  isolated  in 
Rr(u)  =  (Rer(u))e.  Note  that  L  =  L{e,e^). 

Next,  we  proceed  with  the  evaluation  of  the  operator  Rr{w).    We  now 
introduce  a  coupling  constant  A  as  an  expansion  parameter.  Write 


L  =  Lq  +  XL j  , 


(2.55) 


where  L  corresponds  to  the  Hamiltonian: 


2S 


He>r  =  H0  +  XVj  =  Hi  +  A  (Ve,r  -  1  £  a^^.) 


(2.56) 


«j 


The  i^r  corresponds  to  Lq  and  Ve>r  +  Vfr  corresponds  to  Lj.  Here,  Qij  is 
a  function  of  the  radiator's  electron  coordinates  (for  single  or  multielectron 
radiators),  Fe,r  is  the  radiator  and  plasma-electron  interaction,  Vfi  is  the  ion- 
quadrupole  interaction,  and  H®r  =  Hr(e)  +  H®  is  the  combined  non-interacting 
radiator  plasma-electron  Hamiltonian  containing  the  ion-dipole  term.  We  will 
now  implicitly  define  an  electron  broadening  operator  M(u>),  by  setting 


Rr(u)  =  [w  -  Lr(€*,  e^i/)  -  M(w)] 


-1 


(2.57) 


where  Lr(e,  e^v)  is  the  Liouville  operator  corresponds  to  the  Hamiltonian: 


Hr(e,  e^v)  =  Hr(e)  -  -  J]  Qijfoj 


(2.58) 


«J 


Here,  .ffr(e)  is  the  Hamiltonian  of  the  field-free  radiator  with  the  ion-dipole 
term,  and  we  have  used  the  quadrupole  approximation  for  the  remainder  of 
the  interaction  of  the  plasma  ions  with  the  radiator.  Specifically, 


Hr(e) 


Hr  —  d  •  e 


(2.59) 


Note  that  Hr(e)  =  Hr(e,0). 

So,  using  Eqs.  (2. 52), (2.55)  and  (2.57),  we  can  write, 


[u  -  Lr(e,  e^)  -  M{uj)}    l  =  ([u  -  L0  -  AL/]"1^   .  (2.60) 

Next,  we  expand  M(u>)  in  powers  of  the  coupling  constant,  expand  both  sides 
of  Eq.  (2.60),  equate  coefficients  of  equal  powers  of  A,  and  thus  identify  the 
terms  of  the  M{u)  expansion. 
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The  right-hand  side  of  Eq.  (2.60)  can  be  expanded  in  a  Lippmann-Schwinger 
expansion  and  the  left-hand  side  can  be  expanded  in  a  Taylor  series  in  A.  First, 
we  will  expand  the  right-hand  side.  For  simplicity  in  notation  we  will  denote 
the  full  resolvent  operator  by  Re^r(uj)  =  R,  and  define  a  zero  order  resolvent 
operator  as 


Rq  =  [u  -  Lq] 


-l 


We  can  construct  an  iterative  equation  for  R  by  noting  that 


R  —  Rq  =  Rq 


Rq        —  R 


R 


or  explicitly, 


(2.61) 


(2.62) 


[u-Lq-(u-  L)) 


u>  —  L      u  —  Lq       u>  —  Lq 
This  can  be  written  equivalently  as, 


u 


-L 


(2.63) 


1  1  1       .r        1 

+  -  — —  XLr 


u  —  L        u  —  Lq       oj  —  Lq  uj  —  L 


or  simply  as  just 


R  =  Rq  +  XRqLjR 


(2.64) 


(2.65) 


This  is  the  Dyson  equation  for  this  system  in  Liouville  space.    This  implicit 
equation  for  R  can  be  iterated  to  give  a  series  expansion  for  R  in  powers  of  A: 


R  =  £  A^o^/Ao)"  • 
To  second  order  in  A  we  have 


(2.66) 
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R  =  R0  +  XRqLjRq  +  X2R0LiR0LiRq 


In  order  to  examine  what  (•  •  -)e  does  to  this,  look  at 


(2.67) 


Rr(u)  =  (R)e  ■ 
To  second  order  in  the  interaction  Liouville  operator  this  is 


(2.68) 


Rr(e>)  =  (Ro)e  +  \(RoLIR0)e  +  X2(R0LIRoLIRQ)e 


(2.69) 


Here  Rq  is  a  function  of  Lr  and  Le,  the  zero-order- radiator  and  perturbing- 
electron  Liouville  operators  with  Lr  corresponding  to  Hr(e)  and  Le  to  H® . 
Tighe44  has  shown  that  when  Lq  operates  on  objects  of  the  form  B  =  peb, 
where  b  operates  on  radiator  coordinates  only,  then 


(.Lq  )e  —  £" 


(2.70) 


On  expanding  Rq  this  leads  to 


-1  -  d0/ 


{R0)e  =  [u-Lr}-L  =  RJr(u;). 


(2.71) 


For  combinations  of  Rq  with  a  general  operator  C  in  Liouville  space  Tighe44 
has  shown  that 


(fl0£)e  =  R°r(C)e, 


(2.72) 


and 


(CRq)€  =  (C)R°r 


(2.73) 


31 

These  identities  allow  us  to  simplify  the  expression  for  the  resolvent.  We  obtain 

(R)e  =R°r  +  Ai?J(X7)ei??  +  X2R0r(LIR0LI)eR0r  .  (2.74) 

For  the  left-hand  side  of  Eq.  (2.60)  we  postulated  an  equivalent  form 
for  the  resolvent  that  contains  all  of  the  perturbing  electron  statistics  in  an 
implicitly  defined  operator  M(\,u>)  which  is  analytic  in  the  coupling  constant 
A.  Restating  this  we  have 


(R)e=  [u  -  Lr{e^u)  -  M{\,u)] 


-1 


(2.75) 


where  Lr  is  defined  above.  Expanding  M(A,u>)  in  powers  of  A  we  have 


M(A,w)  =  M(°)(w)  +  AMM(w)  +  A2M(2)(w)  + 


(2.76) 


To  expand  the  right-hand  side  of  Eq.    (2.75)  about  A  =  0  we  will  need  the 
following  easily  proved  operator  identity  for  some  operator  A  in  Liouville  space: 


ox  ~  A    d\A 

Expanding  -4-1(A)  about  A  =  0  we  have 


(2.77) 


,&4^(0)   ,   lx2^-1(0)+o(A3) 


^-1(A)  =  A~1(0)  +  A 


d\ 


+  2* 


OX2 


(2.78) 


In  this  case  we  have 


A(X)  =  [u  -  Lr(?,0)  -  XLrfcp,)  -  M(A,w)]   ,  (2.79) 

with  Lrfam,)  corresponding  to  V^.    Using  the  identity  of  Eq.    (2.77),  along 


with 
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dA(0) 


=  Lr^v)  -  M^H  , 


(2.80) 


and 


^  =  -2^»H 


A-\0)  =L-  Lr(e,  0)  -  M(°)(u;) 


-1 


we  can  evaluate  the  expansion  of  A   1  about  A  =  0.  We  obtain 


(2.81) 


(2.82) 


(R)e  =  A'1  =  ^(O)  +  A^-!(0)  [Lrfav)  +  mW(w)]  ^(O) 

-  A2{ A-\0)  [Xr,,-(c^)  +  AfW(«)]  ^(O)  [lr,e(6^)  +  M^H]  A"l(0) 

-  A-^O)^2)^)^-1^)}  +  0(A3)  . 

(2.83) 

Comparison  of  the  two  expansions  for  (R)e,  term  by  term  in  powers  of  A  up 
to  second  order,  allows  us  to  make  the  following  identifications  for  the  terms 
in  the  expansion  of  M(\,u>): 

A(°)  term 


At1)  term 


A(2)  term 


#(«)  =  0  , 


mM(l>)  =  (Lr,e)e  , 


M(2)(W)  =  (LjRqLj)  -  (Lj)eijO(Ij)e 


(2.84) 


(2.85) 


(2.86) 
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This  last  expression  for  M^(u)  can  be  further  simplified  if  we  expand  Lj. 
Recall  that  Lj  =  Xe,r  +  £j>  with  Le^  corresponding  to  Ve,r  and  £,-  r  to  Vf* . 
This  gives 


M^{u)  =  (LrteR0Lr,e)e  -  (Lr,e)eR0r(Lr,e)e 

+  Lr!i((R0  -  R°r)Lr,e)e  +  (Lr<e(R0  -  R°r))eLrji  . 
With  the  use  of  Eqs.  (2.72)  and  (2.73)  the  last  two  terms  cancel  to  give 


(2.87) 


MV\u>)  =  (Lr,eRoLr,e)e  ~  (Lrte)eR0r(Lr,e)< 


(2.88) 


If  the  dipole  approximation  is  made  for  Ve,r,  we  will  find  for  the  corre- 
sponding Liouville  operator,  (£e,r)e  =  0.  This  is  easy  to  see  by  noting  that 
in  order  to  evaluate  this  expression  we  need  to  calculate  (Ee)e  which  must 
average  to  zero  since  there  is  no  preferred  direction  in  the  plasma.  This  holds 
true  for  the  case  of  hydrogenic  or  multielectron  radiators.  The  first  surviving 
term  in  the  M(u>)  expansion  will  then  be  to  second  order.  So  to  lowest  order 
in  A, 


M(a;)  =  M(2)(u;)  =  (i:e,rJR0(a;)Lejr) 
where  Rq((jj)  =  [u  -  Lq]~1.  This  gives  us 


e  i 


(2.89) 


/•OO 

M(u)  =  -i  /      dte^iLese-^LesJe  ■  (2.90) 

JO 

This  operator  is  tetradic37  because  of  the  use  of  the  Liouville  operators  and, 
when  evaluated  as  a  matrix  element,  will  give45 
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M(u)ab,a>b'  =  -o 


hb1  2_j  daa"  '  da"a'G(&'ua"b) 


+  <W  \^  db'b"  '  dbnbG(-Auab„) 

b" 


(2.91) 


-  daa,  ■  db,b{G(  Au>ab,)  +  G(-Aua,b)} 
where  the  first  term  represents  electron  broadening  of  the  upper  state,  the 
second  term  represents  electron  broadening  of  the  lower  state  and  the  third 
term  represents  interference  effects. 

The  many-body  electron  broadening  function  G(Au>ab)  is  given  by 


fOO 

G(Acoab)  =   /      dteiAu^(Ee-Ee(t))ei 
JQ 


(2.92) 


iH?t  j?  (c\\o-iH0A 


with  Et(t)  =  elH^Ee(0)e-lh^  and  Auab  =  u  -  (ua  -  w6).  Thus  the  elec- 
tron perturber  effects  are  given  in  terms  of  an  autocorrelation  function  for  the 
plasma  electric  field  at  the  radiator  due  to  the  plasma  electrons.  We  have 
calculated  only  the  imaginary  part  of  M(u)  which  corresponds  to  line  width 
operator.  The  real  part  corresponds  to  the  line  shift  operator  and  yields  only 
a  small  shift44  which  is  ignored  here. 

In  our  approach  to  this  calculation  we  approximate  the  effects  of  corre- 
lations among  the  perturbing  electrons  by  the  use  of  G(Au>)  =  G(up)  for 
| Aw  |  <  u>p,  where  up  is  the  electron  plasma  frequency.  This  procedure  has 
been  shown  to  be  a  good  approximation  by  a  detailed  treatment  of  dynamic 
electron  correlations.  This  handling  of  electron  correlations  allows  us  to  use 
single  particle  Coulomb  waves  in  taking  the  trace  over  perturbing  electron  co- 
ordinates. The  Coulomb  waves  are  single  particle  solutions  to  Schrodinger's 
equation  with  a  central  1/r  potential  with  effective  radiator  charge  Xi  where 
X  =  Z  —  apt;  aft  being  the  number  of  bound  radiator  electrons.    G(Au)  is 
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then  evaluated  using  the  method  of  Tighe  and  Hooper.46  Other  theories47  also 
exist  for  the  evaluation  of  G(Au>).  It  is  instructive  to  point  out,  however,  that 
we  can  relate  this  electric  field  autocorrelation  function  to  the  charge  density 
autocorrelation  function48  which  gives  us, 

(Ee  ■  Ee(t))e  =  j  <Px  J  d*x' '^ (p(x) ptf ,  t)) e  .  (2.93) 

This  illustrates  the  relationship  of  the  plasma  charge-density  fluctuations  to  the 
spectral  line- broadening  problem.  Note  that  this  charge  density  correlation 
function  is  not  translationally  invariant  for  ionic  radiators  as  is  sometimes 
assumed.49  Translational  invariance  would  imply  a  uniform  perturbing  electron 
gas;  for  highly  charged  radiators,  this  assumption  would  ignore  aspects  of  the 
radiator-perturbing  electron  correlations. 

Altogether,  our  approach  leads  to  a  line  shape  function  of  the  following 
form: 

J(u>)  =       de      de^W^eft^J^u^e^nv)   , 


where 


J(u>,  e,  tpy)  =  --lmTrr  \d-  [w  -  Lr(e,  epV)  -  M(u)]    l  prd\  (2.94) 


or  in  matrix  element  form 


J(w,  ?,  tfu,)  =  --I™  J2  PadVa,  ■  dab  [u  -  Lr(e,  e^)  -  M{uj\  a}h,  ab  ,   (2.95) 


aba'b' 


where  a  and  b  refer  to  initial  and  final  radiator  states  respectively. 
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The  above  equations  for  the  line  intensity  l(u)  are  valid  for  the  general 
case  of  multielectron  radiators  with  combined  upper  and  lower  state  electron 
broadening  effects  and  Stark  shifts.  In  this  work  we  will  be  considering  hy- 
drogenic  ion  resonance  lines  which  have  the  ground  state  as  the  lower  state. 
For  hydrogenic  ions  the  ground  state  is  not  subject  to  electron  broadening  and 
only  very  weakly  affected  by  the  plasma  electric  microfield.  Therefore  we  will 
ignore  lower  state  broadening  and  interference  effects.  This  will  leave  only  the 
upper  state  term  of  the  three  terms  in  the  matrix  element  expression  for  the 
electron  broadening  operator  M(u)  given  by  Eq.  (2.91).  This  gives 


M{u) 


i 

3 


£  daa"  ■  da„a,G(Au>a,,f)  ,  (2.96) 

a" 
where  subscript  /  refers  to  the  single  lower  radiator  state. 

Since  lower  state  ion  produced  Stark  shifts  are  unimportant  here,  the  ra- 
diator Liouville  operator  L^e^)  in  Eq.  (2.95)  can  be  evaluated  to  give 

J(w,  e,  e/u,)  =  --Im  J2  Padaf  ■  dfa,  [u  -  (Hr(e,  e^)  -  uf)  -  M(u)]  "J, 

aa' 

(2.97) 

where  Uf  is  the  lower  state  energy  of  the  radiator.  We  note  that  since  there  is 

only  one  lower  state  the  sum  over  b  is  absent. 

For  our  case,  the  dipole  and  quadrupole  operators  become  single  electron 

quantities  with 


d  —  —ef 


Qi,j  =  -e(3xiXj  -  r28ij) 


(2.98) 
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Here  r  and  a;,-  refer  to  radiator  electron  coordinates.  The  corresponding  mul- 
tielectron  quantities  are  just  sums  of  terms  similar  to  these  for  each  radiator 
electron. 

Any  intrinsic  line  width,  such  as  the  resonance  width  induced  by  the  pres- 
ence of  an  electric  field  (see  Chapter  3),  can  be  approximately  included  in  the 
line  profile  formalism  by  replacing  the  affected  level  energy  £  by  a  complex 
resonance  energy50  E.  For  level  j  this  gives 

Ej  w  Ej  -  i-1  ,  (2.99) 

where  Ej  is  the  energy  at  the  center  of  the  resonance  and  the  presence  of 
iTj/2  assumes  a  Lorentzian  shape  for  the  resonance.  Tj  is  the  full  width  at 
half  maximum  of  resonance  j.  We  have  chosen  the  minus  sign  in  Eq.  (2.99) 
to  be  consistent  with  the  fact  that  the  electron  broadening  operator  M(u>)  is 
proportional  to  —  i.  A  choice  of  Ej  —  i-£  would  correspond  to  a  proportionality 
factor  of  +i.  The  use  of  this  complex  energy  is  only  an  approximation  since  we 
are  attempting  to  describe  the  time  dependent  behavior  of  a  resonant  state  in 
a  stationary  state  formalism.  We  use  this  procedure  to  approximately  include 
the  field  induced  resonance  width  for  the  upper  states  of  the  Lyman  series  in 
hydrogenic  argon. 

When  considering  the  comparison  of  theoretical  line  profiles  to  experimen- 
tal results,  it  is  important  to  remember  the  factor  u4  in  the  relation  for  the 
total  power,  P(u)  =  Const.  u4I{u).  If  the  radiation  energy  is  high  enough  and 
the  line  shape  is  broad  enough,  this  factor  will  not  be  approximately  constant 
across  the  line  shape  and  I(u)  must  be  redefined  to  include  it.   Otherwise  an 
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additional  source  of  line  asymmetry  will  be  neglected.  For  experimental  com- 
parison, it  is  also  important  to  include  instrumental  broadening5  as  well  as 
opacity  effects51  if  significant. 


CHAPTER  III 
HIGHER  ORDER  FIELD  EFFECTS 

In  the  preceeding  chapter  we  have  reviewed  some  fundamentals  of  dense 

plasmas  and  formulated  a  theory  of  spectral  line  shapes.  We  will  now  examine 

some  of  the  processes  and  effects  that  can  make  important  contributions  to 

these  line  shapes  in  the  high  density  regime. 

3.1  The  Ion-Quadrupole  Effect 
Previously  we  introduced  the  ion  coordinate  joint  probability  function 
W(e,  e^j/)  to  describe  the  joint  probability  distribution  for  the  microfield  e  to- 
gether with  the  field  gradients  e^.  Recall  that  for  notational  convenience,  e^, 
represents  all  of  the  field  gradient  terms.  Writing  out  this  probability  function 
explicitly  we  have 

W(e,ep,)  =(8(e-  E)8(exx  -  Exx)8(eyy  -  Eyy)8(ezz  -  Ezz) 

x  8{exy  -  Exy)8(eyz  -  Eyz)8(ezx  -  Ezx) ),-  (3.1) 

=j-c  J  e-PV8(e-  E)8(exx  -  Exx)  ■  ■  ■  d™r  , 

where  <■••>,•  is  the  ensemble  average  over  ion  coordinates,  Zc  is  the  configu- 
rational  partition  function  for  the  ion  subsystem,  N  is  the  number  of  ions,  and 
V  is  the  potential  energy  of  the  ion  subsystem  in  the  presence  of  the  charged 
radiator.  The  quantities  E  and  E^  are  the  many  particle  electric  field  and  field 
gradients  at  the  radiator  and  are  functions  of  the  perturbing  ion  coordinates. 
Using  an  integral  representation  of  the  field  gradient  delta  function, 
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S(exx  -  Exx)  =  i-  /  daxx  e-»^(^-^)  ? 


(3.2) 


gives 


W&  e^)  =  TtAktF  I  d6°  C-^»e"+-)  /  d*Nr  ePV6(e-  M*~E'*+'") 

(27r)°Zc  J  J 

=  tAjQ(«)  f  dQae-i(<T"€^+cryy€yy+-\eiQ)e 
=  Q{e)P{e\e^)  , 


where 


(3.3) 


0    _   <JXXEXX   +  OyyEyy   +  <TZZEZZ   +  GXyEXy  +  OyZEyZ   +  (TZXEZX    .  (3.4) 


We  have  defined  the  conditional  averaged  quantity 


■0,   _/^e-/^-£ye 


(3.5) 


This  is  the  ensemble  average  of  e'0  where  all  included  ion  configurations  are 
constrained  to  have  the  microfield  value  e.  It  is  normalized  to  give  one  if  0  =  0. 
By  defining  the  average  in  this  way  we  have  retained  the  benefits  of  using  the 
ordinary  plasma  microfield  function  Q(e)  to  describe  the  uniform  field  at  the 
radiator  and  we  are  now  free  to  approximate  the  higher  order  field  gradients  if 
we  choose.  The  function  P(e*|e^j/)  is  the  conditional  probability  for  €pV  given 
e.  We  take  the  field  to  be  in  the  z  direction  for  simplicity  and  without  loss  of 
generality. 
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Direct  evaluation  of  this  constrained  average  is  quite  difficult  if  we  wish  to 
retain  all  of  the  ion  interactions.  Consequently,  we  will  make  a  simple  approx- 
imation to  this  function  in  order  to  proceed.  We  will  expand  the  exponential 
in  a  cumulant  expansion^2  by  taking 


•DO 


(eiQ)f  =  exp     £^<?J,  (3.6) 

\n=l  / 

where 

d  =  (0>f 

C2  =  ((0  -  (6), )2)e  (3.7) 

C3  =  ((6  -  (6)e )3)f  . 
Higher  order  terms  are  more  complicated.  In  terms  of  the  expansion  parameter 

8,  which  we  introduced  in  the  multipole  expansion  (see  Eq.    (2.22)),  we  have 

C\  ~  S6  and  C2  ~  8°.  We  see  this  by  noting  that  0  goes  as  dEi/dxj  which  goes 

as  1/r3  and  1/r  is  of  order  8.  For  our  highest  density  of  ne  =  1  x  1025cm-3, 

this  gives  C\  ~  0.043  and  C2  ~  0.0018  for  the  n  =  3  level.  The  C2  term,  which 

is  actually  the  variance  of  the  quantity  Theta,  is  therefore  in  the  neighborhood 

of  4%  of  the  lower  order  term  C\ .  Here,  we  will  keep  only  the  C\  term  in  our 

approximation.  This  gives  us 

(eie}(  u  el'<e><  .  (3.8) 

Inserting  this  into  Eq.  (3.3)  gives 

=  Q(?)6(exx  -  (Exx)()8(eyy  -  (Eyy)€)  ■  •  •  (3-9) 
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This  gives  us  for  the  approximation  of  the  conditional  probability: 

P(?M  «  K^u  -  (E^)f)  .  (3.10) 

We  have  developed  a  systematic  expansion  for  the  conditional  probability  func- 
tion and  have  examined  the  magnitude  of  the  expansion  terms.  We  can  identify 
the  first  term  of  this  expansion  as  the  approximation  used  by  Joyce.12  Return- 
ing to  the  expression  for  the  spectral  line  intensity,  Eq.  (2.50)  now  gives 

/(«)  =  J dt  J dt^Q^P^t^JiuX^u)  •  (3.H) 

The  microfield  function  Q(e)  does  not  depend  on  the  direction  of  e  but 
only  on  its  magnitude,  e.  We  can  therefore  preform  the  angular  integration 
and  define  a  new  microfield  distribution  function  as 

P(e)  =  4Tve2Q(e)  .  (3.12) 

Now,  the  intensity  becomes 

I(u)  =   fde  f  de^  P(e)P(e\eliV)J(u,  e,  e^).  (3.13) 

Using  our  approximation  for  P(e\e^)  we  can  carry  out  the  e^,  integrations  to 
obtain 

I(u)  =  / "ifcP(e)J(a»,«,(E^>c)  .  (3.14) 

The  field  gradient  terms  in  J(u>,e,  e^u)  have  now  been  replaced  with  their 
constrained  averages  that  depend  only  on  e  so  we  can  write  J(uf,e,  (EuV)€)  = 
J  (w,  e),  where  the  prime  indicates  the  replacement  of  the  field  gradient  by 
its  constrained  average.     For  the  remainder  of  this  work,  we  will  drop  the 
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prime  with  the  understanding  that  we  actually  mean  J'{u,  e).  Our  use  of  the 
approximation  for  P{e\e^v)  has  lead  to  a  significant  simplification  of  the  line 
shape  formula  that  still  retains  the  approximate  effect  of  the  field  gradients. 

In  order  to  proceed  we  will  need  to  evaluate  this  constrained  average  of 
the  field  gradient.  From  the  ensemble  of  all  possible  ion  configurations  in  the 
presence  of  the  radiator  we  choose  only  those  configurations  that  give  a  certain 
field  magnitude  e  at  the  origin.  If  we  average  a  quantity  F  over  this  subset 
of  configurations,  we  have  constructed  the  constrained  field  average  that  we 
introduced  above.  We  denote  this  average  by 


(F)e  = 


(3.15) 


where  Q(e)  is  the  plasma  microfield  function  and  (•  •  •)  is  the  complete  ion 
subsystem  ensemble  average.  This  average,  as  we  have  said,  is  such  that  for 
F  =  1,  (F}(  =  1.  We  will  show  that  it  is  possible  to  relate  this  constrained 
average  to  a  functional  derivative  of  the  plasma  microfield  function.  This  will 
allow  us  to  evaluate  the  constrained  average  in  terms  of  available  simplified 
models  for  the  microfield  function  Q(e). 

For  our  purposes  we  will  take  F  to  be  an  additive  function  of  the  ion 
coordinates. 


F  =  Y,f(ri)  =  Jdrf(r)n(r), 


(3.16) 


where 


n(r)  =  J2  6(f  ~  fi) 


(3.17) 
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Note  that  J drn(r)  =  N{,  the  total  number  of  ions  in  the  system.   So  we  can 
now  write 


=  j  drf(r)nig(r;  e)  , 


where  we  have  defined 


^on.*<w-fl 


WO 


On  expanding  the  delta  function  we  have 


(3.18) 


(3.19) 


"  i 


,— i\-€i 


J\E\ 


W<*  t)  =  ^  J  J^  «-""Wr>n  .  (3.20) 

We  would  like  to  relate  the  average  on  the  right-hand  side  of  Eq.  (3.20)  to  a 
functional  derivative  of  the  plasma  microfield  function.  To  do  this  we  introduce 
G(A),  the  generating  function  of  the  microfield  distribution,  where 


i\E\ 


G(\)  =  \nQ(\)  =  \n{etAI") 


If  we  define  a  function  (f>,  such  that 


(ft  =  iX  •  E(r)  , 


then  G(X)  becomes  a  functional  of  <f>,  or 


(3.21) 


(3.22) 


G(\)  =  G[<f>]  . 

— * 

Now  we  can  cast  G(\)  in  the  form  of  a  functional  derivative,  where 


(3.23) 
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1 

G(X)  =  / 
Jo 


dl 


dG{l) 
dl 


'1*1 


dr 


d<f>  8G[<j)) 
dl    Scj> 


dfil  ■  E(r)- 


8G[J>] 


S(il  ■  E(f)) 
where  I  =  IX  and  A  is  the  unit  vector  in  the  direction  of  A. 

Alternately,  from  the  definition  of  G(A)  in  Eq.  (3.21),  we  have 


(3.24) 


r 

G(\)  = 

Jo 


dl 


din  Q(f) 


dl 


A      (iX  •  gcff) 


(eUE) 


il-E\ 


dl 


iX  •  (ge^fj 


Note  that  we  can  write, 


JV 


a-E=    f  driX-E(f)J26(r-fi) 

=    I  driX  ■  E(r)n(r). 
J 

Using  Eqs.  (3.25)  with  (3.26)  gives 


(3.25) 


(3.26) 


lnQ(A)=  /   J/  fdriX  ■£(&&££- 

JO        J  teilE\ 


(ct/.£7) 


(3.27) 


Comparison  of  Eqs.  (3.24)  and  (3.27)  gives  us  the  connection  between  the 
constrained  average  considered  here,  (F)£ ,  and  a  functional  derivative  of  the 
microfield  generating  functional  G[<f>\. 


{n{r)e*-S)  SG[<f>] 


(e«-E)  S(il  ■  E(f)) 

which  on  substitution  into  Eq.  (3.20)  above,  gives 


(3.28) 


4G 


^--mf^-^wiky      (3'29) 

Substitution  into  Eq.  (3.18)  above  gives  us  the  relationship  we  want: 


^-A/'Ww^w^iSfe-     (3-30) 

This  explicitly  shows  that  we  can  relate  the  constrained  average  to  a  functional 
derivative  of  a  microfield  generating  function  G(X)  =  \nQ(\).  We  are  free  to 
use  whatever  functional  model  we  choose  for  Q(\). 

As  an  aid  in  evaluating  these  constrained  averages  we  define  the  transform 
of  g(r;  e)  by 


g(f;A)Sl      *3t]       ■  (3.31) 

From  Eqs.  (3.20)  and  (3.28)  we  see  that  the  transforms  are  related  by 


^  *)  =  ^  J ^  e~iXrQWm  A)-  (3.32) 

Thus  a  knowledge  of  g(r;  A)  obtained  from  Q(\)  will  lead  to  g(f;  e)  through 
the  transform  and  then  to  (F)e  through  Eq.  (3.18). 

To  evaluate  the  ion-quadrupole  term  recall  that  the  perturbed  radiator 
Hamiltonian  is 


Hr(e,  (£^)e)  =  Hr(e)  -  ±  £  Qij(Eij)e  «  (3-33) 

where  the  ion-quadrupole  term  is  just  Vrj.  To  evaluate  the  (E{j)e  terms  note 
that  when  the  angular  integrations  are  performed,  the  y>,  or  azimuthal,  integrals 
lead  explicitly  to42 
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{EXy)e  =  (Eyz)e  =  {Ezx)e  =  0 


(3.34) 


Together  with  the  property  djEj  =  d{Ej,  this  leads  to  the  simplification 


J2  Qij(Eij)e  =  Y,  Qii(Eii)c  ■  (3-35) 

With  the  field  in  the  z  direction  there  is  no  difference  between  the  x  and 
y  directions  so  (Exx)(  =  (Eyy)t.  Additionally,  with  the  tracelessness  of  the 
quadrupole  tensor4*  we  have 


/  t  Qi,j(Eij)e  =  Qzz  {{Ezz)i  -  (Exx)€) 


'J 


=  ~Qzz  (<£«>< -|<V.£)£ 


(3.36) 


If  the  ion  interaction  potentials  are  screened  we  have  V  •  E  ^  0  due  to  the 
continuous  charge  distribution.  However,  the  form  of  the  multipole  expansion 
we  have  employed  made  use  of  V  •  E  =  0  at  the  origin  which  is  not  strictly  true 
for  our  screened  ion  potentials.  We  reconcile  these  two  facts  by  noting  that  if 
we  restrict  the  screening  charge  density  of  the  perturbing  ions  to  lie  outside 
the  Bohr  orbit  of  the  radiator's  bound  electron,  we  will  have  V  •  E  —  0  at  the 
origin.  Tests  of  the  computation  of  the  field  gradient  term  show  that  its  value 
upon  taking  V  •  E  —  0  is  approximately  20%  less  than  that  obtained  on  taking 
V  ■  E  /  0  at  a  field  value  corresponding  to  the  maximum  of  the  microfield 
probability  function  at  ne  =  1  x  1025cm-3  and  kT  =  800  eV.  The  difference 
is  less  for  lower  densities.  In  light  of  the  uncertain  nature  of  the  actual  form 
of  the  ion  potential,  this  does  not  seem  unreasonable.   Consequently,  we  will 

— *  — # 

take  V  ■  E  =  0  in  the  evaluation  of  the  field  gradient  term  in  agreement  with 
the  form  of  our  multipole  expansion.  This  amounts  to  requiring  that  we  have 
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no  penetration  of  the  radiator  by  the  perturbing  ions  or  their  screening  charge 
density.  For  the  ion-quadrupole  term  this  gives  us 

Vr,i  =  --Qzz{Ezz)c,  (3.37) 

with  the  radiator  Hamiltonian  given  by  Hr(e)  +  Vr  ,-.  To  evaluate  Vr  j  we  must 
calculate  the  constrained  average  term  (Ezz)e,  and  to  do  this  we  must  pick  a 
model  for  the  microfield  function  Q(e). 

In  this  dissertation  we  choose  the  adjustable  parameter  exponential  approx- 
imation (APEX)53-54'55-56  for  the  microfield.  It  has  the  advantage  of  incorpo- 
rating correlations  between  the  perturbing  ions  while  retaining  the  functional 
simplicity  of  the  independent  particle  (IP)  model42  of  the  microfield.  Joyce42,12 
used  the  IP  model  to  evaluate  the  ion-quadrupole  effect  which  we  will  see  leads 
to  an  overestimate  of  its  magnitude  due  to  the  neglect  of  the  inter-ion  correla- 
tions. We  will  compare  these  two  models  along  with  a  nearest  neighbor  (NN) 
approximation  for  the  field  gradient. 

The  APEX  model  for  the  microfield  is  based  on  the  independent  particle 
model  but  with  an  effective  inverse  screening  length,  a.  The  field  is  effectively 
renormalized  in  this  way  by  requiring  that  two  conditions54  be  met:  the  so 
called  second  moment  condition  and  the  local  field  condition  (see  Appendix 
B).  The  transformed  conditional  distribution  function  g(r ;  A)  for  the  APEX 
microfield  distribution  is  given  by  (see  Appendix  B) 

g(f ;  A)  =  g(r)eiX-&&  ,  (3.38) 

where  g(r)  is  the  equilibrium  radial  distribution  function  and  E*(f)  is  the 
APEX  renormalized  field.  This  field  is  given  by 
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E*{f)  =  E(r)(l  +  ar)e-Qr,  (3.39) 

where  a  is  determined  by  the  second  moment  of  Q(e),  and  E(f)  is  the  screened 
field  produced  by  the  plasma  ions.  We  calculate  the  radial  distribution  function 
g(r)  by  integral  equation  methods  due  to  Rogers.39  Lado  and  Dufty57  have 
shown  that  this  method  of  calculating  g(f ;  e)  by  using  the  APEX  microfield 
underestimates  the  inter-ion  correlations  in  some  cases.  In  particular,  for  the 
case  of  a  hydrogen  plasma  with  I";;  =  10  at  a  field  value  of  e/e0  =  0.4  (e0  = 
e/ro,i)'  tne  APEX  calculation  underestimates  the  maximum  value  of  g(f ;  e)  by 
a  factor  of  two.  We  retain  the  APEX  model,  however,  because  of  its  simplicity 
and  significant  improvement  over  the  IP  model  of  Joyce12,42  which  ignores 
ion-ion  correlations  entirely. 

With  the  APEX  model  for  g(r ;  e)  obtained  from  Eqs.  (3.32)  and  (3.38)  we 
can  evaluate  the  constrained  average  by  using  Eq.  (3.18)  to  obtain 


(F)e  =  ^  J  drm9(r)Q(i-  E*(r)),  (3.40) 

where  it  is  understood  that  the  microfield  distributions,  Q(e),  are  in  the  APEX 
approximation. 

For  the  ion-quadrupole  effect,  the  constrained  average  we  need  to  evaluate 
is  (Ezz)£.  In  particular,  we  have 


i 

In  general,  the  one  perturber  derivative  term  is  given  by 


(3.41) 


^(r)  =  Ze^(l  +  l)(W-3^e) 


(3.42) 


50 


so  for  our  case  we  have 


,r/A 


e2Z(r)  =  Ze-^(l  +  ^)(l-3cos20)    . 
The  microfield  function  in  the  APEX  approximation  is54 


with 


(3.43) 


I  TOO 

QaPExW  =  ITT  /      dkksm(ke)en<h^  ,  (3.44) 


roo 


E(r) 


hi(k)  =  4tt  jf   drr'g(r)j^r  ( j0  (kE*(r))  -  1)  ,  (3.45) 

where  Jq(x)  =  sin(x)/x  is  the  zero  order  spherical  Bessel  function  of  the  first 
kind.  From  now  on  when  we  speak  of  Q(e)  we  will  implicitly  assume  the  APEX 
approximation  and  suppress  the  subscript.  Using  this  approximation  we  can 
evaluate  Eq.  (3.40)  for  the  case  of  F  —  Ezz  to  obtain 


,-r/A 


(Ezz)€  =  -^-  Hdkkh^W   rdrr2g(r)^-(l-r/X)IQ(r).  (3.46) 
2iriQ{t)  Jq  Jo  H 

The  angular  integral  Iq  has  been  evaluated  by  Joyce42  and  is  given  by 


IQ(r)  =   fdtyl  -  3cos20)jo  (*  |*-  $*{r)\) 
=  S*h(ke)j2(kE*(r))  , 


(3.47) 


with  dQ,  =  d(pd(cos  9). 

To  evaluate  (Ezz)€  numerically,  we  cast  its  variables  in  unitless  form  by 
taking 
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I  =  e/e0 

ft   — -  £q/C 

a  =  r0>i/X 

^  =  r/r0,i  . 


(3.48) 


where 


t0  =  4"  .  (3.49) 

0,1 


We  obtain 


(E2z)c  =  -«  I  4- 


/0°°  d~kPen^Ck)j2(~k~e)  /°°  dxg(x)*=f-(l  +  ax)j2  (kE*{x)) 


rZjJ  f°°  dk~kenih^k)  sm{k~e) 


(3.50) 


with 


*l(*)  =  -  f°  dxx2g{x)^\  (j0(kE*(x))  -  l)   .  (3.51) 

In  this  form,  the  equation  for  (Ezz)(  can  be  directly  evaluated  numerically. 
Following  Joyce  we  implement  this  function  in  the  line  shape  calculation  by 
the  use  of  a  Pade  approximate.  For  a  particular  case,  (EZz)e  is  calculated  for 
a  mesh  of  e  values.  The  Pade  function  is  then  fit  to  the  calculated  values  by 
minimizing  the  squared  difference  between  the  two  while  varying  the  coeffi- 
cients of  the  Pade  function.  The  resulting  approximation  then  yields  (Ezz)( 
for  any  value  of  the  scaled  field  up  to  e  =  20  with  an  accuracy  of  about  5% 
compared  with  the  values  calculated  from  Eq.  (3.50).  This  is  quite  sufficient 
for  calculating  most  line  shapes.  If  more  accurate  values  of  the  field  gradient 
are  required,  the  directly  calculated  values  can  be  used. 
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We  would  like  to  examine  the  behavior  of  the  field  gradient  term  (Ezz)€  in 
the  limit  of  large  and  small  values  of  e.  For  e  — »  0  it  is  easy  to  show  that  the 
field  dependence  follows  the  limit  of  J2(ke)  to  lowest  order  in  e.  This  will  give 


Um(Ezz)(  =  -C~e2, 

€—►0 


(3.52) 


where  C  is  a  constant  that  depends  on  ne,  T  and  Z,-.  Thus,  the  zero  field  limit  is 
zero.  This  is  easy  to  see  physically.  If  e  =  0  there  will  be  no  preferred  direction 
about  the  radiator  and  the  fluctuations  in  the  field  gradient  Ezz  will  sum  to 
zero  upon  averaging  over  the  ensemble  of  ion  configurations.  The  e  — >  oo  limit 
is  somewhat  more  complicated  so  we  will  examine  the  behavior  of  a  simple 
nearest  neighbor  model  for  (Ezz)€ .  The  large  field  behavior  of  the  microfield 
at  a  charged  radiator  is  governed  by  a  single  nearest  neighbor  ion.  °'59  For  the 
strongest  fields,  the  nearest  neighbor  will  be  significantly  inside  the  screening 
length  Xq  e  so  screening  will  be  negligible.  To  approximate  (Ezz)£  let  the  field 
gradient  at  the  radiator  be  that  produced  by  a  single  ion  located  at  the  distance 
from  the  radiator  sufficient  to  give  a  field  value  at  the  radiator  of  e.  This  gives 
us  the  nearest  neighbor  approximation  of 


<*«>e.™  =  "^  (£) 


J-      *3/2 


(3.53) 


for  large  e. 

In  Fig.  (3.1)  we  compare  three  different  models  for  (Ezz)(.  For  purposes  of 
comparison,  we  take  (e/rot-3)  =  1.  We  look  at  the  APEX  model  which  includes 
ion-ion  correlations,  the  independent  particle  model  of  Joyce  which  ignores 
ion-ion  correlations  and  the  nearest  neighbor  model  for  large  fields  as  discussed 
above.  We  see  that  the  APEX  model  gives  a  field  gradient  whose  magnitude  is 
consistently  smaller  than  that  of  the  other  two.  This  is  to  be  expected  since  the 
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ion-ion  correlations  tend  to  cause  the  plasma  ions  to  repel  each  other  and  thus 
to  be  less  likely  to  converge  on  the  radiator  and  create  a  large  gradient.  The 
ion-ion  correlations  will,  in  general,  lessen  plasma  ion  microfields  by  causing 
a  rough  ordering  of  plasma  ions.  The  nearest  neighbor  model  produces  the 
largest  gradient  due  to  its  assumption  of  all  charge  located  on  the  z  axis.  It 
also  includes  no  correlation  effects.  The  IP  model  of  Joyce  falls  in  between 
because  it  includes  radiator-ion  correlations  but  no  ion-ion  correlations. 

In  Fig.  (3.2)  we  examine  the  behavior  of  the  APEX  model  for  the  field  gra- 
dient at  several  values  of  the  electron  number  density.  For  ne  =  1  x  1024cm-3 
the  average  Ezz  gradient's  magnitude  is  less  than  that  of  the  nearest  neighbor 
model  as  was  pointed  out  above.    As  the  density  increases  the  APEX  field 
gradient  approaches  that  of  the  nearest  neighbor  model  for  all  relevant  val- 
ues of  the  field  (in  this  case  we  refer  to  a  range  of  field  values  of  e/e0  =  0  to 
5).  For  a  particular  density,  we  can  see  that  the  gradient  also  approaches  the 
nearest  neighbor  value  as  the  field  increases.  This  behavior  can  be  explained 
by  recalling  that  the  presence  of  large  fields  at  the  radiator  is  generally  due  to 
the  location  of  a  nearest  neighbor  perturbing  ion  close  to  the  radiator.    The 
large  field  behavior  at  constant  density  should  be  due  to  a  nearest  neighbor 
ion  and  ion-ion  correlations  play  little  role  in  the  interaction  of  two  nearest 
neighbors  at  these  densities.    The  probability  is  small  for  having  a  positively 
charged  perturbing  ion  close  to  the  positively  charged  radiator.  It  is  even  less 
probable  to  have  two  ions  close  to  a  positively  charged  radiator;  and  we  need 
two  for  ion-ion  correlations  to  be  important.  The  approach  of  (Ezz)e  toward 
the  nearest  neighbor  value  as  density  increases  should  be  due,  simply,  to  the 
ions  being  forced  closer  to  the  radiator.  With  decreased  interparticle  spacing 
the  field  gradient  naturally  increases. 
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3.2  Atomic  Data  by  Perturbation  Theory  Solution 
To  evaluate  the  radiator  subsystem  matrix  elements  necessary  for  the  cal- 
culation of  the  spectral  line  shape  /(a;),  we  are  free  to  use  any  complete  set  of 
basis  states  that  spans  the  Hilbert  space  of  the  system.  A  true  representation 
of  an  operator  in  this  basis  will  thus  consist  of  a  matrix  whose  dimension  is  the 
same  as  the  number  of  basis  vectors  necessary  to  span  the  particular  Hilbert 
space.  For  our  case  this  dimension  is  a  nondenumerable  infinity.  However,  a 
reasonable  representation  of  the  operator  can  be  obtained  by  using  a  truncated 
set  of  basis  vectors.  For  weak  fields,  a  reasonable  representation  of  the  Hamil- 
tonian  for  a  radiator  in  the  state  with  principal  quantum  number  n  can  be 
obtained  by  using  eigenvectors  of  the  field-free  radiator  Hamiltonian  for  states 
n  and  n  +  1  as  basis  vectors.  "  For  stronger  fields,  the  need  to  include  other 
states,  including  positive  energy  states,  will  be  necessary.  In  fact,  the  reso- 
nance nature  of  the  strong-field  states  can  not  be  described  without  including 
positive  energy  states  to  account  for  the  unbound  character  of  field  dependent 
solutions  that  allow  for  the  possibility  of  tunneling.  In  order  to  include  these 
effects  we  will  employ  an  alternative  method  using  a  basis  set  that  consists  of 
the  eigenvectors  of  a  radiator  Hamiltonian  that  includes  the  ion-dipole  term 
from  the  perturbing-ion  multipole  expansion.  This  Hamiltonian  is 

Hr(e)  =  H?-d-e.  (3.54) 

Hr(e)  is  diagonal  in  this  field  dependent  basis  set.  By  using  this  Hamiltonian, 
we  will  automatically  incorporate  the  characteristics  of  the  resonate  nature 
of  the  states  into  our  basis.  Specifically,  we  will  use  Hr(e)  as  our  zero  order 
radiator  Hamiltonian.  The  first  order  correction  to  this  will  be  the  next  term 
in  the  perturbing-ion  multipole  expansion,  the  ion-quadrupole  term  Vr,\ 


55 


To  calculate  our  field  dependent  basis  we  solve  Schrodinger's  equation  in 
the  position  representation  with  the  Hamiltonian  Hr(e).  The  equation  is 

Hr(e)il>{r)  =  Eif>(r)  .  (3.55) 

The  natural  coordinates  to  use  for  this  Hamiltonian  are  parabolic  coordinates.61,62 
The  wave  equation  will  then  be  separable.  Assuming  the  plasma  microfield  to 
be  in  the  z  direction,  we  write  Eq.    (3.55)  in  atomic  units  (see  Appendix  A) 
and  scaled  to  Z  =  1  as 


V2  +  -  +  2E  -  2Fz)  if>(r)  =  0  , 


(3.56) 

where  F  is  the  magnitude  of  the  plasma  ion  microfield.  In  parabolic  coordinates 
(see  Appendix  C),  Schrodinger's  equation  becomes 


d  f,df\      JLA  (  ?£\       l  ^ 


Z  +  ridZYdS)       Z  +  1  drj  \' dr)  J       fr  dip2 
r^  +  2E  -  F(i  -  ri)\  $  =  0. 


(3.57) 


(€  +  *) 

We  exploit  the  separability  of  the  equation  by  looking  for  a  solution  of  the 


form: 


0(r)  =  fU)g(v)- 


,±im<p 


• ITT 


(3.58) 


We  define  new  scaled  variables  and  transform  the  two  resulting  separated  equa- 
tions. We  take 


x  =  £/n 

V  =  v/n 

A  =  (l/4)n3F 


(3.59) 


and  obtain  the  system  of  equations 


and 


with 


f  d     d 

<  — x 

\dx    dx 


5G 


2  9 

47  +  y^- 


x2  +  a\ 


/(*)  =  0  , 


d     d       to2  .  n2 

dy    dy       Ay 


(3.60) 


+  —  Ey  +  \y2  +  B\  g(y)  =  0  ,  (3.61) 


A  +  B  =  n 


(3.62) 


Field-Free  Solution 

For  the  field-free  case  we  take  A  =  0  and  obtain  the  solutions  (in  Z  =  1 
atomic  units) 


with 


£o  = 


2n2 


Aq  =  ni  + 


m  +  1 


2 

BQ  =  n2  +  — - —  , 


(3.63) 


n  =  ni+n2  +  m  +  l 


(3.64) 


Here  n\  and  n2  are  the  integer  valued  parabolic  quantum  numbers.61  The 
field- free  eigenfunctions  are 


fo(x)  =  u(ni,m,x) 
90(y)  =  u(n2,m,y) 


where 


u(ni,m,x)  = 


"nix 


}xm/2 


,2L™{x) 


[(ni  +  m)!]3/2 
L™iix)  is  an  associated  Laguerre  polynomial  given  by63 


(3.65) 


(3.66) 
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»i 


££(*)  =  [(ni+m)!]2£ 


(-»)» 


p=0 


(m  +  p)!  (ni  —  p)\p\ 


(3.67) 


For  our  case  m  is  restricted  to  integers  for  which  m  >  0.  Other  definitions  of 
L™(x)  are  quite  common. 


Field  Dependent  Solution 

In  order  to  solve  Eqs.  (3.60)  and  (3.61)  for  A  ^  0  we  will  expand  the 
unknowns  f(x),  g(y),  A,  B  and  E  in  powers  of  A.  The  zero  order  terms  will 
be  given  by  the  A  =  0  results  above.  This  perturbation  expansion  technique 
closely  follows  the  work  of  Hoe  et  al.     '"^  The  expansions  are 


1  +  4  Yl  up\P 
P>1 

a  =  y,apxp 

p>0 

B  =  J2  BPXP  (3-68) 

p>0 

f(x)  =  V,  Ap    2_]    apu(nl  +  h,m,x) 
p>0       h^-riy 

p>0       h>—ri2 
The  sums  for  index  h  are  over  all  possible  values  of  n\  and  ri2-  Upon  evaluation 

of  the  coefficients  al  and  bp,  however,  we  will  find  that  all  terms  in  the  sum 

are  zero  except  for  max(— n\,—  2p)  <  h  <  2p.   This  is  also  true  for  the  sum 

involving  ni-   We  immediately  see  that  because  of  the  condition  A  +  B  =  n, 

we  have 


A0  +  B0=n 


(3.69) 


and 
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Bp  =  -Ap       for    p  >  0  .  (3.70) 

The  other  coefficients  in  the  expansion  can  be  obtained  by  solving  for  the  p'th 
coefficients  in  terms  of  the  p- 1  and  lower  coefficients.  Thus,  a  knowledge  of  the 
zero  order  terms  will  allow  us  to  construct  all  higher  order  terms  iteratively.  To 
accomplish  this  we  substitute  the  expansions  for  f(x)  and  g(y)  into  Eqs.  (3.60) 
and  (3.61),  multiply  by  u(x)  and  u(y)  and  integrate  over  x  and  y,  respectively. 
For  the  x  coordinate  we  obtain 

/•OO 

/      dxu(nl  +  k,m,x){H0(x)  +  A  +  ^(l  +  2n2E)-\x2\f(x)  =  0  ,   (3.71) 
where 

tt  /   \        d      d       mr      x 

*w  -  t*xtx  ~  H  -  3  •  <3-72> 

Following  Hoe  et  al64'65  we  call  Hq(x)  the  field-free  effective  Hamiltonian  for 
the  x  coordinate.  It  gives 

H0(x)u(ni  +  h,  m,  x)  =  -A$u(n\  +  h,  m,  x)  ,  (3.73) 

with 


m  +  1  . 

Aq  =  ni  +  h  +  — ^—  .  (3.74) 

Carrying  out  the  integration  gives 


2P 


/  ^P  /]  aP\  dxu(ni +k,m,x)HQ(x)u(ni  + h,m,x) 

P>0      max(-ni-2p)        ^J°  (375) 

+  A6kth  +  (1/4)(1  +  2n2E)  {x)k>h  -  X  (z\h\  =  0  , 


where 


and 
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/•OO 

(xTl)k  h  —   I      dxxnu(n\  +  k,m,x)u(ni  +  h,m,x)  (3.76) 

Jo 


/•oo 

/      dxu(n\  +  k,m,x)u(n\  +  h,m,x)  =  6^  /,  .  (3.77) 

JO 

Upon  expanding  E  and  A  in  powers  of  A,  we  see  that  the  p'th  coefficient  is 
given  by 


p-1  2p  (  p-l 

kaP  =  E  AP~i  ai  ~  E  \  (x2)*,/»ap-l  +  ^k,h  E  "p-i^i 

i=0  h>max(-ni  ,-2p)   {  i=0 

(3.78) 

If  the  maximum  value  of  the  sum  index  is  less  than  the  starting  value  of  the 

index,  then  the  sum  is  defined  to  be  zero.  We  will  return  to  the  special  case  of 

k  —  0  below. 

Likewise,  we  obtain  for  the  y  coordinate 


P-l  2p  r  p-l  ] 

kbP  =  E  BP~i  bi+  E  {  (v\hbhp-i  ~  (y)t,h  E  vP-$     ■ 

i-0  h>max(-n2-2p)   [  »'=0  J 

(3.79) 

For  p  =  0  we  must  have 


a0  =  bQ  =  60/l 


(3.80) 


To  obtain  the  energy  coefficients  up  we  take  the  k  =  0  case  of  Eqs.  (3.78) 
and  (3.79).  We  obtain  two  equations  for  Ap.  Equating  these  and  solving  for 
up  gives 


GO 


»'=1 


p-l 


+  Tn  E 1  (y2)o,fc^i  -  (*2)o,/^-i  -  £  *-<  (*)o,ft«?  +  (V)0,hti 


i=l 


>    . 


(3.81) 


The  sum  over  h  is  indexed  from  max(-n\,  -2p)  or  max(-ri2,  -2p)  to  2p  as 
appropriate. 

We  still  need  to  evaluate  the  coefficients  a*  and  b^  for  fc  =  0.  To  accomplish 


this  note  that 


/•oo 

/      dxf2(x)  =  l, 
JO 


(3.82) 


and 


yoo 

/o 


<fyy  (y)  =  i 


(3.83) 


Evaluating  these  by  using  the  series  expansions  for  f(x)  and  #(y)  and  exam- 
ining the  coefficients  of  a  particular  term  of  order  Ap  will  give 


(C  =  — — 


2P  p-l 

2      E      £4«k. 

h>max{— ni,— 2p)  J=l 


(3.84) 


and 


p         2 


2p  p-l 

E      E  hH-i  •  (3-85) 

h>max(— ri2,— 2p)  i=l 
We  now  have  the  coefficients  for  calculating  the  energy  E  and  the  components 

of  the  wavefunctions  f(x)  and  g(y).  To  evaluate  these  we  need  the  integrals 
(x)k,h  and  (x2)k,h  as  well  as  those  for  (y)k)h  and  {y2)k,h-  These  can  be  straight- 
forwardly evaluated  to  give 


CI 


(x)M  =  (2r»i+2*  +  m  +  l)*M 

-  y/(ni  +  k  +  l)(r»i  +k  +  m  +  l)6k+lfh  (3.86) 

-  y/(ni  +  k)(ni  +  k  +  m)6k_lh, 


and 


(*\h  =  (6(«1  +  k)(nY  +  k  +  m  +  l)  +  (m  +  l)(m  +  2)}  Skfh 

-  2y/(ni  +  k  +  l)(ni  +  k  +  m  +  l)(2ni  +  2fc  +  m  +  2)6k+1>h 

-  2y/(n\  +  k)(n\  +  k  +  m)(2n\  +2k  +  m)8k_lfl 

+  y/(nl  +k  +  2)(r»i  +  k  +  m  +  l)(n\  +  k  +  l)(ni  +  fc  +  m  +  2)6k+2th 

+  y/(nl  +  ^)(nl  +  k  +  m-  l)(ni  +  fc  +  m)(ni  +  k  -  l)^_2,/i  • 

(3.87) 

The  relations  for  (y)j.  /j  and  (y2)^  /,  are  obtained  from  these  formulae  by  taking 

x  — >  y  and  nj  — ►  n2- 

Normalization 

In  order  to  compute  matrix  elements  from  the  field  dependent  wave  func- 
tions given  in  Eq.  (3.58),  we  also  need  to  determine  their  normalization.  We 
take  the  normalization  constant  to  be 


n3    roc    yoo 
N(ni,n2,m,\)  =  —  /      dxdy  (x  +  y)[f(x)g(y)]2  .  (3.88) 

4  Jo    Jo 

Expanding  N(ni ,  n,2,  m,  A)  in  powers  of  A  we  have 

N(ni ,  n2,  m,  A)  =  ^  iVpAp  .  (3.89) 

p>0 

To  evaluate  the  coefficients  Np,  we  expand  the  functions  f(x)  and  y(y)  in 
powers  of  A,  matching  terms  on  the  right-  and  left-hand  sides  to  obtain 


G2 


3    P 


i=0  /»,Jb 


(3.90) 


The  second  sum  is  indexed  from  max(— ni,—  2p)  or  ma:r(— ^2,  —  2p)  to  2p  as 
appropriate.  Note  that 


iV0  =  \nA  . 


(3.91) 


Dipole  Matrix  Elements 

Now  that  we  have  developed  the  formulae  for  the  perturbation  solution  of 
the  Schrodinger  equation  for  a  radiator  in  the  presence  of  a  uniform  electric 
field,  we  can  go  on  to  calculate  the  dipole  and  quadrupole  matrix  elements 
needed  for  the  calculation  of  the  spectral  line  shape  I(u>).  Since  we  know 
the  wave  functions  only  in  terms  of  their  perturbation  expansions,  it  will  be 
necessary  to  develop  perturbation  expansions  for  the  matrix  elements  as  well. 

We  will  follow  the  calculation  of  the  atomic  matrix  element  of  the  coor- 
dinate z  as  an  example.  The  matrix  elements  for  x  and  y,  as  well  as  the 
quadrupole  matrix  elements  are  calculated  in  the  same  way  and  only  the  re- 
sults are  presented  here.  Further  details  concerning  the  technique  can  be  found 
in  the  references  of  Hoe  et  al.    '"" 

We  wish  to  calculate  the  atomic  matrix  element 

(t/>(n\ ,  n'2,  m',  A')|z|0(ni,  n2,  m,  A))  . 
The  normalized  wave  function  is  given  by 


if>  =  %l>(ni,n2,m,\)  = 


y/2TrN(ni,n2,m,\) 


fn1m\(t)9n,m\{l)JmV  .      (3-92) 
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so  the  matrix  element  of  the  z  coordinate  of  the  radiator's  bound  electron  is 

{^\^\^)=lJd3xU-r,)r(n'hn'2,m',X')^(nhn2,m,\) 
1      roo     />oo     r2w 

•J  U       «/  U       «/  u 

TOO       /-OO 
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x  fn'1m'X'U/n')fnlmX((/n)gn>m/x>(ri/n')gn2mX(i1/n)  , 


w 


here 


In  order  to  evaluate  $,  consider 


We  have  defined 


(3.93) 


$  =  yjN(n'ln'2m'X')N(nin2m\)  .  (3.94) 

As  a  notational  simplification,  take 


ZJS3S  =  (0'W)  .  (3.95) 


A/"  =  /Y(ni ,  n'2,  m',  A')i\r(ni ,  n2,  m,  A) 

=  E  E  Q'P'QvFp'+p 

p'>0p>0  (3.96) 

=  E  E  *<?p-.-^ 

p>0  »=0 

=  JTMPFP. 

p>0 


3\  P 


G4 


and 


^sEWh-  (3-98) 

fc=0 

This  gives  us  an  expansion  for  Af  in  powers  of  the  field  F.  Next,  we  use  this 

to  obtain  the  expansion  of  $.  Since  $  =  y/Jf,  we  have 

$2  =  M"  ■  (3.99) 

If  we  expand  $  and  A/"  in  powers  of  A  and  multiply  it  out,  we  obtain 


EE*«v^  =  Ev 

P>0  i=0  p>0 


So, 


Note  that  p  =  0  gives 


»=0 

P-1 


$o  =  Vtfi  =  Jl%No 


Solving  for  $p  gives 


(3.100) 


(3.101) 


(3.102) 


*P  = 


Xp  -  Ed  *i*p-i 

2^ATP 
or  in  terms  of  the  more  basic  quantities 


for    p  >  1  , 


(3.103) 


$n  = 


n'  nz 


p-1 


5  E("/3),'(n3)p-,>'«Ap-i  -  £  MH     .  (3-104) 

*'=0  ,=i 
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for  p  >  1.  For  the  p  =  0  case,  we  have 

$0  =  9n'2"2  •  (3.105) 

To  proceed  with  the  evaluation  of  Zn\n$n  ,  return  to  Eq.  (3.93).  We  can 
separate  the  functions  of  the  two  coordinates  to  obtain 

*%8Z  =  %r  {f(2)g(0)  "  f(0)g(2)}  •  (3-106) 

where 

fM  =   rdttafn,im,x,(t/n')fnimX(Z/n)  ,  (3.107) 

J  u 

and 

G(°)  =   /      driTiOg^X'il/^gnimxWn)  •  (3-108) 

The  orthogonality  condition  will  not  apply  since  we  are  dealing  with  integrals 
of  the  form  jj°  dtf({/n')f(£/n)  and  in  general  n'  /  n. 

To  evaluate  Eqs.    (3.107)  and  (3.108),  we  expand  their  left-hand  sides  in 
powers  of  F  to  obtain 


F(a)  =  J2F^FP 


(3.109) 


P>0 


Expanding  the  right-hand  sides  of  the  integrals  and  rearranging  as  a  series  in 
increasing  powers  of  F  allows  us  to  equate  terms  of  like  powers  in  F  from  both 
sides  to  obtain 


GG 


h,k  i=0 

h,k  1=0 


(3.110) 


where  the  sums  over  /i  and  fc  are  as  before,  and  where 


Tkah=  I      ^^(«l  +  A;,m^/n)U(n'1  +  /l,m,,e/n'),  (3.111) 

j  u 


and 


(  \         t 
Jk°h  =  drl  r)au(n2  +  k,  m,  r]/n)  u(n'2  +  h,  m' ,  77/n')  .  (3.111) 

•/  u 

^fc,A  can  ^e  evamated  from  the  definition  of  u(ni,m,r}/n)  to  obtain 


r(/)       /  2nn'  Y+1     r 

%h  =  (^^^7  J      v(ni  +  fcWni  +  k  +  m)-'K  +  ^K  +  ft  +  m)! 

/    2n>    yf    2n    \4  »f*  ^  /  -2„'  \f  /  -2n  V 

(m  +  q)\{ni  +k-  q)lq\(ml  +  6)!(n'j  +h-  s)ls\  * 

For  JJ£,  just  take  n\  ->  n2  in  Eq.  (3.112). 

To  evaluate  Zn,1^2™   we  use  the  expansions  for  F(a)  and  G(a\  We  obtain 


with 


G7 


'■nin2m    —  2_^  1Vr 
p>0 


where 


We  can  also  introduce  the  expansion 


=  Ezpfp 

p>0 


Rewriting  Eq.  (3.113)  as 


pn\n2m 
■nin2m    — 


*SS 


n\n„m 


1"2 

n2m    i 


and  using  the  expansions  for  each  side  gives 


Tp  =  £  Zf$p_i 

z'=0 


Solving  this  for  Zp  gives 


(3.114) 


*'fc*t(>P<%*->P<&)  .  (3.115) 

i=0 


(3.116) 


(3.117) 


(3.118) 


Zp  - 


(3.119) 


If  p  =  0  we  take  the  sum  to  be  zero.  This  expression  gives  the  terms  in  the 
expansion  for  (xp'\  z  \ip)  given  by  Eq.  (3.116).  We  have  therefore  successfully 
calculated  the  perturbation  expansion  for  the  atomic  matrix  element  of  the 
radiator's  z  coordinate. 

We  have  followed  the  calculation  of  the  z  matrix  element  in  detail  in  order 
to  demonstrate  the  perturbation  theory  method.  For  the  remaining  matrix 
elements  we  present  only  the  final  results.    The  matrix  elements  for  x  and  y 
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are  calculated  by  the  same  procedure  except  for  the  angular  integrals  which 
we  evaluate  explicitly.  For  x  =  V^/  cos  tfi  we  have 


f2*  0 

/      dip  cos ipe-<m' ~m^  =  ~\8   j,     ,+8      i     A 
JO  2    L  m+!>m   ^  °m-l,m'\ 


2 

2       m,m' 


(3.120) 


For  y  —  y/£ij  sin  if  we  have 


fdv  sin^-l™'-*  =  |  [6m+1<m,  -  Sm_iM] 


?2      m,m' 


(3.121) 


The  factor  of  27r  will  cancel  with  a  similar  factor  from  the  wave  function  nor- 
malization. 

The  x  matrix  element  is  given  by 


y"in2™'  -X^  Y    jrP 


where 


p>0 


rP-l 


x     RP  -  Ego  x&p-i 

p  ~ ^ 

$0 


with 


The  y  matrix  element  is  given  by 


*EYPFP> 

p>0 


where 


(3.122) 


(3.123) 


R*  -  X-  E  (*f  ^  +  *?*><£/?>)  .  (3.124) 

i=0 


(3.125) 
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,P-l 


with 


1> ^—  (3-126) 


^  -  §<L  E  (^(3/2)^i)  +  F^G^)   •  (3.127) 

Quadrupole  Matrix  Elements 

The  last  atomic  matrix  element  we  need  is  the  quadrupole  moment  tensor 
component  Qzz.  Its  operator  in  terms  of  radiator  electron  coordinates  is  given 
by  (taking  e  =  1) 


Qzz  =  _  (3*2  -  r2) 

The  matrix  elements  of  Qzz  are  then  given  by 


(3.128) 


where 


(Qzz)n\nlm=J2QpFP'  (3-129) 

p>0 


Q  _sP-ZUQi*P-i  mqm 

Qp  -  Q- .  (3.130) 


with 


i=0 

(3.131) 

The  calculation  of  Qzz  goes  beyond  the  work  of  Hoe  et  al64>65  but  the  method 
is  the  same  and  is  as  outlined  above  in  the  calculation  of  the  z  matrix  element. 
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For  the  perturbation  theory  calculations  in  this  work,  we  will  go  to  sixth 
order,  taking  p  =  6.  This  will  be  more  than  sufficient  for  our  purposes. 

3.3  Atomic  Data  by  Numerical  Solution 
Next,  we  give  a  qualitative  discussion  of  the  numerical  solution  of  Schrodinger's 
equation  for  a  hydrogenic  radiator  in  a  uniform  electric  field.  We  have  already 
developed  the  perturbation  solution  useful  for  relatively  small  field  values.  Di- 
rect numerical  solution  is  applicable,  in  principal,  for  fields  of  any  magnitude, 
though,  in  practice,  it  turns  out  to  be  inconvenient  for  very  small  field  val- 
ues because  it  is  difficult  to  numerically  handle  the  resultant  extremely  sharp 
resonances.  We  will  discuss  this  point  further  below.  Consequently,  the  pertur- 
bation and  numerical  solution  techniques  are  complementary  and  additionally 
should  serve  as  consistency  checks  since  they  should  produce  matching  results 
for  a  suitable  intermediate  range  of  fields.  For  large  field  values,  the  resonance 
nature  of  the  radiator  states  becomes  an  essential  part  of  their  description  and 
the  perturbation  theory  used  here  is  no  longer  useful  for  finding  the  field  depen- 
dent eigenvalues  and  eigenfunctions  of  Schrodinger's  equation.  We  also  discuss 
the  calculation  of  the  width  of  these  resonances  and  the  numerical  evaluation 
of  the  necessary  atomic  matrix  elements. 

The  presence  of  the  uniform  field  at  the  radiator  changes  the  very  nature 
of  bound  atomic  states  by  turning  each  discrete  energy  level  into  a  shape  reso- 
nance; in  other  words,  there  will  be  a  solution  for  a  particular  set  of  quantum 
numbers  that  includes  a  continuum  of  energy  values.66'67  For  a  field  strength 
given  by  A,  the  states  can  be  characterised  uniquely  by  the  quantum  numbers 
m/,  A  and  the  energy  E.  The  quantum  numbers  m/  and  A  remain  discrete 
but  E  now  has  a  continuous  spectrum.    For  weak  field  values,  the  density  of 


71 
states  is  greatest  for  a  radiator  electron  with  an  energy  that  closely  corre- 
sponds to  the  discrete  states  given  by  the  perturbation  theory  solution  to  the 
problem.    As  the  field  value  increases,  however,  the  density  of  states  broad- 
ens and  the  probability  of  finding  the  electron  with  an  energy  more  widely 
spaced  from  the  resonance  center  increases.    This  is  a  result  of  the  field  in- 
duced broadening  of  the  resonance.    The  width,  T,  of  the  resonance  is  also 
roughly  inversely  proportional  to  the  lifetime  of  the  electron  inside  the  radia- 
tor.  This  is  most  easily  seen  from  the  time-energy  uncertainty  principle  but 
also  follows  from  a  detailed  WKB  treatment  of  resonance  decay.68  In  fact,  the 
WKB  treatment68'69  of  the  field  ionization  problem  leads  to  a  probability  of 
finding  the  quasibound  electron  inside  the  potential  barrier  given  by  the  time 
dependence  factor  e~rt.    Thus,  when  t  =  1/T  the  probability  has  decreased 
from  one  at  t  =  0  to  1/e  =  0.37.  The  presence  of  the  field  lowers  the  Coulomb 
potential  along  the  upfield  direction  thus  producing  a  potential  well  defined  by 
a  finite-size  potential  barrier.    The  electron  is  then  able  to  tunnel  out  of  the 
radiator.    Indeed,  in  the  presence  of  the  field,  there  are  now  no  truly  bound 
states  in  the  sense  of  the  electron  being  permanently  associated  with  a  partic- 
ular atom  baring  radiative  or  collisional  ionization.  However,  if  the  lifetime  of 
the  electron  inside  the  radiator  potential  barrier  is  long  compared  to  the  time 
of  interest  for  the  radiator,  it  is  effectively  bound  to  the  radiator.    For  even 
larger  field  values,  the  plasma  electric  field  potential  will  exceed  the  attractive 
central  Coulomb  potential;  the  energy  level  will  be  above  the  top  of  the  poten- 
tial barrier.   In  this  case,  if  the  energy  level  is  close  to  the  top  of  the  barrier, 
the  electron  can  still  spend  a  significant  amount  of  time  in  the  vicinity  of  the 
radiator  nucleus  and  will  retain  some  of  the  character  of  a  bound  state.  As  the 
field  is  increased  further  this  bound  character  will  gradually  be  reduced.  This 
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phenomenon  results  in  the  smooth  broadening  of  the  resonances  as  the  field 
increases  until  they  overlap.  At  that  point  the  resonances  have  merged  into  a 
relatively  smooth  background  continuum. 

We  need  to  numerically  solve  the  system  of  coupled  second  order  ordinary 
differential  equations  (ODE's)  given  by  Eqs.  (3.60)  and  (3.61).  In  order  to  em- 
ploy the  solutions  of  these  equations  for  the  calculation  of  the  matrix  elements 
and  line  shapes,  and  since  the  eigenvalues  are  not  discrete,  we  will  use  the 
Breit-Wigner62  Lorentzian  parameterization  of  resonances  and  take  the  eigen- 
functions  and  eigenvalues  at  the  center  of  each  resonance  as  a  representative 
value  over  the  entire  resonance.  As  long  as  the  resonance  is  fairly  narrow  and 
remains  distinct,  this  is  a  reasonable  approximation.  It  breaks  down  for  strong 
fields  where  the  resonance  becomes  strongly  asymmetric.   Fig.   (3.3)  gives  an 
example  of  this  phenomenon.  The  reason  for  this  asymmetry  lies  in  the  shape 
of  the  potential  barrier;  as  the  energy  level  goes  higher  it  sees  a  thiner  potential 
barrier.  The  thiner  the  potential  barrier,  the  shorter  the  lifetime  of  the  state 
and  the  greater  the  uncertainty  in  its  energy.  Consequently,  the  broader  width 
on  the  high  energy  side  of  the  resonance  is  due  to  the  greater  uncertainty  in 
the  energy  value  at  that  point.  These  resonance  asymmetries  have  been  inves- 
tigated experimentally  by  Harmin.70  We  do  not  need  to  describe  this  behavior 
in  detail  since  its  contribution  to  the  line  shape  is  greatly  attenuated  due  to 
the  low  probability  of  occurrence  of  the  relevant  high  field  values  as  reflected 
in  a  small  value  of  the  microfield  probability  function. 

We  employ  the  numerical  solutions  of  Eqs.  (3.60)  and  (3.61)  by  R.  Mancini 
and  C.  Hooper71  who  follow  the  method  of  E.  Luc-Koenig  and  A.  Bachelier.66'67 
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We  give  here  a  general  qualitative  procedural  description  of  the  solution  tech- 
nique (see  the  above  references  for  further  details).  The  solutions  for  the  equa- 
tion in  the  scaled  variable  x  are  effectively  bound  and,  therefore  for  large  x 
are  exponentially  damped.  The  number  of  nodes  in  the  wavefunction's  x  com- 
ponent is  given  by  the  parabolic  quantum  number  n\.  The  solutions  for  the 
equation  in  the  scaled  variable  y  are  effectively  unbound;  that  is,  the  solutions 
for  large  y  outside  the  potential  well  are  oscillatory  in  nature.  On  the  other 
hand,  for  x  and  y  ->  0  we  have  /(0)  =  #(0)  =  0. 

To  solve  the  the  two  equations,  Eq.  (3.60)  and  Eq.  (3.61)  consistently, 
first,  we  solve  the  effectively  bound  equation  in  x  for  a  particular  case  given 
by  specific  values  of  F,  E,  nj  and  m,.  This  will  give  a  value  for  the  constant 
A  by  imposing  the  exponentially  damped  behavior  of  F(x)  for  large  x  (this  is 
an  eigenvalue  problem  for  A).  We  use  this  to  determine  the  constant  B  and 
proceed  with  the  solution  of  the  effectively  unbound  equation  in  y.  We  change 
Eqs.  (3.60)  and  (3.61)  into  a  form  more  convenient  for  numerical  solution  by 
eliminating  the  first  order  derivative  with  the  transformations 


F{x)  =  y/xf(x)  , 


(3.132) 


and 


G(v)  =  \/yg(y) 

Our  two  equations  will  then  be  transform  to 


(3.133) 


d2T{x) 
dx^ 


+  TxF{x)  =  0  , 


and 


(3.134) 


where 


and 
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d2G(y)  ,Trf,     n 


(3.135) 


(3.136) 


B      n 


+  —E- 


m' 


+  Ay 


(3.137) 


y        2  "  4y2 

The  motivation  for  making  this  transformation  is  that  now  we  have  to  deal 
with  second  order  ODE's  with  missing  first  derivatives  which  can  be  efficiently 
integrated  using  the  Numerov  algorithm. 

To  solve  Eq.  (3.134)  for  an  energy  value  E,  a  scaled  field  value  A,  a  value 
for  the  z  component  of  the  orbital  angular  momentum  m/  and  a  value  for 
the  parabolic  quantum  number  r»i,  we  pick  a  trial  value  for  A.   We  solve  the 
differential  equation  starting  at  x  =  0  using  a  power  series  expansion;  this 
is  used  to  initialize  an  outward  numerical  integration  using  a  fourth  order 
Numerov  algorithm.66-72  This  numerical  integration  is  continued  up  to  the 
outer  turning  point  of  the  effective  potential  barrier.    (If  the  energy  level  is 
above  the  top  of  the  potential  barrier  so  that  there  is  no  outward  turning 
point,  we  use  the  maximum  of  the  potential  barrier.)    For  suitable  large  x 
values,  an  exponentially  damped  asymptotic  solution  is  used  to  begin  an  inward 
numerical  integration,  again  using  the  Numerov  method.    The  two  solutions 
overlap  at  the  outer  turning  point  of  the  potential  barrier  (or  the  potential 
barrier  maximum  if  the  energy  level  is  above  it).    At  the  meeting  point,  we 
match  F(x)  by  adjusting  a  multiplicative  constant,  and  from  the  matching 
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condition  for  djr(x)/dx  we  compute  a  correction  to  our  initial  guess  for  A. 
This  defines  an  iterative  procedure  which  is  continued  until  the  initial  guess 
for  A  agrees  within  a  given  tolerance  with  the  value  found  from  matching  the 
derivatives.  * 

To  solve  Eq.  (3.135),  the  unbound  equation  in  the  scaled  variable  y,  we 
note  that  once  we  have  the  value  of  the  constant  A  we  have  B  —  n  —  A.  This 
gives  us  E,  A,  n\,  mi,  A  and  B;  therefore,  it  is  not  necessary  to  make  any  initial 
guesses.  We  start  the  outward  integration  as  we  have  done  for  the  variable  x 
using  a  power  series  solution,  and  then  switch  to  the  numerical  integration 
using  the  Numerov  method  as  before.  Since  the  constant  B  has  already  been 
determined,  we  continue  the  integration  outward  to  large  y  beyond  the  effective 
potential  barrier  where  we  fit  a  large-y  analytic  oscillatory  asymptotic  solution 
to  the  numerical  results.  This  will  give  the  value  of  the  amplitude  of  Q(y)  for 
large  y.  By  studying  the  distribution  of  amplitudes  with  respect  to  the  energy, 
the  width  and  center-location  of  the  resonance  can  be  estimated. 

Since  the  values  of  the  energy  can  have  a  continuous  spectrum,  we  have 
the  normalization  condition  per  unit  energy  given  by 


r27T 


1       rlTT  /-00  />oo 

4  Jo     dipJo      Vo    d*ti  +  ri^rnt,A,E*rt,A'tE'  =  6m,!*lsA^6(E-E'). 

(3.138) 
it  can  be  shown62  this  is  equivalent  to  the  condition  Iout  =  1/27T,  where  Iout 
is  the  outward  flux  in  the  77  direction  at  large  rj.  The  outward  flux  Iout  can 
be  given  in  terms  of  the  outgoing  wave  if>out  where  we  have  decomposed  the 
stationary  wavefunctions  into  outgoing  and  incoming  components  in  the  77  di- 
rection. 
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The  density  of  states  can  be  obtained  by  studying  the  behavior  of  the 
electron  probability  density  near  the  nucleus  for  the  component  of  the  wave- 
function  associated  with  a  constant  value  of  the  outward  flux.  As  stated  before, 
we  can  approximate  the  density  of  states  for  a  particular  resonance  using  the 
Breit-Wigner  parameterization.      We  obtain 

D»(£)  =  !%(*)— J*-p-  ,  (3.139) 

where  Er  corresponds  to  the  energy  value  at  the  maximum  of  the  density  of 
states  of  the  resonance,  and  T  is  the  full  resonance  width  at  half  maximum. 
This  is  a  reasonable  parameterization  as  long  as  the  resonance  remains  distinct. 
The  average  lifetime  Tres,  of  the  resonance  can  be  associated  with  the  resonance 
width  through65 


Tres  ~  '•/* 


(3.140) 


This  is  just  a  manifestation  of  the  time-energy  uncertainty  principal. 

Once  we  have  numerically  evaluated  F(x)  and  Q(y),  we  are  free  to  calculate 
the  matrix  elements  needed  for  the  line  shape.  By  straight  forward  numerical 
integration  we  evaluate  matrix  elements  for  x,  z  and  Qzz.  We  cut-off  the  in- 
tegration over  the  scaled  variable  y  at  the  outer  turning  point  of  the  effective 
potential  barrier  for  energy  values  below  the  maximum  of  the  potential  bar- 
rier. (For  energies  greater  then  the  maximum  of  the  potential  barrier,  we  stop 
the  integration  at  the  value  of  y  corresponding  to  the  maximum  of  the  poten- 
tial barrier.)  Due  to  their  resonance  nature,  the  wavefunctions  extend  overall 
space  but  we  are  only  interested  in  the  portion  located  near  the  radiator.  The 
numerical  values  of  the  matrix  elements  are  not  significantly  sensitive  to  small 
changes  in  the  location  of  this  integration  cut-off.  Since  we  now  have  properly 
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normalized  wavefunctions,  we  can  evaluate  the  field  dependent  matrix  element 
of  a  real  operator  O  between  two  resonance  states  denoted  by  subscripts  "i" 
and  "j"  by  integrating  over  the  energy  range  of  each  resonance.  The  square  of 
the  matrix  element  for  an  operator  O  can  be  written  as, 

\Oij(F)\2  =  J dE{  J dEj  \{MEi,F)\G\^{E^F))\2  .  (3.141) 

For  narrow  resonances  we  can  approximate  this  results  by 


Oij 


.    1    .        .   1 

1    ^i(Erti,F)\0\^(ErJlF))  , 


(3.142) 
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where  we  have  used  the  value  of  the  eigenfunctions  at  resonance  center.  This 
is  a  reasonable  approximation  because  the  main  field  dependence  of  the  wave 
functions,  for  narrow  resonances,  is  contained  in  the  normalization  factor  which 
is  a  function  of  the  square  root  of  the  density  of  states.  Our  evaluation  of  the 
square  of  the  matrix  element  allows  us  to  integrate  over  the  density  of  states 
to  obtain  the  7iT,y2  factors.   Since  these  factors  are  always  positive  and  real, 
the  sign  (and  any  factors  of  t)  of  Oij  is  reliably  given  by  the  evaluation  of 
(^i(Ei,F)\0\^j(Ej,F))  at  Ei  =  Er%i.     Fig.     3.4  gives  an  example  of  the 
calculation  of  atomic  matrix  elements  by  both  perturbation  theory  and  direct 
numerical  solution  of  Schrodinger's  equation.    We  see  that  the  two  solutions 
join  smoothly  for  intermediate  field  values  and  diverge  for  larger  field  values 
as  the  accuracy  of  the  perturbation  theory  solution  breaks  down. 

As  the  field  F  goes  to  zero,  Eq.   (3.141)  will  reduce  to  the  usual  no  field 
limit  where 


Oij(0)  =  U 


0 


4) 


(3.143) 
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Here, 


+i  =  u(nhmh^)u(n2,mi,r]) 


0tm,ip 


2tt 


(3.144) 


which,  as  expected,  is  just  the  solution  to  Schrodinger's  equation  for  the  field 
free  case.  Here,  the  state  "i"  is  specified  by  E,  n\,  mi  and  F.  Since  for 
a  solution  with  a  discrete  values  of  E  there  is  a  unique  value  of  the  principal 
quantum  number  n,  we  can  obtain  n2  through  the  relation  n  =  ni+n2  +  \mi\  +  l 
and  our  limiting  state  can  be  characterized  by  n\,  n2  and  m/  as  expected. 

3.4  Electron  Derealization  and  Field  Ionization 
We  have  been  using,  as  a  model  for  the  radiator-plasma  interaction,  the 
picture  of  a  radiator  in  the  presence  of  a  uniform  plasma  microfield  produced 
by  the  averaging  of  the  fields  from  all  the  plasma  particles.    We  have  added 
the  ion-quadrupole  correction  to  account  for  a  relatively  small  field  gradient  at 
the  radiator.   This  picture  is  only  a  model,  of  course.    As  the  plasma  density 
increases,  it  becomes  unreasonable  to  ignore  the  presence  of  highly  charged 
nuclei  close  to  the  radiator.  This  fact  creates  a  qualitatively  different  picture. 
In  the  uniform  field  case,  a  quasibound  electron  will  tunnel  out  of  the  radiator 
and  continue  to  accelerate  toward  the  field  source  at  infinity.  In  reality,  large 
fields  are  produced  primarily  by  nearby  charges.     The  quasibound  electron 
will  then  orbit  about  the  ensuing  multicenter  ionic  potential  in  a  molecular 
orbital      instead  of  tunneling  into  the  interparticle  space;  the  electrons  will 
become  delocalized.    A  complication  to  this  molecular  orbital  picture  is  the 
large  number  of  free  plasma  electrons  in  the  internuclear  volume.  It  is  unrea- 
sonable to  assume  that  the  quasibound  electrons  will  not  strongly  interact  with 
these  free  electrons  as  they  resonate  between  various  nuclei.  Perhaps  a  more 
accurate  picture  for  these  conditions  would  be  the  self  consistent  treatment  of 
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both  bound  and  free  electrons  in  the  presence  of  the  multicenter  potentials. 
Some  work  has  been  done  along  these  lines  from  an  analytic  point  of  view  by 
Rogers,75  and  from  a  numerical  simulation  point  of  view  by  Younger  et  al.7" 
Additionally,  a  model  consisting  of  a  cluster  of  nuclei  with  bound  electrons  has 
been  considered  by  Collins  and  Merts77  but  for  Te  =  0.  Much  work  remains 
to  be  done  before  these  models  are  capable  of  accurately  treating  spectral  line 
transitions  for  plasmas  of  high  density  and  temperature. 

In  addition  to  the  question  of  whether  field  induced  ionization  or  resonance 
exchange  of  electrons  between  nuclei  is  a  better  description  of  what  strong 
fields  do  to  quasibound  electrons,  there  is  the  question  of  the  level  broadening 
produced  by  these  two  mechanisms.  The  field  ionization  picture  contains  level 
broadening  due  to  the  resonance  character  of  the  energy  levels  in  the  presence  of 
the  nearly  uniform  electric  field.  This  is  most  easily  seen  from  the  uncertainty 
principle  where  the  width  is  related  to  the  resonance  lifetime  by  AE  «  h./Tres. 
In  the  molecular  orbital  picture,  an  electron  bound  to  two  nuclei  will  still  have 
a  discrete  energy.  However,  as  the  number  of  nuclei  in  the  molecule  or  cluster 
increases,  splitting  is  induced.  If  each  level  is  m-fold  degenerate  and  there 
are  N  nuclei  in  the  cluster,  the  splitting  will  be  on  the  order  of  miV-fold.77 
This  large  degree  of  splitting  will  be  qualitatively  similar  to  the  intrinsic  level 
broadening  produced  in  the  field  ionization  case.  We  see,  therefore,  that  the 
level  broadening  from  the  two  different  models  should  produce  similar  results 
under  similar  conditions,  at  least  as  far  as  the  level  width  is  concerned. 

At  this  point,  it  appears  that  the  molecular  orbital  picture  may  be  the 
most  capable  of  describing  our  dense  plasma  since  it  is  a  truly  many-body 
theory.  However,  there  is  the  complicating  factor  of  the  intercollision  of  free 
and  bound  electrons.    As  the  bound  electrons  exchange  between  the  nuclei, 
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they  will  be  subject  to  collisions  with  the  free  plasma  electrons.    If  they  are 
then  thermalized  into  the  plasma,  the  situation  will  be  more  like  the  field 
ionization  case.   A  realistic  picture  probably  lies  somewhere  between  the  two 
over  simplified  models. 

In  the  molecular  orbital  model,  radiator  levels  are  not  depopulated.  On  the 
other  hand,  in  the  field  ionization  model,  all  levels  will  be  depopulated  on  time 
scales  relevant  to  radiative  decay,  for  sufficiently  strong  fields.  If  the  actual 
phenomenon  lies  somewhere  between  these  two  pictures,  then  the  resultant 
effect  will  range  somewhere  between  no  change  in  the  line  shape  due  to  no 
level  depopulation,  to  a  substantial  decrease  in  line  radiation  for  large  fields 
due  to  level  depletion  from  field  ionization.  The  levels  will  not  be  completely 
depopulated,  of  course,  due  to  collisional  and  radiative  repopulation.  What 
effect  could  these  two  models  have  on  the  spectral  line  shape? 

In  order  to  estimate  the  importance  of  this  possible  ionization  effect  on 
level  populations,  we  will  examine  the  limiting  case  of  field  ionization  due  to 
tunneling  with  no  reverse  process  to  repopulate  the  radiator  level  except  the 
usual  collisional  and  radiative  processes.  We  will  construct  a  simple  kinetic 
population  model  to  describe  the  influence  of  field  induced  ionization  on  the 
relative  population  of  the  energy  levels  of  hydrogen-like  argon.  We  will  present 
only  the  extreme  case  of  no  back  tunneling  in  order  to  ascertain  the  maximum 
possible  population  depletion.  As  outlined  above,  the  real  situation  could 
lead  to  results  ranging  somewhere  from  no  tunneling  depletion  at  all  (the  pure 
molecular  orbital  picture)  to  maximum  depletion  due  to  one  way  tunneling  (the 
field  ionization  picture).  We  will  not  look  for  departures  from  absolute  LTE 
populations  but  only  relative  departures  from  LTE  among  the  upper  states 
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involved  in  the  spectral  transitions;  relative  population  differences  are  all  we 
will  be  able  to  observe  from  experimental  line  spectra. 

We  will  examine  the  relative  distribution  of  populations  of  the  n  =  3  and 
4  levels  of  hydrogenic  argon.  The  model  includes  connections  between  the 
hydrogenic  excited  state  levels  to  the  ground  state  and  the  fully  stripped  ion 
while  omitting  the  n  =  2  levels  which  are  little  affected  by  field  ionization.  At 
the  high  densities  we  are  concerned  with  in  this  work,  collisional  processes  will 
dominate  the  other  population  altering  processes  and  relative  LTE  populations 
will  be  very  probable  if  we  do  not  include  the  field  ionization.  We  consider  only 
these  few  levels  because  we  are  mainly  interested  in  the  question  of  whether  the 
field  ionization  effect  can  depopulate  the  levels  enough  to  produce  a  significant 
change  in  the  line  shape. 

We  will  incorporate  into  the  kinetic  model  collisional  excitation,  collisional 
de-excitation,  collisional  ionization,  three  body  recombination,  spontaneous 
emission,  radiative  recombination  and  field  ionization.  We  can  easily  formulate 
a  rate  equation  for  the  population  N{  of  a  level  denoted  by  the  subscript  "i". 
The  number  of  excited  quasibound  levels  is  denoted  by  N,  the  fully  stripped  ion 
by  the  subscript  "N+l"  and  the  ground  state  level  by  "0".  Our  rate  equation 


is 


,N         -1  N+l 

j=Q  j=i+l 

N+l  i-l 

-    £    NineCfj-^NineCfj  (3.145) 

j-i+l  ;=0 

N  +  l  i-l 

+   £  NjAjti  -  J2  NiAij  -  NiTi  , 
j=i+l  ;=0 
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for  i  =  1  to  N  +  1.  Note  that  since  the  fully  stripped  ion  has  no  electron  to 
escape,  I>+1  =  0.  We  will  take  the  populations  to  be  relative  to  the  ground 
state  population,  N0.  This  will  give  Nj  -»  Nj/Nq  and  N0  -»  1  for  j  =  1  to 
N  +  l.  Since  we  are  only  interested  in  the  steady  state  solution  to  these  kinetic 
equations,  we  examine  the  case  with 


dNi 

dt 


=  0, 


(3.146) 


for  all  i.  Next,  we  define  the  rate  coefficients  in  Eq.  (3.145): 
The  collisional  excitation  rate  is  given  by  Cf  •,  where78 

Cfj  =  3.75  x  10-5  T-1/2!  (i\  f  \j)  |2  i-^j/kT  cm3  sec-l  _  (3147) 

Here,  T  is  in  degrees  K,  AEij  is  the  energy  difference  between  levels  i  and  j, 
and  (t|  r  [;')  is  the  atomic  matrix  element  of  r  in  atomic  units. 

The  collisional  de-excitation  rate  is  given  by  Cf  ■.   This  can  be  related  to 
the  collisional  excitation  rate  by  the  principle  of  detail  balance.79  This  gives 


(3.148) 


for  states  of  equal  statistical  weight. 

The  collisional  ionization  rate  is  given  by  C? N  ,  ,,  where80 

Cf>N+1  =  3.46  x  10-5  (j^jp)  T-1/2  e-**  {1  -  x{  e^E^)}  cm3  sec"1  . 

(3.149) 
Here,  T  is  in  degrees  K,  Ij  is  the  ionization  potential  of  level  i  in  eV  and 


Xj  —  111  Kl    , 


(3.150) 
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and 


»oo        „—xt 


f°°      e_a 

El(x)=  dt—  (3.151) 

is  the  first  order  exponential  integral. 

The  three  body  recombination  rate  is  given  by  C^+1  ,.  This  can  be  related 

to  the  collisional  ionization  rate  by  the  principle  of  detail  balance. 

The  spontaneous  emission  rate  is  given  by  the  Einstein  coefficient  A{j, 
where^ 


Aij  =  2.U2xl010(AEiJf\(i\r\j)\2  sec"1   . 

Here,  AEjj  and  («|  r  \j)  are  in  atomic  units. 

The  radiative  recombination  rate  is  given  by  An+\  ;,  where81 


(3.152) 


AN+1>i  =  5.197  x  10-UneZx6i'eXiEi(xi)  sec"1 


(3.153) 


Here,  Z  is  the  radiator  charge,  ne  is  the  plasma  electron  number  density  given 
in  cm   °  and  X{  is  given  by  Eq.  (3.150). 

For  the  steady  state  case  and  because  we  have  taken  Nq  =  1,  we  can  write 
all  of  the  N  +  1  rate  equations  in  the  form 


a0  +  aiJVi  +  a2N2  +  ■■■  +  aNNN  +  aN+1NN+i  =  0 


(3.154) 


where  the  coefficients  a;  do  not  depend  on  the  level  populations.  These  equa- 
tions may  be  cast  in  matrix  form  to  obtain 


C  =  RN  . 


(3.155) 
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Here,  C  is  a  column  vector  of  dimension  N  +  1  containing  the  ag  rate  coeffi- 
cients, N  is  a  column  vector  of  dimension  iV  +  1  containing  the  level  populations 
N{  and  i2isaiV  +  lx./V  +  l  matrix  containing  the  remaining  rate  constants. 
This  equation  is  easily  solved  to  give  the  N  +  I  relative  populations  iVt\  We 
have 


N  =  R~lC 


(3.156) 


This  effectively  solves  the  relative  populations  problem  for  the  N  quasi- 
bound  excited  state  levels  plus  the  the  fully  stripped  ion.  For  the  total  number 
of  quasibound  levels,  we  use  the  n  =  3  and  n  =  4  levels  of  a  hydrogen-like 
ion.  This  gives  a  total  of  N  =  32  +  42  =  25  levels.  The  effect  of  the  lack  of 
inclusion  of  the  n  =  2  levels  was  examined  by  comparison  with  a  more  detailed 
kinetic  population  model.  This  showed  no  significant  change  in  the  relative 
populations  upon  inclusion  of  these  levels. 
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FIGURE  (3.1)  Field  gradient  term  (Ezz)€  as  a  function  of  the  field  for  the 

APEX  ( ),  IP  ( )  and  NN  ( )  models  at  ne  =  1  x  1024  cm-3  and 

kT  =  800  eV.    The  field  gradient  is  in  units  of  e/r3  ^    The  three  models  are 
discussed  at  the  end  of  section  3.1. 
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^  FIGURE  (3.2)  Approach  of  the  constrained  average  field  gradient  term 

{Ezz)e  to  the  nearest  neighbor  limit  for  increasing  density.   ( )  refers  to 

=  1  x  1024  cm"3,  ( )  refers  to  ne  =  5  x  1024  cm"3,  ( )  refers  to 

=  1  x  10     cm   3  and  ( )  refers  to  the  nearest  neighbor  model  discussed 

in  the  text.    All  temperatures  are  kT  =  800  eV.    The  constrained  average  is 
evaluated  using  the  APEX  model  given  by  Eq.  (3.50). 
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FIGURE  (3.3)  Example  of  the  asymmetry  of  a  resonate  state  discussed 
in  section  3.3.    This  is  the  resonance  with  quantum  numbers  n,ni,U2,m  = 

4, 0, 2, 1  or  n,  q,  m  =  4,  —2, 1.  The  resonance  center  is  given  by  ( )  and  the 

width  at  half  maximum  is  given  by  ( ).  The  field  strength  F  is  in  units 

of  a.u./l  x  10       and  the  energy  E  is  in  units  of  a.u./l  x  10      .   a.u.   denotes 
atomic  units  (see  Appendix  A). 
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FIGURE  (3.4)  Comparison  of  numerical  and  perturbation  theory  results 
for  the  calculation  of  a  quadrupole  atomic  matrix  element.  The  example  given 
here  is  a  diagonal  matrix  element  with  n\   =  2,  n2  =  0  and  m  =  0.    The 

numerical  solution  as  discussed  in  section  3.3  is  given  by  ( )  and  the  sixth 

order  perturbation  theory  solution  calculated  from  Eq.     (3.129)  is  given  by 

( )•  Au  quantities  are  expressed  in  atomic  units  scaled  for  radiator  charge, 

Zi  =  1. 


CHAPTER  IV 
RESULTS  AND  DISCUSSION 

In  the  previous  chapters  we  have  formulated  a  general  theory  of  plasma 
spectral  line  broadening  and  discussed  many  of  the  approximations  used  to 
arrive  at  a  calculable  result.  The  theory  that  we  use  to  make  these  calculations 
has  been  generalized  beyond  previous  formulations  in  that  it  includes  higher- 
order  microfield  effects  on:  atomic  matrix  elements,  radiator  energy  levels, 
state  lifetimes,  and  level  populations.  We  will  now  present  and  discuss  the 
results  of  these  calculations. 

The  presence  of  these  higher  order  field  effects  can  be  expected  to  lead  to 
several  discernable  changes  in  the  spectral  line  shapes  of  radiators  in  dense 
plasmas.  For  the  physical  conditions  we  are  examining  the  electron  number 
density,  ne,  varies  from  1  x  1024  to  1  x  1025cm-3  and  at  a  temperature,  kT, 
of  800  eV.  In  this  range,  the  lowest  order  corrections  of  the  ion-quadrupole 
and  the  quadratic  Stark  effects  will  generally  give  rise  to  a  blue  asymmetry  of 
the  spectral  line  shape.  This  means  that  the  intensity  of  the  high  energy  (or 
blue)  side  of  the  spectral  line  will  be  enhanced  over  that  of  the  low  energy  (or 
red)  side  of  the  line.  This  comes  about  by  the  preferential  shifting  to  lower 
energy  of  components  comprising  the  manifold  of  energy  levels  associated  with 
a  principal  quantum  number  n.  If  all  of  the  components  were  shifted  by  an 
equal  amount,  there  would  be  no  discernable  change  in  the  spectral  line  shape; 
only  an  overall  line  shift.  The  n  =  1  ground  state  for  these  transitions  in 
highly  ionized  hydrogenic  ions  is  little  affected  by  the  field  because  of  the 
much  stronger  binding  potential  for  the  n  =  1  state.  Additionally,  the  ground 

89 


90 

state  experiences  no  linear  Stark  effect.  Energy  levels  that  have  been  Stark 
shifted  to  the  low  energy  side  of  the  unperturbed  upper  state  of  the  transition 
correspond  to  quasibound  electrons  that  are  found  to  have  the  maximum  value 
of  their  probability  amplitude  on  the  upheld  side  of  the  radiator  potential 
well.  This  corresponds  to  the  side  of  the  origin  where  the  potential  barrier  has 
a  maximum;  the  other  side,  the  downfield  side,  corresponds  to  the  potential 
barrier  increasing  without  bound.  The  upheld  electrons  are  more  easily  affected 
by  the  field  and  hence,  their  energy  levels  are  shifted  more.  This  preferential 
shifting  produces  a  spreading  out  of  intensity  on  the  low  energy  side  of  the  line, 
and  thus,  an  increase  in  the  peak  height  of  the  high  energy  side  of  the  line. 
The  addition  of  the  ion-quadrupole  effect  enhances  this  trend  but  also  generally 
produces  a  bunching,  or  lessening  of  the  splitting  amongst  the  energy  levels 
on  the  blue  side.  This  further  adds  to  the  blue  asymmetry  of  the  line.  The 
consequences  of  the  field  dependence  in  the  wave  functions  and  their  resulting 
matrix  elements  are  harder  to  characterize  due  to  less  systematic  results  on  the 
line  shape.  Therefore  we  study  this  effect  only  in  combination  with  the  other 
field  effects.  In  this  dissertation,  we  assume  that  the  areas  of  all  line  shapes 
are  normalized  to  one. 

The  possible  field  ionization  depletion  of  the  radiator  upper  level  popula- 
tions will  lead  to  a  lessening  of  intensity  in  the  line  wings.  Again,  the  red  wing 
should  be  more  strongly  affected  due  the  lower  potential  barrier  height  on  the 
upfield  side  of  the  potential  well.  For  sufficiently  strong  fields,  The  Stark  effect 
will  cause  the  energy  levels  from  states  with  adjacent  principal  quantum  num- 
bers to  overlap.  Inclusion  of  this  phenomenon  is  important  for  the  accurate 
representation  of  line  merging. 


91 
4.1  The  Lvman  a  Line 
The  La  calculation  contains  the  effect  of  the  microfield  on  the  atomic 
physics  through  field  dependent  matrix  elements  and  energy  levels.  This  gives 
rise  to  a  field  dependent  fine  structure  correction  to  the  radiator  Hamiltonian 
(see  Appendix  E).  The  field  dependence  appears  in  the  dipole  matrix  elements, 
and  because  we  include  the  ion-quadrupole  effect,  also  in  the  quadrupole  matrix 
elements.  We  do  not  study  the  effect  of  possible  field  ionization  on  the  levels  of 
the  LQ  transition  because,  for  this  case,  the  resonance  width  T  <  0.1  eV.  Since 
for  this  transition  the  resonance  nature  of  the  states  is  not  important,  perturba- 
tion theory  is  adequate  for  the  calculation  of  the  atomic  physics.  At  the  lower 
end  of  the  density-temperature  range  we  are  interested  in  (ne  =  1  x  1024  cm-3, 
kT  =  800  eV),  the  Doppler  effect  is  an  important  source  of  broadening  so 
its  inclusion  is  essential.    For  the  densities  and  temperatures  we  are  dealing 
with,  the  n  =  2  and  3  levels  are  well  separated  in  energy  and  do  not  overlap 
for  any  relevant  field  strengths.    Consequently,  we  will  not  need  to  consider 
this  overlap  for  the  LQ  line  calculation.   In  Fig.   (4.1)  we  have  the  La  line  at 
ne  =  1  x  1024  cm"3  and  kT  =  800  eV.  The  higher  order  field  effects  produce 
only  a  slight  blue  asymmetry  in  the  line  which  attenuates  the  red  peak  by 
less  than  5%  in  magnitude.    The  fine  structure  splitting  is  readily  apparent 
with  Doppler  broadening  being  responsible  for  more  than  half  the  width  of 
the  peaks.    We  estimate  the  Doppler  width  from34  hwD  ~  tiujijv/kT/mrc'2, 
where  fruij  is  the  transition  energy  and  mr  is  the  radiator  mass.  In  this  case 
huiD  ~  3  eV.  As  the  density  goes  up  the  Stark  broadening  will  increase  and  the 
Doppler  broadening  will  become  less  important  as  its  fraction  of  the  total  width 
decreases.  In  Fig.   (4.2)  the  density  has  increased  to  ne  =  5  x  1024cm-3;  the 
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Stark  broadening  has  increased  and  begun  to  obscure  the  fine  structure  split- 
ting. However,  the  higher  order  field  effects  are  still  of  only  minor  importance, 
accounting  for  changes  in  peak  heights  of  less  than  4%.  For  ne  =  1  x  1025  cm-3 
there  is  still  no  significant  attenuation  of  the  red  peak  due  to  the  higher  or- 
der effects.   We  conclude  that,  for  the  density  range  of  this  work,  the  plasma 
microfield  is  not  strong  enough  for  these  higher  order  effects  to  become  impor- 
tant.  As  the  average  plasma  microfield  strength  continues  to  increase  due  to 
increasing  density,  it  will  eventually  become  of  the  same  order  as  the  central 
Coulomb  field.  At  that  point  the  higher  order  effects  will  be  much  more  likely 
to  be  manifest.    If  we  go  much  beyond  1  x  1025  cm-3,  however,  the  electron 
degeneracy  will  have  to  be  incorporated  into  the  theory  as  we  pointed  out  in 
section  2.1. 

4.2  The  Lyman  0  Line 

The  Lfj  line  has  its  excited  state  electron  in  the  n  =  3  state.  This  makes 
it  much  more  susceptible  to  the  effect  of  the  plasma  microfield  than  in  the 
La  case.  To  investigate  the  higher  order  effects  on  this  line,  we  will  include 
the  ion-quadrupole  effect,  field  dependent  atomic  physics,  the  presence  of  a 
slight  resonance  width,  T,  for  each  energy  level,  and  the  overlap  of  the  levels 
of  principal  quantum  number  n  =  3  and  4. 

This  last  process  allows  for  atomic  matrix  elements  connecting  the  two 
principal  quantum  number  manifolds.  It  also  contributes  more  terms  to  the 
sum  over  intermediate  states  that  is  performed  when  evaluating  the  electron 
broadening  operator  Af(Aw).  In  the  field-free  atomic  physics  calculation,  this 
is  also  possible  but  the  energy  levels  are  widely  separated  at  their  unperturbed 
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energy  values.    Consequently,  their  contribution  to  the  sum  over  intermedi- 
ate states  will  be  minimal.  The  decreasing  value  for  large  Aw  of  the  electron 
broadening  many-body  function  G(Aw)  will  attenuate  each  added  term  if  the 
energy  separation  between  levels  is  greater  than  the  plasma  frequency  ujPjC. 
In  the  field  dependent  atomic  physics  picture,  the  zero-order  energy  levels  are 
perturbed  by  the  Stark  effect  and  can  overlap.  When  this  occurs,  the  G(Au) 
function  takes  its  maximum  value  for  each  intermediate  state  and  the  total 
electron  broadening  term  will  be  larger.    In  general,  the  red  and  blue  levels 
of  each  principal  quantum  number  manifold  interact  most  strongly  with  each 
other,  while  the  red  levels  of  one  principal  quantum  number  and  the  blue  lev- 
els corresponding  to  another  interact  hardly  at  all.61  This  quasi-selection  rule 
causes  the  electron  broadening  of  the  n  =  3  level  to  be  most  influenced  by  the 
red  levels  of  the  n  =  4  level.  This  mechanism  contributes  a  further  blue  asym- 
metry to  the  Lp  line  by  way  of  a  broader  red  wing  caused  by  this  additional 
electron  broadening.   These  electron  broadening  effects  are  illustrated  in  Fig. 
(4.4)  where  we  show  the  Lp  line  at  ne  =  1  x  1025cm~3  and  kT  =  800  eV  for 
the  case  of  full  field  dependent  atomic  physics  plus  the  ion-quadrupole  effect. 
This  is  compared  to  the  same  line  but  with  all  connections  to  the  n  =  4  levels 
excluded  from  the  electron  broadening  of  the  n  =  3  to  1  transition.   This  is 
accomplished  by  limiting  the  sum  over  intermediate  states  in  the  evaluation 
of  M(Au)  to  dipole  matrix  elements  connecting  n  =  3  levels  only  with  other 
n  =  3  levels.  The  sums  for  the  n  =  4  to  1  transition,  however,  contain  all  n  =  3 
and  4  levels.   The  n  =  3  to  1  line  with  the  normal  sum  over  all  n  =  3  and  4 
intermediate  states  is  approximately  10%  broader  at  half  maximum  intensity 
than  the  line  with  the  restriction  on  the  electron  broadening  sum. 
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Fine  structure  splitting  for  this  line  is  on  the  order  of  2  eV  (see  Appendix 
E).  For  ne  =  1  x  1024cm-3,  the  overall  line  width  is  about  40  eV.  We  will 
ignore  the  much  smaller  fine  structure  splitting.  It  has  been  shown83'42  that 
this  splitting  introduces  a  slight  red  line  asymmetry  which  is  obscured  by  the 
blue  asymmetry  effects  as  the  density  increases.  The  overall  shift  of  the  line 
due  to  the  fine  structure  is  not  insignificant  but  is  ignored  here  because  we  are 
primarily  interested  in  the  line  shape.  At  these  densities  it  is  difficult  to  de- 
termine exact  line  positions  due  to  uncertainties  in  the  experimental  methods. 
We  will  also  examine  the  possible  influence  of  field  induced  ionization  on  the 
line  shape. 

Before  we  look  at  the  profiles  of  the  Lp  line  in  detail,  we  look  at  some  energy 
level  diagrams  for  the  n  =  3  and  4  levels  as  function  of  the  field  strength.  Fig. 
(4.5)  gives  the  energy  levels  as  a  function  of  the  field  for  the  case  of  field 
dependent  atomic  physics  but  no  ion-quadrupole  effect.  The  manifolds  for  the 
two  levels  begins  to  overlap  at  e  values  slightly  greater  than  one.  Recall  that 
I  is  our  unitless  scaled  field  magnitude.  According  to  the  Inglis-Teller  limit,84 
this  is  the  point  where  the  two  spectral  lines  merge  and  become  indistinct.  This 
will  only  occur,  however,  if  both  lines  are  of  equal  intensity.  In  a  spectral  series 
the  intensity  of  the  lines  generally  decreases  as  you  go  up  in  the  series  (as  the  n 
value  of  the  upper  state  increases).  As  this  occurs,  we  will  find  that  the  L7  line 
merges  with  the  Lp  line  but  the  Lp  line  itself  will  still  be  plainly  visible  due 
to  its  much  higher  intensity.  The  Inglis-Teller  limit  is  more  accurately  applied 
to  lines  corresponding  to  large  values  of  the  principal  quantum  number,  n. 
For  that  case  the  intensities  of  adjacent  lines  are  similar  and  their  merger  will 
produce  an  indistinguishable  whole. 
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In  Fig.  (4.5)  we  can  see  the  deviation  from  the  linear  Stark  effect  caused  by 
the  higher  order  corrections  to  the  energy.  This  is  most  prominent  in  the  n  =  4 
levels.  In  Fig.  (4.6)  the  Hamiltonian  is  diagonalized  with  the  ion-quadrupole 
terms  added.    Now  the  levels  from  the  two  manifolds  can  interact  strongly 
through  the  field  dependent  quadrupole  moment  matrix  elements  of  the  ion- 
quadrupole  term.    There  is  a  downward  shifting  of  red  energy  levels  of  both 
manifolds  and  a  slight  bunching  together  of  the  blue  levels.  When  the  levels 
approach  each  other  between  these  two  manifolds,  there  is  strong  interaction. 
Level  mixing  and  repulsion  is  plainly  evident.  The  highest  energy  (bluest)  n  = 
3  level  experiences  a  strong  avoided  crossing  with  the  lowest  energy  (reddest) 
n  =  4  level.  This  diagram  was  calculated  from  a  direct  diagonalization  of  the 
Hamiltonian  for  the  n  =  3  and  4  levels  using  field  dependent  atomic  physics. 
The  previous  figure  is  composed  of  the  radiator  energy  levels  calculated  directly 
from  Schrodinger's  equation  containing  the  ion-dipole  term.    In  a  plot  of  the 
line  shape,  the  actual  positions  of  these  transition  energies  become  somewhat 
blurred  due  to  the  added  imaginary  electron  broadening  and  resonance  width 
terms.    Figs.    (4.5)  and  (4.6)  do  give  us,  however,  a  good  indication  of  how 
the  levels  shift  and  interact  as  a  function  of  the  field.    The  presence  of  the 
ion-quadrupole  effect  not  only  provides  for  direct  energy  shifts  of  the  levels, 
but  it  also  allows  the  levels  from  different  principal  quantum  number  manifolds 
to  interact  resulting  in  complicated  level  mixing  and  repulsion. 

The  Lp  line  without  any  higher  order  field  effects  or  fine  structure  is  an 
almost  completely  symmetric,  double  peaked  line.  The  slight  asymmetry  is 
due  to  the  frequency  dependence  of  the  electron  broadening  operator  G(Au). 
The  two  levels  in  the  center  of  the  n  =  3  manifold  have  electric  quantum 
number  q  =  0  and  are  dipole  forbidden  (see  Fig.   (4.5)).  This  leaves  only  the 
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outer  components  which  split  linearly  with  the  field  thus  producing  the  two 
nearly  symmetric  peaks.  When  the  ion-quadrupole  effect  is  added,  the  high 
energy  levels  bunch  together  slightly  while  the  low  energy  levels  are  deflected 
downward  imparting  a  blue  asymmetry  to  the  line.  This  asymmetry  can  be 
seen  to  increase  as  the  density  increases  in  Figs.  (4.7),  (4.8)  and  (4.9).  In  these 
figures  the  effect  of  the  n  =  4  levels  is  also  apparent.  The  Ly  line  is  centered 
at  4133  eV.  In  Fig.  (4.7),  at  an  electron  density  of  ne  =  1  x  1024cm-3  we  see 
a  slight  increase  in  the  intensity  of  the  blue  wing,  but  as  the  density  increases 
a  broad,  flat,  blue  shoulder  develops.  At  5  x  1024cm-3,  the  Ly  line  is  flat 
having  merged  with  the  La  to  become  the  prominent  shoulder  visible  in  Fig. 
(4.8).  This  illustrates  the  point  we  made  earlier  about  the  Inglis-Teller  limit 
applying  only  to  the  less  intense  member  of  a  pair  of  merging  lines. 

Having  established  the  ion-quadrupole  effect  as  a  source  of  line  asymmetry, 
we  will  return  to  it  later  to  examine  the  difference  in  theoretical  approaches 
to  its  calculation.  Now,  we  examine  the  effect  of  the  field  dependent  atomic 
physics  on  the  L@  line  along  with  the  ion-quadrupole  effect.  From  the  energy 
level  diagram,  we  expect  the  added  field  dependence  to  produce  further  blue 
asymmetry  of  the  lineshape.  In  Fig.  (4.10)  we  see  this  beginning  for  ne  = 
1  x  1024cm-3.  This  trend  continues  in  Fig.  (4.11)  at  ne  =  5  x  1024cm-3. 
A  check  of  these  additional  asymmetries  against  an  approximate  quadratic 
Stark  effect  calculation83'42  shows  that  this  asymmetry  effect  is  due  primarily 
to  the  quadratic  correction  to  the  energy.  At  1  x  10  cm-3,  however,  the 
added  asymmetry  goes  beyond  the  quadratic  correction  and  is  due  to  higher 
order  energy  shifts  and  field  dependent  changes  in  the  atomic  matrix  elements. 
This  strong  additional  asymmetry  can  be  seen  in  Fig.  (4.12).  In  Figs.  (4.13), 
(4.14)  and  (4.15)  we  compare  the  full  field  dependent  atomic  physics  line  shape 
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calculations  for  a  range  of  densities,  with  the  field  independent  atomic  physics 
calculations  with  the  quadratic  Stark  effect  correction  to  the  level  energies. 
Both  cases  also  have  the  ion-quadrupole  effect.    The  quadratic  Stark  effect 
correction  consists  of  adding  the  second  order  perturbation  theory  term  from 
the  treatment  of  an  ion  in  a  uniform  electric  field,  to  the  upper  state  radiator 
energy  levels.    This  amounts  to  throwing  out  all  field  dependence  from  the 
atomic  physics  except  this  diagonal  energy  correction  second  order  in  the  field. 
It  can  be  assumed  that  remaining  differences  between  the  two  line  shapes  in 
the  illustrations  is  due  to  additional  field  dependence  of  the  energies,  matrix 
elements  and  intensity  factors.  We  conclude  that  at  ne  =  5  x  1024  cm-3  and 
below,  the  field  independent  atomic  physics  plus  the  approximate  quadratic 
Stark  effect  correction  is  a  good  approximation  to  the  full  field  dependent 
results  for  the  Lp  line  of  hydrogenic  argon. 

In  section  3.4  we  discussed  the  effect  of  the  possible  field  ionization  pro- 
duced when  the  plasma  microfield  becomes  comparable  to  the  strength  of  the 
Coulomb  field  experienced  by  the  radiator  electron.  We  pointed  out  that  the 
real  situation  was  probably  somewhere  between  the  extremes  of  field  induced 
ionization  in  a  uniform  field  and  no  ionization  at  all  due  to  inverse  processes 
including  exchange  of  bound  electrons  between  neighboring  ions  with  quasi- 
molecular  orbitals.   We  will  examine  the  limiting  case  of  field  ionization  in  a 
uniform  field  to  estimate  the  maximum  possible  reduction  in  line  intensity  due 
to  population  depletion.   In  Fig.   (4.16)  we  see  a  reduction  in  intensity  in  the 
line  wings  that  gives  rise  to  a  slightly  sharper  line  with  no  change  in  the  line 
asymmetry.    As  density  increases  we  see  in  Fig.    (4.17)  a  further  narrowing 
of  the  line  and  an  additional  reduction  in  the  intensity  of  the  red  wing.   This 
follows  from  the  fact  that  the  lower  energy  states  are  more  strongly  affected 
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by  the  field  and  hence,  ionize  first.  Again  the  peak  asymmetry  is  not  affected. 
This  trend  continues  in  Fig.  (4.18)  at  ne  =  1  x  1025cm-3.  This  possible 
line  narrowing,  if  real,  could  have  consequences  for  density  diagnostics.  The 
narrower  lines  could  be  interpreted  as  indicative  of  a  somewhat  lower  density. 
However,  the  strong  asymmetries  would  still  be  present  indicating  the  true 
higher  density  of  the  plasma. 

The  resonance  widths  T  that  are  included  in  these  calculations  do  not  add 
significant  width  or  distortion  to  the  line  shape.  This  is  due  to  the  resonance 
width  becoming  appreciable  only  where  the  microfield  probability  function 
becomes  negligibly  small.  For  example,  for  the  most  red  n  =  3  level  (the  level 
most  strongly  affected  by  the  field),  the  field  dependent  resonance  width  is 
T  ~  10~5  eV  for  the  field  value  e  w  1  and  T  ~  1  eV  for  e  w  2  compared  to  a 
Lp  line  width  of  approximately  150  eV.  For  our  argon  plasma  at  1  x  1025  cm-3 
and  800  eV,  when  I  =  2  the  microfield  function  P(e)  has  fallen  to  only  25% 
of  its  peak  value  which  occurred  at  e  «  0.8.  The  resonance  width,  and  hence 
resonance  lifetime,  has  a  larger  possible  effect  on  the  level  populations  as  we 
have  seen. 

Now  we  return  to  the  consideration  of  different  treatments  of  the  ion- 
quadrupole  effect.  Previous  calculations12  of  the  ion-quadrupole  effect  con- 
sidered only  ion-radiator  correlations  while  neglecting  ion-ion  correlations.  We 
examined  the  importance  of  these  correlations  by  comparing  the  results  of 
our  calculations  that  included  such  ion-ion  correlations  with  a  simple  nearest 
neighbor  approximation  which  did  not.  In  Fig.  (4.19)  we  see  the  results  for 
ne  =  1  x  10  4cm~3.  As  expected,  the  calculation  with  no  correlations  shows 
the  strongest  asymmetry  as  explained  in  section  3.1.  This  trend  continues  in 
Fig.   (4.20)  for  ne  =  5  x  1024cm-3.  At  ne  =  1  x  1025cm-3,  however,  we  see 
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in  Fig.    (4.21)  that  the  two  calculations  of  the  ion-quadrupole  term  begin  to 
approach  each  other.    This  is  a  reflection  of  the  high  field  limit  of  the  field 
gradient  average  as  illustrated  in  Figs.  (3.1)  and  (3.2). 

This  concludes  our  examination  of  the  calculational  results  of  our  inves- 
tigation of  certain  higher  order  field  effects  on  spectral  line  shapes  in  dense 
plasmas. 

4.3  Conclusions 

We  conclude  this  investigation  by  summarizing  our  findings,  pointing  out 
some  of  their  consequences  and  suggesting  the  next  steps  needed  to  extend  this 
work  to  higher  densities. 

For  our  range  of  physical  conditions  (ne  =  1  x  1024  cm-3  to  1  x  1025  cm-3 
at  kT  =  800  eV)  and  for  the  pure  argon  plasma  we  found  that  the  higher  order 
field  effects  begin  to  become  important  for  the  n  =  3  to  1  transition  (Lp).  The 
n  =  2  to  1  transition  (LQ)  showed  only  very  slight  perturbations  from  the  linear 
Stark  effect  with  fine  structure  splitting  and  Doppler  broadening.  While  the 
fine  structure  does  show  a  slight  field  dependence,  it  produces  no  significant 
effect  on  the  line  splitting  or  shape  for  the  LQ  line.  The  Lp  line  did  show  many 
of  the  additional  field  effects.  The  most  pronounced  was  the  effect  of  the  field 
dependent  atomic  physics  (matrix  elements  and  energy  levels  that  include  an 
exact  treatment  of  the  ion-dipole  term  in  the  radiator  Hamiltonian).  This  field 
dependence  can  be  approximated  by  simply  adding  the  quadratic  correction  to 
the  level  energy  for  densities  of  1  x  1024  to  5  x  1024  cm-3.  At  ne  =  1  x  1025  cm-3, 
however,  the  other  field  effects  on  the  atomic  physics  becomes  important  and 
can  not  be  ignored. 
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We  have  also  seen  that  a  more  accurate  treatment  of  the  ion-ion  correlations 
in  the  ion-quadrupole  effect  reduces  the  line  asymmetry  usually  produced.  At 
the  very  highest  densities  that  we  examined,  a  simple  nearest  neighbor  model 
produces  a  fair  approximation  to  the  more  accurate  results.  For  densities  higher 
than  these  the  nearest  neighbor  approximation  may  be  quite  sufficient. 

Our  examination  of  possible  field  ionization  effects  has  shown  that  this 
could  have  effects  on  the  spectral  line  shape,  producing  a  minimal  narrowing 
of  the  line  at  high  density.  This  point  is  highly  speculative,  however,  and  no 
firm  conclusions  can  be  drawn  until  a  much  more  detailed  study  is  carried  out. 
This  would  probably  entail  a  many-body  theory  treating  plasma  electrons  and 
quasibound  electrons  on  a  more  equal  footing. 

We  have  also  seen  that  it  is  possible  to  carry  the  line  shape  calculation 
to  densities  beyond  the  Inglis-Teller  merging  point.  We  included  the  n  =  4 
levels  in  our  calculation  of  the  n  =  3  to  1  transition  and  found  that  transitions 
from  the  n  =  4  level  formed  a  prominent  shoulder  on  the  high  energy  side  of 
the  Lp  line.  We  concluded,  at  least  for  this  case,  that  the  Inglis-Teller  limit 
applies  to  the  significantly  lower  intensity  member  of  a  neighboring  pair  of  line 
transitions  and  does  not  indicate  the  complete  merging  of  the  two  lines.  The 
Inglis-Teller  limit  is  more  appropriate  for  higher  series  members  where  the  two 
line  intensities  are  nearly  equal. 

To  carry  these  line  shape  calculations  to  higher  densities  for  the  same  tem- 
perature range,  it  will  be  necessary  to  include  a  treatment  of  degenerate  plasma 
electrons.  We  showed  in  chapter  2  that  for  densities  of  ne  =  1  x  1026  cm-3  the 
Fermi  energy  is  of  the  same  order  as  the  thermal  energy.  When  the  plasma 
electrons  become  degenerate  they  may  not  be  available  for  electron  collisional 
broadening.85  This  may  reduce  the  electron  broadening  at  high  densities.  To 
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treat  higher  series  members  at  these  densities,  it  will  be  necessary  to  develop 
a  more  accurate  way  to  deal  with  the  width  of  the  resonance  state.  This  may 
be  possible  through  the  use  of  a  line  shape  formalism  that  treats  bound  and 
continuum  states  on  an  equal  footing. 

At  the  present  time  the  accurate  experimental  verification  of  these  higher 
order  field  effects  is  problematical,  but  this  may  change  as  X-ray  spectroscopic 
techniques  improve.  In  conclusion,  there  appears  to  be  no  fundamental  reason 
why  the  spectroscopy  of  high  density  plasmas  can  not  be  pushed  to  even  more 
extreme  conditions  than  are  successfully  dealt  with  today.  In  these  regimes 
of  ultra-high  densities,  we  can  expect  to  continue  to  find  new  and  interesting 
physical  phenomena  in  need  of  explanation. 
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FIGURE  (4.1)  The  La  line  of  hydrogenic  argon  for  ne  =  1  x  1024  cm-3  and 
kT  =  800  eV  as  discussed  in  section  4.1.   The  profile  with  field  independent 

atomic  physics  and  no  ion  quadrupole  effect  is  given  by  ( ).  The  profile  with 

field  dependent  atomic  physics  and  the  APEX  ion- quadrupole  effect  is  given 

by  ( ).  Both  profiles  have  the  fine  structure  and  are  Doppler  convoluted. 

They  are  calculated  from  Eqs.  (2.97)  and  (3.14). 
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FIGURE  (4.2)  The  La  line  of  hydrogenic  argon  for  ne  =  5x  1024 cm-3  and 
kT  =  800  eV.    The  profile  with  field  independent  atomic  physics  and  no  ion 

quadrupole  effect  is  given  by  ( ).  The  profile  with  field  dependent  atomic 

physics  and  the  APEX  ion-quadrupole  effect  is  given  by  ( ).  Both  profiles 

have  the  fine  structure  and  are  Doppler  convoluted. 
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FIGURE  (4.3)  The  LQ  line  of  hydrogenic  argon  for  ne  =  1  x  1025  cm-3  and 
kT  =  800  eV.   The  profile  with  field  independent  atomic  physics  and  no  ion 

quadrupole  effect  is  given  by  ( ).  The  profile  with  field  dependent  atomic 

physics  and  the  APEX  ion-quadrupole  effect  is  given  by  ( ).  Both  profiles 

have  the  fine  structure  and  are  Doppler  convoluted. 
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FIGURE  (4.4)  Comparison  of  the  Lp  line  of  hydrogenic  argon  with  and 
without  electron  broadening  connections  to  the  n  =  4  levels.  The  normal 
line  shape  with  field  dependent  atomic  physics  and  the  ion-quadrupole  effect 

is  given  by  ( ).    The  same  line  shape  but  with  all  dipole  matrix  elements 

connecting  the  n  =  3  and  4  levels  deleted  from  the  sum  in  Eq.  (2.96)  which  is 

used  to  evaluate  the  electron  broadening  operator  M(Auj)  is  given  by  ( ). 

Here  ne  =  1  x  1025  cm"3  and  kT  =  800  eV. 
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FIGURE  (4.5)  Energy  level  diagram  for  the  n  =  3  and  4  manifolds  of  levels 
as  a  function  of  the  plasma  microfield  using  field  dependent  atomic  physics  but 
no  ion-quadrupole  effect.  Figure  discussed  in  section  4.2.  The  n  =  3  levels  are 

given  by  ( )  and  the  n  =  4  levels  are  given  by  ( ).  The  field  is  given 

by  the  unitless  e  =  e/e0  and  the  energy  is  in  Rydbergs. 
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FIGURE  (4.6)  Energy  level  diagram  for  the  n  =  3  and  4  manifolds  of  levels 
as  a  function  of  the  plasma  microfield  using  field  dependent  atomic  physics  and 
the  ion-quadrupole  effect.  Figure  discussed  in  section  4.2.  The  n  =  3  levels  are 

given  by  ( )  and  the  n  =  4  levels  are  given  by  ( ).  The  field  is  given 

by  the  unitless  e  =  e/e0  and  the  energy  is  in  Rydbergs.  The  ion-quadrupole 
corrections  are  calculated  in  the  APEX  approximation  for  ne  =  1  x  10  5  cm-'' 
and  kT  =  800  eV. 
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FIGURE  (4.7)  The  Lp  line  of  hydrogenic  argon  for  ne  =  1  x  1024  cm-3  and 
kT  =  800  eV.   The  profile  with  field  independent  atomic  physics  and  no  ion 

quadrupole  effect  is  given  by  ( ).  The  profile  with  field  independent  atomic 

physics  and  the  APEX  ion-quadrupole  effect  is  given  by  ( ).  Both  profiles 

are  Doppler  convoluted.  They  are  calculated  from  Eqs.  (2.97)  and  (3.14)  and 
discussed  in  section  4.2. 
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FIGURE  (4.8)  The  Lp  line  of  hydrogenic  argon  for  ne  =  5  x  1024  cm-3 
and  kT  =  800  eV.   The  profile  with  field  independent  atomic  physics  and  no 

ion  quadrupole  effect  is  given  by  ( ).    The  profile  with  field  independent 

atomic  physics  and  the  APEX  ion-quadrupole  effect  is  given  by  ( ).  Both 

profiles  are  Doppler  convoluted. 
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FIGURE  (4.9)  The  Lp  line  of  hydrogenic  argon  for  ne  =  1  x  1025  cm-3 
and  kT  =  800  eV.   The  profile  with  field  independent  atomic  physics  and  no 

ion  quadrupole  effect  is  given  by  ( ).    The  profile  with  field  independent 

atomic  physics  and  the  APEX  ion- quadrupole  effect  is  given  by  ( ).  Both 

profiles  are  Doppler  convoluted. 
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FIGURE  (4.10)  The  Lp  line  of  hydrogenic  argon  for  ne  =  1  x  10 
and  kT  —  800  eV.     The  profile  with  field  independent  atomic  physics  and 

the  APEX  ion  quadrupole  effect  is  given  by  ( ).     The  profile  with  field 

dependent  atomic  physics  and  the  APEX  ion-quadrupole  effect  is  given  by 
( ).  Both  profiles  are  Doppler  convoluted. 
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FIGURE  (4.11)  The  Lp  line  of  hydrogenic  argon  for  ne  =  5  x  1024cm-3 
and  kT  =  800  eV.     The  profile  with  field  independent  atomic  physics  and 

the  APEX  ion  quadrupole  effect  is  given  by  ( ).    The  profile  with  field 

dependent  atomic  physics  and  the  APEX  ion- quadrupole  effect  is  given  by 
( )•  Both  profiles  are  Doppler  convoluted. 
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FIGURE  (4.12)  The  Lp  line  of  hydrogenic  argon  for  ne  =  1  x  1025cm-3 
and  kT  =  800  eV.     The  profile  with  field  independent  atomic  physics  and 

the  APEX  ion  quadrupole  effect  is  given  by  ( ).    The  profile  with  field 

dependent  atomic  physics  and  the  APEX  ion- quadrupole  effect  is  given  by 
( ).  Both  profiles  are  Doppler  convoluted. 
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FIGURE  (4.13)  The  Lp  line  of  hydrogenic  argon  for  ne  =  1  x  10 
and  kT  =  800  eV.  The  profile  with  field  dependent  atomic  physics  and  the  ion 

quadrupole  effect  is  given  by  ( ).  The  profile  with  field  independent  atomic 

physics,  the  ion- quadrupole  effect  and  the  quadratic  energy  correction  is  given 
by  ( ).  Both  profiles  are  Doppler  convoluted. 
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FIGURE  (4.14)  The  Lp  line  of  hydrogenic  argon  for  nt  =  5  x  1024cm~3 
and  kT  =  800  eV.  The  profile  with  field  dependent  atomic  physics  and  the  ion 

quadrupole  effect  is  given  by  ( ).  The  profile  with  field  independent  atomic 

physics,  the  ion-quadrupole  effect  and  the  quadratic  energy  correction  is  given 
by  ( ).  Both  profiles  are  Doppler  convoluted. 
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FIGURE  (4.15)  The  Lp  line  of  hydrogenic  argon  for  ne  =  1  x  1025cm-3 
and  kT  =  800  eV.  The  profile  with  field  dependent  atomic  physics  and  the  ion 

quadrupole  effect  is  given  by  ( ).  The  profile  with  field  independent  atomic 

physics,  the  ion- quadrupole  effect  and  the  quadratic  energy  correction  is  given 
by  ( ).  Both  profiles  are  Doppler  convoluted. 
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FIGURE  (4.16)  The  Lp  line  of  hydrogenic  argon  for  ne  =  1  x  1024cm-3 
and  kT  =  800  eV.  The  profile  with  no  field  ionization,  field  dependent  atomic 

physics  and  the  ion  quadrupole  effect  is  given  by  ( ).    The  profile  with 

field  ionization,  field  dependent  atomic  physics  and  the  ion-quadrupole  effect 
is  given  by  ( ).  Both  profiles  are  Doppler  convoluted. 
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FIGURE  (4.17)  The  Lp  line  of  hydrogenic  argon  for  ne  =  5  x  1024cm-3 
and  kT  =  800  eV.  The  profile  with  no  field  ionization,  field  dependent  atomic 

physics  and  the  ion-quadrupole  effect  is  given  by  ( ).     The  profile  with 

field  ionization,  field  dependent  atomic  physics  and  the  ion-quadrupole  effect 
is  given  by  ( ).  Both  profiles  are  Doppler  convoluted. 
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FIGURE  (4.18)  The  Lp  line  of  hydrogenic  argon  for  ne  =  lx  1025cm-3 
and  kT  =  800  eV.  The  profile  with  no  field  ionization,  field  dependent  atomic 

physics  and  the  ion  quadrupole  effect  is  given  by  ( ).    The  profile  with 

field  ionization,  field  dependent  atomic  physics  and  the  ion- quadrupole  effect 

is  given  by  ( ).  Both  profiles  are  Doppler  convoluted  and  are  discussed  in 

section  4.2. 
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FIGURE  (4.19)  The  Lp  line  of  hydrogenic  argon  for  ne  =  1  x  1024cm~3 
and  kT  =  800  eV.  The  profile  with  field  dependent  atomic  physics  and  the 
APEX  ion  quadrupole  effect  is  given  by  ( ).  The  profile  with  field  depen- 
dent atomic  physics  and  the  nearest  neighbor  ion-quadrupole  effect  is  given  by 
( )•  Both  profiles  are  Doppler  convoluted. 
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FIGURE  (4.20)  The  Lp  line  of  hydrogenic  argon  for  ne  =  5  x  1024  cm-3 
^  ^T  =  800  eV.  The  profile  with  field  dependent  atomic  physics  and  the 
APEX  ion  quadrupole  effect  is  given  by  ( ).  The  profile  with  field  depen- 
dent atomic  physics  and  the  nearest  neighbor  ion-quadrupole  effect  is  given  by 
( )•  Both  profiles  are  Doppler  convoluted. 
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FIGURE  (4.21)  The  Lp  line  of  hydrogenic  argon  for  ne  =  1  X  1025cm~3 
and  kT  =  800  eV.  The  profile  with  field  dependent  atomic  physics  and  the 
APEX  ion  quadrupole  effect  is  given  by  ( ).  The  profile  with  field  depen- 
dent atomic  physics  and  the  nearest  neighbor  ion- quadrupole  effect  is  given  by 
( ).  Both  profiles  are  Doppler  convoluted. 


APPENDIX  A 
SOME  EXPRESSIONS  IN  USEFUL  UNITS 

We  will  present  here  some  frequently  encountered  physical  quantities  ex- 
pressed in  various  useful  units.  The  reason  for  this  is  purely  one  of  convenience. 
I  have  found  that  I  constantly  recalculate  these  things  because  I  did  not  think 
it  important  enough  to  write  them  down  where  I  could  find  them  again.  We 
will  also  discuss  various  units  used  in  atomic  physics  and  line  broadening  calcu- 
lations because  a  long  history  and  a  mixing  of  several  disciplines  has  produced 
a  state  of  utter  confusion  in  some  places. 

We  give  several  useful  quantities  from  chapter  2  in  simple  numerical  form. 
In  these  expressions,  rij  is  given  in  cm-3,  kT  is  given  in  eV  and  Zj  is  the  unit 
of  charge  for  species  j  (i.e.  for  argon  Z  =  18  and  for  electrons  Z  =  1).  We 
give  lengths  in  cm  but  sometimes  a  more  convenient  unit  for  atomic  scales  is 
the  Angstrom  with  1  A  =  10~8  cm.  Also,  the  correspondence  between  energy 
and  temperature  is  1  eV  «-»  1.1604  x  104  deg  K  or  approximately  11,000  deg 
K.  Other  relations  are  given  in  the  following  list. 

Ion  sphere  radius: 


roJ  =  0.62035  n-  1/3  cm 


(A.l) 


Debye  sphere  radius: 


\DJ  =  1.6293  x  10-2Zr\kT)l/2nj1/2cm. 


(A.2) 


Electron  thermal  wavelength: 


123 


124 


Ae  =  6.9195  x  lO-^JbT)"1/2  cm  .  (A.3) 

Atomic  Bohr  radius: 

rn  =  5.2918  x  10-9  Z~ln2  cm  .  (A.4) 

Electron  plasma  frequency: 

hup,e  =  3.7133  x  10-11  nj2  eV  .  (A.5) 

Electron-electron  coupling  constant: 

re,e  =  2.3212  x  10~7 \kT)-ln\,Z  .  (A.6) 

Plasma  mass  density: 

p  =  1.6605  x  10-24^negcm-3  ,  (A.7) 

where  Ay/  is  the  ion  atomic  mass  in  atomic  mass  units  and  Z  is  the  average 
effective  ionic  charge. 

For  formulae  expressed  in  atomic  units,  we  take  h  =  e  =  me  =  1  and  energy 
is  in  units  of  Hartrees.  We  have  1  Rydberg  =  l/2Hartree  =  13.606  eV.  Com- 
mon conversion  factors  are  1  eV  =  1.6022  x  10_12erg  =  8.0655  x  103  cm-1.  To 
convert  the  electric  field  strength  F,  in  atomic  units  with  Z  =  1,  to  the  unit- 
less  electric  field  I  (e  is  commonly  used  to  give  the  plasma  microfield  function 
P(e)),  we  have 

e  =  1.37425  x  1016Z3n~2/3  F  .  (A.8) 

Many  other  useful  factors  and  relations  can  be  found  in  Ref.    86. 


APPENDIX  B 
THE  APEX  CONDITIONAL  DISTRIBUTION  FUNCTION 

In  this  appendix  we  will  derive  an  expression  for  the  transform  of  the 
conditional  distribution  function57  in  the  APEX  approximation.  We  will  also 
discuss  the  APEX  approximation  to  the  low  frequency  microfield  distribution 
function  on  which  the  conditional  distribution  depends. 

The  conditional  distribution  function  can  be  expressed  as 


g(n  A)  = 


1   6  (in  0(A)) 


n ; 


(B.l) 


E*(r)  =  Z,e— (1  +  ar)e 
and  the  plasma-ion  Debye-Huckel  field  is 


—ar 


(B.3) 


1 

E{f)  =  Z,e— (1  +  Kr)e 


— «r 


(B.4) 


H  8  (i\  ■  E(f))  ' 

where  n,-  is  the  ion  number  density  and  Q{\)  is  the  Fourier-transformed  mi- 
crofield function.  We  can  use  a  model  for  Q{\)  to  obtain  a  corresponding 
model  for  g(f ;  A)  which  can  be  transformed  to  give  the  conditional  distribu- 
tion function  g(f ;  e). 

With  the  APEX  model,54 

InQ(A)  =  ni  J  drg(r)jj&  (•**«  -  l)  ,  (B.2) 

where  the  APEX  renormalized  field,  E*(r),  is 
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12C 


Here,  a  is  the  effective  inverse  screening  length5    and  k  =  l/Xj)  e  is  the  inverse 
Debye-Hiickel  screening  length.  These  expressions  are  obviously  related  by 


^-mm' 


(B.5) 


where 


R(r) 


(1  +  /cr)e 


— nr 


(B.6) 


(1  +  ar)e-ar 

Thus  E*(f)  is  a  functional  of  E(f):  i.e.   E*  =  E*[E\.   Now  we  can  construct 
the  APEX  version  of  g(f;  A)  by  using  Eq.  (B.2)  in  Eq.  (B.l).  We  obtain 


W\  A) 


J  i\-E*(r')\  > 


8(i\-E(r)}J  'i\-E*(r') 

J  ar  ^r>E*{r')  8E(r) 

+   [dr>       9{r<)        (&'*V)  -  l)  ^E{r>)        E{r>)  ^^ 
J         i\-E*(r')^  I 


SE(r)       E*(r')  8E(r) 


The  second  term  gives  zero  and  leaves  us  with  the  result, 


(B.7) 


E(r>) 


g(r;  A)  =   /  dr '  g(r' 
Carrying  out  the  integration,  we  obtain 


}E*{r>)  {  'Rir') 


(B.8) 


9APEx(n  A)  =  g(r)eiXS'^  , 


(B.9) 


which  is  the  expression  we  desired.  In  the  above  development,  we  have  made 
use  of  the  expression 


8E(f') 
8E(f) 


=  8(r-r'), 


(B.10) 
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which  gives  us 


8E\r>)  _8E\r')8E{r>)  1 

8E(r)      '   8E(r>)  8E{r)         [V      V  } R(r') 


(B.ll) 


Eq.  (B.9)  will  allow  us  to  calculate  the  constrained  average  (F)€  in  the  simple 
APEX  approximation. 

This  approximation  is  based  on  the  APEX  model  for  the  low  frequency 
microfield  distribution  function.  We  now  briefly  review  the  key  features  of 
this  approximation.  A  simple  physical  picture  of  the  APEX  model  consists  of 
replacement  of  the  usual  screened-ion  plasma  by  a  plasma  of  noninteracting 
quasiparticles.  Each  of  which  will  produce  a  field  at  the  radiator.  However, 
this  field  has  an  adjustable  parameter  a  which  is  determined  by  the  imposi- 
tion of  two  constraints:53-54  the  second  moment  condition  and  the  local-field 
constraint.  The  single  particle  APEX  electric  field  is  given  by  Eq.  (B.3). 

The  local-field  constraint  requires  that  the  field  produced  at  the  radiator 
due  to  the  sum  of  the  adjustable  parameter  fields  from  the  quasiparticles  in  a 
volume  element  d3r  at  r  be  equal  to  the  field  produced  by  the  actual  plasma 
ions  interacting  with  the  chosen  inter-ion  potential.  In  our  case  we  choose  the 
Debye-Hiickel  screened  potential  to  describe  the  actual  ion  interactions  relevant 
to  the  low  frequency  plasma  microfield.54  If  we  define  the  radial  distribution 
function  for  the  quasiparticles  surrounding  the  radiator  as  G(r),  the  local  field 
constraint  requires  that 


niG(r)E*(r)d3,*  =  nig(r)E(r)d3r  , 


or 


(B.12) 
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G(r)  =  g(r) 


E(r) 
E*(r) 


(B.13) 


The  second  moment  condition  requires  that  the  APEX  approximation  to 
Q(e)  produce  the  exact  second  moment.  Specifically,  that 


/' 


i*J2 


dee'QAPEX{e)  =  (&)  , 


(B.14) 


where  E  is  the  total  field  produced  by  the  ions.  Again,  the  ions  are  assumed 
to  interact  with  Debye-Huckel  potentials. 

The  right-hand  side  of  Eq.  (B.14)  can  be  rewritten  as54 


m  = 


(Xe)2 

1 


(VV  •  W) 
(V2V) 


A:7T 


where 


(B.15) 


0(A)  =  -^  j^ '  dvrg{r)e-rlX. 


(B.16) 


Here  fi  =  1/kT  and  x  1S  the  effective  radiator  charge. 

The  left-hand  side  is  evaluated  using  the  local  field  constraint  to  give 

J dee2QAPEX(e)  =  J dee2 (6  (e-  E*))G 

=  (e*2)g 

r°°      «  o  (B-17) 

=  4tt  /      drr2G{r)E*\r) 

JO 
/•oo 

=  4tt  /      drr2g(r)E(r)E*(r)  . 

JO 

Here  the  subscript  G  on  (•  •  -)q  means  that  we  use  G(r)  rather  than  g(r)  when 
evaluating  the  ensemble  average  of  the  single  particle  additive  function  E*. 
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The  adjustable  parameter  a  is  then  determined  by  Eqs.    (B.15)  and  (B.17) 
which  gives 


ZiUi 


rOO  too 

/      drrg(r)e-r/X  =/3\2         drr2g(r)E(r)E*(a,r)  .  (B.18) 

JO  Jo 


We  evaluate  this  equation  using  an  HNC  calculation39  for  g(r)  and  the 
value  of  a  calculated  above.  We  are  then  free  to  evaluate  the  APEX  microfield 
distribution  function  given  by 


Qapex(^) 


f     d) 

J     (27T 


d\_ 


-Xe 


Q(A), 


(B.19) 


where  the  APEX  approximation  for  Q(\)  is  given  by  Eq.  (B.2).  We  will 
use  this  approximation  for  the  microfield  probability  function  to  perform  the 
integration  over  field  values  in  the  expression  for  I(uj)  given  by  Eq.  (3.11). 


APPENDIX  C 
PARABOLIC  COORDINATES 

Parabolic  coordinates61  are  denned  in  terms  of  the  cartesian  coordinates 


by 


x  =  y/^T]  cos  (p 

y  =  \fin  sin  y 


z  =  ^a-r)), 


and  in  terms  of  spherical  coordinates  by 


(C.l) 


£  =  r(l  +cos0) 

■q  =  r(l  —  cos0) 

<p  =  tp  . 
The  parabolic  coordinates  are  defined  on  the  domains 

0<£  <oo 

0  <  rj  <  oo 

0  <  tp  <  2tt  . 
The  volume  element  in  parabolic  coordinates  is  given  by 


(C.2) 


(C.3) 


(fix  =  -(£  +  7i)dtdT}dv  , 


(C.4) 


and  the  Laplacian  is  given  by 


Ud 


+ 


Vir)  + 


{t  +  n)dtvdtj    U+ij)dT}\  'dnJ    tod? 


1  d2 


(C.5) 
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Once  the  Schrodinger  equation  in  parabolic  coordinates,  Eq.  (3.57),  is 
separated  and  solved,  the  solutions  for  the  coordinates  £  and  77  have  some 
simple  characteristics.71  Motion  along  the  downfield  coordinate  £  is  bound 
since  it  corresponds  to  the  part  of  the  potential  that  increases  without  bound. 
By  downfield  we  mean  in  the  direction  of  the  electric  field  vector.  Likewise, 
upfield  means  the  opposite  direction.  The  separation  constant  A  can  take  only 
discrete  values  related  to  ni  which  gives  the  number  of  wavefunction  nodes 
in  the  (  direction.  Motion  along  the  upfield  coordinate  77  is  unbound  since  it 
corresponds  to  the  part  of  the  potential  that  forms  a  potential  barrier.    The 

number  of  wavefunction  nodes  in  the  77  direction  is  given  by  n2  for  the  range 

between  the  origin  and  the  classical  turning  point. 

We  have  been  using  the  three  parabolic  quantum  numbers  nj,  n2  and  m/ 

to  describe  the  field  dependent  states.   Another  set  is  also  in  common  usage; 

these  are  the  electric  quantum  numbers62  n,  q  and  m/  where  q  =  n2  -  n:  and 

n  is  the  usual  principal  quantum  number  given  by  n  =  n2  +  n2  +  |m/|  +  1.  The 

electron  spin  can  be  included  in  this  description  by  the  addition  of  the  spin 

quantum  number  ms  which  takes  the  values  ±1/2. 


APPENDIX  D 
IMPORTANCE  OF  THE  ION-QUADRUPOLE  EFFECT 

In  this  appendix,  we  present  an  estimate  of  the  conditions  under  which  the 
ion-quadrupole  effect  should  become  important  relative  to  that  of  the  ion-dipole 
for  the  case  of  multielectron  radiators.15  The  ion-quadrupole  effect  attempts 
to  account  for  the  fact  that  perturbing  plasma  ions  do  not  produce  a  uniform 
electric  field  at  the  radiator.  As  long  as  the  perturbers  are  sufficiently  far  away, 
the  uniform  field  approximation  is  not  unreasonable.  As  the  density  increases, 
however,  these  perturbers,  on  the  average,  draw  closer  to  the  radiator  and  can 
produce  a  significant  field  gradient. 

To  help  us  understand  when  the  ion-quadrupole  effect  will  become  impor- 
tant, we  will  examine  in  some  detail  a  simple  model  system  that  consists  of 
two  close  lying  energy  levels  with  the  same  principal  quantum  number.  We 
will  look  at  how  the  dipole  and  quadrupole  effects  cause  the  states  associated 
with  these  two  levels  to  mix. 

Consider  the  radiator  Hamiltonian  with  terms  up  to  the  quadrupole  in  the 
multipole  expansion  of  the  radiator  perturbing-ion  interaction: 

H  =  H°  +  D  +  Q,  (D.l) 

where  H°  is  the  unperturbed  radiator  Hamiltonian  and  D  and  Q  are  the  dipole 
and  quadrupole  corrections  respectively.  We  assume  that  D  ^>  Q  for  the 
expansion  to  be  valid.  Both  D  and  Q  will  depend  on  the  plasma  electric 
microfield.    The  parity  of  D  is  odd  while  that  of  Q  is  even.    Because  of  this 


132 


133 

fact  D  can  have  non-zero  matrix  elements  only  between  states  of  unlike  parity 
while  Q  can  have  non-zero  matrix  elements  only  between  states  of  like  parity. 
This  means  all  diagonal  matrix  elements  of  D  vanish  for  basis  states  of  definite 
parity  such  as  those  produced  by  our  zero-order  Hamiltonian  if  .  For  this 
case,  we  will  take  the  states  a  and  b  to  be  of  different  parity.  This  will  give 
non-zero  off  diagonal  matrix  elements  for  the  dipole  operator  D  but  zero  off 
diagonal  matrix  elements  for  the  quadrupole  operator  Q.  We  will  represent  the 
Hamiltonian,  H,  for  this  two  state  system  in  terms  of  the  eigenstates  of  i?°, 


Hl 


a°)=E°a 


°0) 


Thus  Dab  =  (a0|D|&°)  and  Qa,a  =  (a°\Q\a°)  while  Da,a 
Qab  =  Qba  —^  because  of  the  parity  constraint. 

In  matrix  form,  the  Hamiltonian  is  then  represented  by, 


(D.2) 
Dbb  =  0  and 


H 


(D.3) 


E°a  +  Qa,a         Datb 
Db,a         EQb  +  QbM 

Construction  of  the  secular  equation  for  this  Hamiltonian  from  the  determi- 
nantal  equation  \H  -  EI\  =  0,  where  i"  is  the  identity  matrix,  will  yield  the 
system  eigenvalues: 


E  =  \  ($  +  E°b  +  Qa,a  +  Qbtb  ±  ^(A£0  +  AQa;6)2  +  4|i)a)6|2)  ,      (D.4) 

where  AEq  =  E®  —  E?  is  the  zero-order  energy  splitting  between  levels  with 
the  same  principal  quantum  number  and  AQab  =  Qa,a  —  Qb,b-  The  energy 
difference,  or  splitting,  between  the  two  levels  is  then 


AE 


y/(AEo  +  AQa.b)*+4\DaJb\*, 


(D.5) 
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where  AE  =  E+  —  E-.  This  can  be  expanded  to  give 


AE  =  yjAEl  +  AQlb  +  2AE0AQa.b  +  4|£>a;6|2  .  (D.6) 

We  take  the  eigenfunctions  of  H®  to  be  real  and  with  the  assumption  that 


\Da,b\  >  \&Qa;b\  , 


(D.7) 


we  can  ascertain  the  significance  of  the  quadrupole  term  on  AE.  The  validity  of 
this  assumption  can  be  checked  for  the  conditions  of  interest  with  the  estimates 
for  Da  i  and  AQa.b  given  below.  For  our  physical  conditions  of  interest,  it  holds 
true  but  it  should  be  checked  for  additional  cases. 

If  the  two  eigenvalues  of  if0  are  degenerate  we  will  have  AEq  =  0,  and 
the  dipole  term  will  dominate.  This  is  the  case  for  hydrogen-like  radiators. 
For  non-hydrogenic  radiators,  however,  this  is  no  longer  true.  We  then  find, 
on  inspection  of  Eq.  (D.6),  that  the  quadrupole  term  will  have  an  appreciable 
effect  when 


\AE0AQa.b\  >  2\DaJk\ 


(D.8) 


This  is  an  important  result  that  we  can  use  to  derive  approximate  inequalities 
that  will  tell  us  when  to  expect  the  ion- quadrupole  effect  to  significantly  per- 
turb adjacent  energy  levels  as  compared  to  the  dipole  effect  and  hence  affect 
spectral  line  shapes. 

We  will  make  approximations  to  evaluate  Eq.  (D.8)  for  the  case  of  highly 
ionized  radiators  with  large  values  of  nuclear  charge,  Z.  Consider  an  optical 
electron  in  the  field  of  an  effective  nuclear  charge,  Zeff  =  Z  —  Nb  +  1,  where 
Nf,  is  the  total  number  of  bound  electrons.    Using  hydrogenic  wave  functions 
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for  this  optical  electron,  we  take  for  the  quantum  numbers  in  the  spherical 
representation: 


(D.9) 


a  — >  (n,  l,m)  =  (n,n  —  1,0) 

6-^(n',//,m')  =  (n,n-2,0)  . 
These  satisfy  the  dipole  selection  rules  for  the  z  direction  (i.e.     A/  =  ±1, 

Am  =  0)  and  correspond  to  eigenstates  appropriate  for  connection  to  low 

lying  ground  states  in  resonance-like  line  transitions.    We  use  the  field  from 

a  particle  at  the  average  interparticle  spacing  r,-  as  an  approximation  for  the 

plasma  electric  field: 


Fi  = 


eZ; 


where  r,-  is  approximated  by  the  ion  sphere  radius  given  by: 


47T       Q 


(D.10) 


and  Z{  is  the  charge  of  the  perturbing  ion.  Note  that  the  ion  and  electron 
number  density  are  related  by  ne  =  Z{ni.  For  the  dipole  term  in  Eq.  (D.8)  we 
.61 


use 


where, 


Datb  =  eFi(a\z\b)  =  eFiZ^l 


(D.ll) 


,n,l,m 


■Rnt 


(D.12) 


Zn',/-l,m  -  y  (2/  +  1)(27  -  l)""'''-1 
For  the  case  of  interest  here,  we  have  n'  =  n  and  the  radial  integral  reduces  to 


r: 


»-/-1      2Z.tt 


(D.13) 
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For  the  quadrupole  term  (again,  using  the  field  from  an  ion  at  the  average 
interparticle  spacing)  we  have. 


Qa.a 


e2Zi  /-i  1   2/0^-2, 


,  (a\-r\2coslB-l)\a) 
r,-  J         2 


7">  •     O 


1  /  na0  \2 


(|^)   [5n2  +  1  -  3/(/  +  1)] 


2  1 


2VZ 


c// 


/(/  +  !) -3m 
(2/  +  3)(2/-l) 


Evaluating  these  for  our  particular  case  gives  us 


and 


3  /       2n-l        /47T\2/3e2a0     1/3   2/3 


*  2v^  nW^zVWJ* 


Jeff 


(D.14) 


(D.15) 


|AQa;6|=47rn2(n-l)[ 


4n3  -  12n2  +  4n  +  5 


8n3  -  36n2  +  46n  -  15 


ea0  \z 
7 )    Ue 


KZeffJ 


(D.16) 


The  use  of  these  estimates,  along  with  Eq.  (D.8),  will  give  us  a  condition  for 
the  appearance  of  significant  ion-quadrupole  effects. 
Upon  substitution  we  find  that  the  condition  gives 


(D.17) 


lAP,s8   2  72/3    1/3 
\AEq\  >  -e   Z-'    ne' 

=  3.7xlO-7zf/3ne1/3eV  , 
where  ne  is  in  cm-3.  For  Zj  =  18  and  ne  =  1  x  1024  to  1  x  1025  cm-3  we 
find  that  \AEq\  must  be  larger  than  252  eV  to  547  eV  respectively  for  the  ion- 
quadrupole  effect  to  dominate  the  dipole  effect.  These  energy  splittings  are 
much  too  large  to  be  significant  for  the  natural  splittings  found  in  helium-like 
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and  lithium-like  argon.  However,  for  proton  perturbers  where  Z%  =  1,  we  have 
AJ5q  >  37  eV  for  ne  =  1  x  1024  cm-3.  This  is  more  reasonable,  and  the  ion- 
quadrupole  effect  may  be  important  for  this  case.  Remember  that  this  example 
is  for  levels  with  the  same  principal  quantum  number,  with  opposite  parity  and 
which  satisfy  Eq.  (D.7).  In  conclusion,  we  do  not  expect  the  ion-quadrupole 
effect  to  dominate  the  spectra  of  multielectron  emitters  for  conditions  under 
consideration  in  this  dissertation. 


APPENDIX  E 
FINE  STRUCTURE  CORRECTIONS 

The  fine  structure  correction  to  the  radiator  Hamiltonian  can  be  deter- 
mined from  a  nonrelativistic  expansion  of  the  Dirac  equation  to  order  a2, 
where  a  is  the  fine  structure  constant.87  The  correction  is  given  by 


Hf.s.  = 


P< 


+ 


1  d 


£*>(**) 


(E.l) 


HA  = 


2m2c2 


*Fz 


(Sxp) 


(E.2) 


When  Hfs  is  evaluated  in  the  field-free  total  angular  momentum  basis  \nljjz), 
the  first  three  terms  are  diagonal.  The  fourth  term  can  be  shown  to  be  equiv- 
alent to 


H  H* m  =  W"^  * ' Sa^  x  ^  ' 


(E.3) 


8m3c2       2m2c2 

where  F  is  the  uniform  ion  microfield  at  the  radiator  in  the  z  direction  and 
(p(r)  is  the  radiator  central  potential.  The  electron  spin  angular  momentum 
operator  is  given  by  S,  the  orbital  angular  momentum  operator  is  given  by 
L  and  p  is  the  radiator  electron  momentum  operator.  Shortly,  we  will  be 
interested  in  the  total  angular  momentum  operator  given  by  J  =  L  +  S.  The 
first  three  terms  lead  to  the  usual  fine  structure  energy  shift  term.  The  fourth 
term  is 
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where  a  and  0  are  the  indices  of  the  basis  elements  and  ua  a  is  the  energy 
difference  between  these  two  states.  For  the  field-free  basis,  the  states  with 
the  same  principal  quantum  number  are  degenerate  giving  u>a  p  =  0.  As  a 
consequence,  the  matrix  element  will  be  zero.  However,  this  is  not  in  general 
the  case  for  the  field  dependent  basis  where  this  degeneracy  is  broken  and 
wtt||j  t^  0.  For  the  case  of  the  field  dependent  basis,  H\  can  be  included  in  the 
evaluation  of  Hts.  If  we  do  this  we  find  that  its  matrix  elements  only  appear 
as  off-diagonal  elements  of  Ht  s  .  It  can  be  shown  that 

(a\H4\P)~0(z2a?)  ,  (E.4) 

which  is  to  the  same  order  in  a  as  the  usual  field  independent  fine  structure 
term.  In  this  case,  however,  since  it  only  appears  off- diagonally  in  the  matrix 
representation  of  Ht  s  ,  it  will  only  contribute  corrections  to  the  energy  of 
order  higher  than  Z2a  .  This  can  be  seen  from  looking  at  the  simple  example 
of  a  2  x  2  matrix  eigenvalue  problem  with  only  off-diagonal  corrections  to 
an  already  diagonal  matrix.  Thus,  to  order  a  we  will  ignore  H\  as  a  field 
dependent  addition  to  the  fine  structure  term. 

The  remaining  three  terms  of  the  fine  structure  Hamiltonian  will  give  the 
usual  fine  structure  energy  shift"1'""*'""  of 


En,j  = 


e*2Z4 


n 


(E.5) 


n4     [j  + 1/2      4. 

in  the  total  angular  momentum  basis.  Here,  j  is  the  quantum  number  associ- 
ated with  the  total  angular  momentum  J. 

If  we  ignore  the  field  dependence  of  the  fine  structure  for  a  moment,  we 
can  look  at  the  magnitude  of  the  fine  structure  shift  and  splitting.  For  the 
La  case  of  hydrogenic  argon,  we  find  the  total  transition  energy  AE^  ±  to 
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be  AE^f  =  3319.1  eV  for  j  =  1/2  and  AE^?  =  3323.8  eV  for  ;  =  3/2 
compared  to  the  value  without  fine  structure  of  A-E^i  =  3306.0  eV.  We  see 
a  splitting  of  4.7  eV  and  a  maximum  shift  of  17.8  eV.  For  the  conditions  we 
are  considering  these  are  easily  observable.  We  will  therefore  want  to  examine 
the  influence  of  the  fine  structure  effect  in  our  calculation  of  the  La  line.  For 
the  Lp  line  the  splitting  will  be  reduced  due  to  the  1/n4  dependence  of  the 
fine  structure  term.  Comparing  splitting  of  the  extreme  components  of  the  Lp 
transition;  for  j  =  1/2  we  have  &e\'±  =  3935.1  eV  and  for  j  =  5/2  we  have 
A.E3Y  =  3937.0  eV.  This  gives  a  total  splitting  of  1.9  eV  and  a  maximum 
shift  of  18.5  eV.  Given  the  relatively  broad  Lp  line  at  these  densities,  it  is  safe 
to  ignore  its  fine  structure  splitting.  We  also  neglect  the  overall  shift  of  the  Lp 
line  since  we  are  primarily  concerned  only  with  the  line  shape. 

In  order  to  calculate  the  field  dependent  fine  structure  term  for  the  La 
case,  we  will  transform  the  matrix  elements  of  Hrs  from  the  field  free  basis 
\n,l,j,jz),  in  which  it  is  diagonal,  to  the  field  dependent  parabolic  coordinate 
basis  that  we  are  using  to  calculate  the  other  matrix  elements  in  the  line  shape 
problem.  To  transform  an  element  of  one  basis  to  that  of  another  representa- 
tion, we  introduce  a  unitary  transformation  operator  T  such  that 


T\il>  >=  W  >  , 


(E.6) 


where  \xj)  >  is  in  one  representation  and  \tp'  >  is  in  the  other.  The  corre- 
sponding transformation  of  an  operator  A  from  one  representation  to  another 
is  given  by 


A  =  TTA'  T  , 


(E.7) 
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where  the  matrix  elements  of  T  axe  given  by 

Tij  =<  ftty}  >  ,  (E.8) 

and  T  is  the  transpose  of  T.  In  our  case,  we  will  transform  Hj  s  from 
the  total  angular  momentum  basis,  \nljjz),  to  the  |n/m/ms)  basis.  Next,  we 
transform  to  the  parabolic  basis  |ngm/ms).  Finally,  we  transform  to  the  field 
dependent  parabolic  basis  |ngm/ms(A)).  We  follow  this  somewhat  complicated 
route  because  we  can  compute  analytic  expressions  for  each  step  if  we  use 
perturbation  theory  for  the  field  dependent  basis. 

The  components  of  the  transformation  matrix  for  the  first  two  transforma- 
tions can  be  combined  to  give 


<  nljjz\nqmims  >  =      >J      <  n'iiz|n''m/ms  ><  nl1  m\m's\nqmims  > 

l',rtifvm't 

=  (_i)-/+l/2-i,+(l+m/-«?-")/2v/(2j  +  1)(2/  +  1) 
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(E.9) 


where  the  notation  (•  •  •)  denotes  a  Wigner  3-j  symbol.44'73'89  The  final  trans- 
formation to  the  field  dependent  basis  must  be  worked  out  using  the  perturba- 
tion theory  we  have  explained  in  Chapter  III.  Because  the  method  used  for  the 
calculation  of  these  transformation  matrix  elements  is  the  same  as  the  pertur- 
bation calculation  we  have  presented  in  detail  in  Chapter  III,  we  will  present 
only  the  final  results  here. 

The  remaining  transformation  matrix  elements,  as  a  function  of  the  scaled 
electric  field  A,  are  given  by 


with 


and 


Here 
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<n,q,  mt,  m'8(\)\n,  q,  TO/,  ms  >  =  ^  Cp  \p, 
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'P  "  4n2 


(E.IO) 
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A=mai(-n'2 ,— 2p) 


Additionally,  we  have 


Q)  ~  *m,m' ^m.ni^.n'j  » 
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The  remaining  factors  are 
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with 
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/3    P 

N'p  =  T^<T,  [aiai-i(*)h,k  +  bW-i(y)h,k]   ■  (E.18) 

i=0  h,k 

The  transformation  matrix  elements  given  by  Eqs.    (E.9)  and  (E.10)  are 

what  is  required  to  transform  the  fine  structure  Hamiltonian  from  the  total 

angular  momentum  basis  to  the  field  dependent  parabolic  coordinate  basis. 

This  transformation  is  used  to  add  the  fine  structure  correction  to  the  La  line 

in  our  calculation.   Since  the  n  =  2  levels  of  the  radiator  will  be  only  weakly 

affected  by  the  plasma  microfield  for  the  physical  conditions  we  are  examining, 

perturbation  theory  will  be  sufficient  for  treating  the  field  effects.     For  the 

Lp  line,  however,  the  n  =  3  levels  are  strongly  affected  by  the  microfield  and 

direct  numerical  solution  of  Schrodinger's  equation  is  necessary  but  the  fine 

structure  is  only  a  very  weak  effect  and  can  be  neglected.  Consequently,  only  a 

perturbation  theory  calculation  of  the  fine  structure  corrections  are  necessary 

for  our  work. 

The  complete  representation  transformation  can  now  be  carried  out  using 


(nq'm'lm's(\)\Hu\nq'"m'l"m'l'(\))  = 

£       J2   <  nq,m,lm's(X)\nq"m'l'm^  ><  nq" m'/ m'l\nl> 'j> %  > 

q"m'/m'J  l'j'j'2 
qm,ms      \jjz 

x  (nl'j'j'z\  Hf.s.  \nljjz)  <  nljjz\nqmims  ><  nqmims\nq'"m'l"m'g(\)  >    . 

(E.19) 

This  expansion  gives  us  the  perturbation  theory  calculation  of  the  fine  structure 
correction  to  the  radiator  Hamiltonian  Hr(e)  in  the  field  dependent  parabolic 
coordinate  representation.  As  discussed  in  section  4.1,  we  find  that  this  field 
dependence  has  no  significant  effect  on  the  fine  structure  splittings  of  the  La 
line  of  hydrogenic  argon  at  the  conditions  studied  in  this  dissertation. 


APPENDIX  F 
COMPUTER  CODE  DOCUMENTATION 

In  this  appendix,  we  present  information  on  the  use  of  the  computer  pro- 
grams developed  to  carry  out  the  calculations  discussed  in  this  dissertation. 
We  assume  the  reader  to  be  familiar  with  the  above  material  concerning  the 
theory  of  the  calculation.  The  codes  themselves  are  listed  in  the  University  of 
Florida  Plasma  Physics  Group  Technical  Report.  This  is  available  from  C.  F. 
Hooper,  Jr.  or  R.  Mancini  of  the  Department  of  Physics,  University  of  Florida, 
Gainesville,  FL,  32611. 

Field  Gradient  Term  Calculation 
The  program  EZZAPEX  was  used  to  calculate  the  field  gradient  term 
(Ezz)e  given  by  Eq.   (3.50).   This  code  is  written  in  FORTRAN  for  the  IBM 
3090  mainframe  computer  but  can  easily  be  modified  for  other  systems.  Cur- 
rently, the  program  runs  in  batch  mode  with  a  JCL  file  called  EZZAPEX 
JOB.  When  this  JCL  file  is  submitted  to  the  batch  queue,  four  FORTRAN 
files  must  be  included.  They  are  EZZAPEX,  XMESH,  GOFR  and  RDF.  The 
code  is  written  in  single  precision  and  is  run  with  the  double  precision  option. 
The  input  data  for  a  particular  case  consists  of  five  quantities:   the  electron 
number  density  in  cm"3,  the  temperature  in   eV,  the  adjustable  parameter 
a  for  the  radiator-ion  interaction  from  the  APEX  microfield54,  the  radiator 
charge,  and  the  perturbing  ion  charge.  This  input  is  read  in  free  format  from 
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one  data  line.   As  an  example,  consider  the  case  of  Li-like  argon  in  100%  hy- 
drogen at  ne  =  1  x  1024  cm-3  and  kT  =  1000  eV  with  the  APEX  parameter 
a  =  0.3734.  This  gives  the  input  line: 

1.D24  1000.  .3734  15.  1.      .  (F.l) 

The  details  of  the  numerical  methods  roughly  follow  the  calculation  of  Joyce42 
with  some  modification  as  discussed  in  section  3.1  of  this  work  and  below. 

The  main  subroutine  of  the  code  is  RDF.  It  calculates  the  radial  distribu- 
tion function  g(r)  by  means  of  solution  of  the  HNC  equation.  This  subroutine 
was  developed  by  F.  Rogers.39  The  main  purpose  of  subroutine  GOFR  is  to 
prepare  input  for  and  call  RDF.  The  subroutine  XMESH  sets  the  mesh  for 
the  k  integration  in  Eq.  (3.50).  This  mesh  can  be  refined  if  more  accuracy  is 
required  for  larger  field  values. 

The  output  from  EZZAPEX  is  diverted  to  two  files.  The  main  program 
output  consists  of  the  calculated  values  of  the  field  gradient  as  a  function  of 
the  field  e.  This  output  has  three  parts.  The  first  part  consists  of  a  listing 
of  (Ezz)€/(Ezz)nn  e.  The  second  part  consists  of  the  five  Pade  coefficients 
obtained  from  a  fit  of  the  output  from  the  first  part  to  the  Pade  approximation 
employed  by  Joyce.42  Part  three  consists  of  values  of  (Ezz)e  as  a  function  of  e 
as  generated  from  this  Pade  approximation.  The  Pade  coefficients  are  used  by 
the  field  dependent  atomic  physics  line  shape  code  FDQLINE  to  be  discussed 
below.  A  second  output  file  is  also  generated  for  purposes  of  plotting.  All 
values  of  the  constrained  average  of  the  field  gradient  are  in  units  of  e/r3  z-. 
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Field  Dependent  Line  Shape  Calculation 
The  program  FDQLINE  calculates  the  field  dependent  atomic  physics  line 
shape  function  J(«)  for  the  Lp  line  with  the  higher  order  field  effects  discussed 
in  this  dissertation.    The  different  effects  are  turned  on  and  off  by  means  of 
numerical  switches  in  which  0  means  off  and  1  or  2  means  on.    This  code  is 
also  written  in  FORTRAN  for  the  IBM  3090  but  additionally  uses  the  vector 
processors.  The  JCL  file  for  this  job  is  called  FDQLINE  JOB.  It  must  include 
four  FORTRAN  files,  five  data  files  and  a  line  of  numerical  input  data.   The 
basic  input  data  for  this  code  consists  of  two  data  bases  of  field  dependent 
atomic  physics.    This  consists  of  the  file  STARK-X2  DATA  which  contains 
the  field  dependent  energies,  level  widths  and  dipole  matrix  elements.    The 
data  file  STARQ-X2  contains  the  field  dependent  quadrupole  matrix  elements. 
Additional  data  files  contain  the  field  gradient  Pade  coefficients,  the  APEX  mi- 
crofield  function  and  the  roots  of  the  Laguerre  polynomials  for  Gauss- Laguerre 
quadrature.    The  main  program  is  FDQLINEX  which  calls  all  other  subrou- 
tines. The  subroutines  in  file  FIDATX2  FORTRAN  reads  STARK-X2  DATA, 
and  FIDATQX  FORTRAN  reads  STARQ-X2  DATA.  When  this  data  is  read, 
these  subroutines  process  it  for  use  in  the  line  shape  calculation.   All  atomic 
physics  data  in  files  STARK-X2  and  STARQ-X2  are  in  atomic  units  with  Z  =  1. 
The  subroutine  in  POP4  FORTRAN  generates  the  level  populations  from  ei- 
ther a  Boltzmann  factor  or  the  kinetics  model  discussed  in  section  3.4.    The 
input  parameter  line  consists  of  the  physical  conditions,  the  desired  line  shape 
characteristics  and  the  numerical  switches  discussed  above.    The  first  three 
numbers  are  the  electron  number  density  in  cm-3,  the  temperature  in  eV  and 
the  nuclear  charge  Z.    Next,  comes  an  unused  parameter  set  equal  to  zero. 
The  next  four  numbers  give  the  lower  and  higher  end  of  the  spectral  range  in 
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eV,  the  number  of  Doppler  integration  points  and  the  number  of  data  points 
in  the  final  line  shape.  The  next  two  numbers  give  the  number  of  points  and 
the  step  size  for  the  microfield  integration.  The  remaining  seven  numbers  are 
the  numerical  switches  discussed  above.  We  will  refer  to  them  in  order  and  by 
their  variable  names  as  used  in  the  code.  Switch  NLTE  is  set  equal  to  zero  for 
LTE  populations  and  one  for  NLTE  populations  calculated  from  the  kinetic 
model  discussed  in  the  text.  Switch  IFIDAT  is  set  to  zero  for  field  free  atomic 
physics  and  one  for  field  dependent  atomic  physics.  Switch  NOR  is  set  to  zero 
for  an  unnormalized  line  shape  and  one  for  a  line  shape  area  normalized  to 
one.  Switch  IQD  is  set  to  zero  for  no  quadratic  Stark  effect  approximation  and 
one  for  inclusion  of  the  approximate  quadratic  Stark  effect  energy  corrections. 
Switch  IIQP  is  set  to  zero  for  no  ion-quadrupole  effect,  one  for  inclusion  of 
the  APEX  model  for  the  ion-quadrupole  effect  and  two  for  inclusion  of  the 
nearest  neighbor  model  for  the  ion-quadrupole  effect.  Switch  IGAMA  is  set 
to  zero  for  no  resonance  width  and  one  for  inclusion  of  the  resonance  width 
as  an  imaginary  component  of  the  level  energy.  Switch  IDOP  is  set  to  zero 
for  no  Doppler  convolution  and  one  for  inclusion  of  the  Doppler  convolution 
of  the  line  shape.  The  line  shape  output  is  in  terms  of  intensity  as  a  function 
of  spectrum  energy  in  eV.  Additional  input  data  include  files  with  the  Pade 
coefficients  for  the  field  gradient  term  and  the  plasma  microfield  function.  All 
of  these  files  must  be  included  even  if  the  data  is  not  needed  for  a  particular 
calculation.  The  resulting  spectrum  is  listed  in  the  main  output  file.  Data 
for  the  atomic  physics  data  base  required  as  input  is  generated  as  discussed  in 
sections  3.2  and  3.3. 

The  program  for  the  calculation  of  the  LQ  line  is  FDQLPERA  and  is  very 
similar  to  the  above  with  some  minor  modifications.  There  is  only  one  atomic 
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physics  data  file  needed.  It  is  called  STARPERA  DATA  and  contains  the  per- 
turbation coefficients  for  the  level  energy,  dipole  matrix  elements,  quadrupole 
matrix  elements  and  the  fine  structure  matrix  elements  discussed  in  appendix 
E.  Since  the  electric  field  perturbation  is  relatively  small  for  the  La  cases  stud- 
ied in  this  dissertation,  the  entire  range  of  field  values  used  can  be  treated  by 
perturbation  theory.  The  level  width  is  not  included  so  the  switch  IGAMA  is 
absent.  Since  fine  structure  is  included,  a  new  switch  is  added.  Switch  IFS  is 
set  equal  to  zero  for  no  fine  structure  effect  and  one  for  inclusion  of  fine  struc- 
ture shift  and  splitting.  We  also  specify  the  order  of  the  perturbation  theory 
to  be  used.  Switch  NPMAX  gives  the  order  of  the  perturbation  theory.  Its 
maximum  value  is  six  for  the  data  used  here.  The  main  program  is  contained 
in  the  file  FDQLPER  FORTRAN.  Otherwise  the  code  is  almost  identical  to 
FDQLINE. 

In  carrying  out  the  line  shape  calculation,  it  is  necessary  to  invert  a  complex 
matrix  many  times.  This  is  carried  out  using  the  subroutine  DLINCG  from 
the  IMSL  vector  subroutine  library. 
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